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ef Energies of Cosmic-Ray Particles 
‘ By Cart D, ANDERSON 
4 California Institute of Technology 
a | 
me (Received June 28, 1932) 
ae Cloud chamber photographs of cosmic-ray tracks in a magnetic field up to 17,000 
q gauss are shown. On the assumption that the particles producing the tracks are 
Fe travelling downward through the chamber rather than upward, particles of positive 
i charge appear as well as electrons. From the specific ionization along the track it is 
| concluded that the positives are protons, and are not nuclei of charge greater than 
4 unity. No evidence is uncovered demanding the introduction of a neutron for cosmic- 
et ray phenomena. Eight examples of associated tracks are shown. Energies range from 


below 10° electron-volts to values in a few cases of the order of 10° electron-volts. 





Energy values for 70 tracks are listed. The scattering of cosmic particles in traversing 
a a 6.0 mm lead plate is measured. 


a N AUTOMATIC, vertical Wilson expansion chamber operating in a 
strong magnetic field up to 21,000 gauss was designed for a study of the 
high-energy corpuscles associated with cosmic rays. The expansion chamber 
itself is 15 cm in diameter and has a depth of 2 cm. The axis of the piston lies 
in a horizontal plane in order to effect the most favorable position for photo- 
' graphing the cosmic-ray tracks. The magnetic field was found by direct 
measurement to be homogeneous to within 10 percent throughout the volume 
of the chamber. Photographs are taken through a hole in the pole piece of the 
magnet along the lines of force, thus revealing a particle deflected by the mag- 
netic field as an arc of a circle on the photographic film. A description of the 
g apparatus will follow in a later publication. 
i.) A brief discussion of the results has been published jointly with Professor 
; R. A. Millikan,! to whom the writer is indebted for suggesting the investiga- 
q tion and cooperating in planning it, and whose keen interest throughout the 
course of the work has been of the greatest value. | 
Of the 3000 photographs taken to date, 62 show cosmic-ray tracks of | 
length sufficient for energy measurements, 19 of which are here reproduced 


ws i ee nici Ml as 


t with addition of 3 photographs taken for test purposes. The photographs are 
7/10 full size. The dark back ground in several of the photographs is due to 
; light scattered by the back wall of the expansion chamber. 


? Millikan and Anderson, Phys. Rev. 40, 325 (1932). 
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In the test photograph of Fig. 2 appear tracks in air of the secondary elec- 
trons ejected by gamma-rays from Ra. Figs. 3 and 4 show alpha-particle 
trajectories. 

The remainder of the tracks shown, Figs. 5-23, are ascribed to cosmic rays 
because of their high energy. Cosmic-ray tracks are in all cases readily dis- 
tinguishable from alpha-particle tracks due to the very much greater spe- 
cific ionization of the latter. 


POSITIVELY CHARGED PARTICLES 


A charged particle will be deviated by the magnetic field into an are of a 
circle. The sense of rotation in the chamber as viewed in the photographs 
will be clockwise for a particle of negative charge, and counter-clockwise for 
one of positive charge. The sign of the charge can be ascertained only if the 
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Fig. 1. The number of tracks per unit solid angle as a function of the angle with the ver- 
tical. Only those tracks are included whose curvature in the magnetic field is sufficiently small 
to allow a determination of the direction to be made. Therefore no electrons of energy less than 
100 X 10° volts are included. This space-Jistribution in showing a large percentage of nearly 
vertical tracks differs from that reported by Skolbeltzyn.* 


direction of motion of the particle is known. It is assumed here that the par- 
ticles are traveling downward through the chamber. The small degree of 
scattering to be expected for high-energy particles, combined with the known 
fact that the rays come in from above, appears to justify this assumption. 
In Fig. 1 it is seen that the cosmic-ray tracks prefer the vertical direction over 
the horizontal indicating that backward or large angle scattering is infre- 
quent. 

In many instances on the above assumption as to the direction of motion, 
the tracks are deviated in a sense to indicate the presence of positively 
charged particles as well as electrons. 


* Skobeltzyn, C. R. 194, 118 (1932). 
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A test photograph showing tracks in air of secondary electrons ejected by gamma-rays from Ra, filtered through 


. 12,000 gauss. 
The energies of the electrons range from 10° electron-volts down. 


) 


Fig. 


steel. 
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An a-particle trajectory of 3.85 cm range photographed in air. 


12,000 gauss. 
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A proton of 40 10° volts energy. 


12,000 gauss. 


19, 


ino 


A proton of 450 X 10° volts energy. 


000 gauss. 


/ 


18, 


Fig. 
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SPECIFIC LONIZATION 


A few tracks photographed in an atmosphere of helium with 5 percent air 
show about 14 ion pairs per cm at standard pressure, a value close to that 
found for electrons of about 10° volts energy from Ra gamma-rays, which 
show about 13 ion pairs per em in the same gas. The specific ionization for an 
electron remains practically independent of its energy for energies ranging 
from about 300,000 volts to several hundred million volts. For a given ve- 
locity, protons and electrons ionize the same, and for high energies where the 
velocities of both the electrons and the protons are of the order of the velocity 
of light it becomes impossible to distinguish electrons from protons by their 
ionization. Only at lower energies where the proton velocities are appreciably 
less than the velocity of light will protons show an appreciably greater spe- 
cific ionization than electrons of the same energy. Nuclei of higher atomic 
charge would, for a given velocity, produce many more ions per cm, the 
specific ionization being to a first approximation proportional to the square 
of the charge on the nucleus. 

The specific ionization along the tracks showing positives is in most in- 
stances not much greater than that for the electrons. It is concluded, there- 
for; that the positives can only be protons, and cannot themselves represent 
nuclei of higher atomic number than unity. 

The projection of whole nuclei by the penetrating radiation produced in 
the bombardment of beryllium with alpha-particles has been reported in 
recent experiments. For the explanation of this fact, on the basis of the con- 
servation laws, Chadwick*® postulates a neutron. For the interpretation of the 
cosmic-ray effects so far observed such a neutron is not demanded on the 
basis of the energy-momentum arguments which apply in the experiments of 
Chadwick. Further work will show if the associated tracks of cosmic rays 
represent an effect similar, but on a higher energy scale to the disintegration 
tracks photographed by Feather in the neutron experiments, though the fre- 
quent occurrence of electron tracks in the cosmic-ray experiments seems to 
indicate a different type of phenomenon. 

A possibility to be borne in mind is that, in rare cases, the tracks of curva- 
ture that indicate positives might be in reality electrons scattered backwards 
by the material underneath the chamber and are traversing it from bottom to 
top. Precise data on the specific ionization of the low-energy positives will 
distinguish, however, between downward positives and upward negatives. 


AssociATED TRACKS 
A well-known characteristic of the cosmic tracks is their tendency to oc- 


cur in groups..®? Of the 55 photographs showing cosmic tracks, 7 show 
double tracks and 1 shows three tracks. 


3’ Chadwick, Nature 129, 312 (1932). 

' Feather, Proc. Roy. Soc. A136, 709 (1932). 
5 Skobelzyn, Zeits. f. Physik 54, 686 (1929). 
6 Auger and Skobelzyn, C. R. 189, 55 (1929). 
7 Locher, Phys. Rev. 39, 883 (1932). 
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In general, for paired tracks, the energy of one of the associated pair is 
considerably less than that of the other, in some instances 10° volts and less. 
One of the associated pair is also in all cases definitely an electron. 

The associated tracks have been assumed to be due to the simultaneous 
ejection by a photon of two particles from an atomic nucleus.! 

Another effect which may give rise to associated tracks is a close en- 
counter between a cosmic particle and an electron. Fig. 20 is an example of 
an encounter of this type, the encounter taking place in the wall. For such 
an encounter where an electron of high energy (energy>mc*) produces a 
secondary track, giving to the secondary electron an energy E, the angle @ 
between the primary and secondary electron is given by tan @=(2mc?/E)"?. 
The two tracks of Fig. 10 might represent an effect of this type, the relation 
above being satisfied within experimental uncertainty. The possibility also 
exists that a proton may by direct impact give to an electron energy sufficient 
to produce a secondary track. It is poirted out, however, that on the basis 
of the conservation laws, due to the relatively large proper mass of a proton, 
it is difficult to explain the associated tracks as a binary collision between a 
proton and an electron because of the prohibitively high energy which would 
have to be assigned to the proton in this case to account for the electron 
energies observed. To produce an electron of 100 X10® volts would require a 
proton energy of 10! volts. In both Figs. 8 and 11 one of the associated pair 
has a curvature to indicate a positive, and therefore, if these curvatures are 
correctly measured these cases cannot represent such encounters. Even if 
these are incorrectly measured and are in reality electrons then for Fig. 11 
the value of 6 calculated for an energy of 27 X10® volts is 11°, which would 
agree with the angle indicated by the photograph within experimental un- 
certainty. But for the case of Fig. 8, the calculated value of 17° and the 
measured value of 25° are in conflict. Again the associated pair of Fig. 5 can 
riot possibly be interpreted as due to a binary collision. The hypothesis of 
simple binary collisions is inadequate to account for all the associated tracks. 
Furthermore, encounters of this type in which a large amount of energy is 
transferred to an electron are not to be expected frequently, and the abund- 
ance of associated tracks coupled with the fact that positives as well as elec- 
trons appear, is strong evidence that the associated tracks represent a quite 
different phenomenon, i.e., the ejection of two particles from a nucleus. 


SCATTERING OF THE COSMIC PARTICLES 


Certain of the tracks, Figs. 11 and 17, show sudden though very small 
deflections identical in appearance with the deflections observed in alpha- 
particle tracks due to nuclear encounters, but which are to be expected in the 
gas only rarely on the basis of the present scattering laws for high energy 
electrons or protons. The deflections in some instances represent scattering 
from the walls of the chamber. An example of large angle scattering from 
a lead surface is shown in Fig. 8. In Fig. 10 another instance of large angle 
scattering is shown, the electron being scattered by the glass wall of the 
chamber. 
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Experiments are now in progress to study the scattering of cosmic par- 
ticles in lead. Figs. 22 and 23 show particles traversing 6.0 mm of lead, the 
angle of scattering being in each case readily measurable. 


ENERGIES 


The energies of the cosmic particles as determined from the radii of curva- 
ture of the tracks range from values below 10° electron-volts to, in a few cases, 
values of the order of 10° electron-volts. The greater part of them, however, 
have energies below 500 X 10° volts. 

Precautions were taken to reduce to a minimum the effects of air move- 
ments in the expansion chamber which tend to distort the tracks. The ener- 
gies of the higher energy tracks, i.e., protons of energies of the order of 
500 X 10° volts and electrons of the order of 10° volts, can be determined only 
roughly due to the small curvature in the magnetic field. 


*- 
TABLE I. Energy distribution. 


Number of Number of Energy range in 
electrons protons electron-volts 

2 0 Below 10° 

4 0 From 10°to 10108 
6 1 Ks 10K10®°to 20105 
4 1 20108 to 30108 
3 3 _ 30X10® to 50108 
1 3 Pe 50 X10® to 100108 
5 5 * 100 X10® to 200X108 
4 2 “ 200X108 to 300X108 
5 2 “ 300 K10® to 400 108 
1 4 “400 X10® to 500X108 
3 3 “$00 10® to 700X108 
; 


l “ 700 X 10® to 1000 x 10° 


Group 1 Proton 130 XK 10° See Fig. 5 
Electron 120 x 10° 
Electron 30 X 108 

Group 2 Probably proton 150 x 10° See Fig. 6 
Probably proton 400 x 10° 

Group 3 Probably proton 450 x 105 See Fig. 11 
Electron 27 X 10° 

Group 4 Proton 20 X 10° See Fig. 7 
Electron 4x 10° 

Group 5 Probably proton ~100 X 10° See Fig. 8 
Electron 11 X 10° 

Group 6 Electron 170 x 10° See Fig. 10 
Electron 11106 

Group 7 Electron 180 x 10° See Fig. 9 
Electron 210° 

Group 8 Electron 410° See Fig. 20 


Electron or proton ~400 X 10° 
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Table I gives the distribution in energy of the electrons and protons. The 
number of protons and the number of electrons in various energy ranges are 
listed. It is to be noted that there may be many more electrons in the energy 
range below 10° volts than those listed. Since this includes the energy region 
of radio-activity, it is impossible to distinguish between electrons from radio- 
active sources and low-energy electrons due to cosmic rays. Therefore only 
those very low-energy electrons which are associated with other cosmic-ray 
tracks, and are definitely to be attributed to cosmic rays are listed in Table I. 

There are in addition 5 tracks whose sign of charge is doubtful due to the 
lack of an appreciable curvature in the magnetic field. For these it is possible 
only to assign a lower limit to the energy, i.e., 450 10° volts on the supposi- 
tion that they are electrons and 100 X 10° volts on the supposition that they 
are protons. There is one track, Fig. 21 to which an energy in excess of 10° 
volts must be assigned whether it is assumed a proton or electron. 

In Table II are given the energies of the 8 groups of associated tracks. 

I wish to thank Mr. Seth H. Neddermeyer for assistance in the measure- 
ment of the photographs. 
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Scattering of X-Rays from Solids 


By ALLEN W. Coven 
Ryerson Laboratory, The University of Chicago 


(Received July 5, 1932) 


The total scattering of the Ka radiation of molybdenum from powdered crystals 
of KCI, NaF, and MgO has been observed with an ionization chamber whose window 
subtended several degrees of angle at the scatterers. The scattering from NaF and i 
MgO was found to be approximately the same as that from neon. The scattering from 
KXCl was found to be approximately the same as that from argon. The intensities of 
the diffracted portion of the scattered rays from the crystals have been calculated 
and found to be distributed along the #*/Z curves, thus agreeing with Jauncey’s re- 
lation between gas and crystal scattering, S,=S.+F*/Z. The sums of the diffracted 
intensities plus diffuse intensities have been found to follow the gas scattering curves, 
in agreement with Compton's formula for the scattering from a monatomic gas, 
S,=(P/Z)+01—f?/Z?)/R®. The method of measurement may be used to determine 
values of the structure factor for the atom at rest. 


I. INTRODUCTION 


H. COMPTON! has shown that the intensity of x-rays scattered from 
¢ a monatomic gas or vapor is given by 


So=fP/Z+ (1 — f°/2*)/R* (1) 





where S, is the scattered intensity relative to the value given by the Thomson 
theory of scattering.” f is an average atomic structure factor for the atom at 
rest, Z is the number of electrons in the atom, and 1/ R® is the Breit-Dirac 
factor. R has the value 1+ (h/mcd)vers ¢, where ¢ is the angle of scattering. 
G. E. M. Jauncey* has developed a theory of the scattering of x-rays from 
solids in which Eq. (1) is included as a factor in the scattered intensity. In the 
special case of scattering from a simple cubic crystal consisting of atoms of 
one kind Jauncey’s theory gives* 


f? _— Ff? 1 _— | i 3? F? 


p 
<2) 
tw 
ae 





where F is the atomic structure factor at the temperature of the crystal. Eq. 
(2) connects the diffuse scattering from a crystal with the scattering from a 
monatomic gas of the same kind of atoms as those in the crystal. The term 
F°/Z represents a correction for the effect of interference between atoms in 
the crystal. 

The energy per second emitted in a diffraction line by one of the particles 
in an aggregation of powdered crystals is given by® 









! A. H. Compton, Phys. Rev. 35, 925 (1930). 
2 J. J. Thomson, Conduction of Electricity through Gases, 2d Ed., p. 326. 
°G. E. M. Jauncey, Phys. Rev. 37, 1193 (1931). 

4G. E. M. Jauncey and G. G. Harvey, Phys. Rev. 37, 1203 (1931). 

5 A. H. Compton, X-Rays and Electrons, p. 130. 









422 











SCATTERING OF X-RAYS FROM SOLIDS 423 


, _ 1 rspFet(1 + cos? 26)dV 


8(sin 0)m?c4 





, (3) 


where @ is the glancing angle of incidence on a crystal plane, p is the number 
of planes, F’ is the structure factor for the crystal, is the number of atoms 
per cm’, J is the intensity of the primary beam, and dV is the volume of the 
particle. If the diffracted rays are measured at a distance r with an ionization 
chamber whose window has a length L which is short compared to r sin 8@, 
then the energy per second entering the ionization chamber is 


P, = PL/2zxr sin 20. 


Thomson's classical theory* gives the intensity scattered from the nZdV 
electrons in one of the powdered crystal particles to be 

TnZe*(1 + cos? d)dV 

A => SS Se ae = ’ (4) 


2m*ctr? 





and therefore the energy entering the ionization chamber is equal to A/,, 
where A is the area of the chamber window. To get the energy per second in 
a diffraction line in the units corresponding to S, and S, of Eqs. (1) and (2) 
we take the ratio of P, to AJ,. We then have 


P, w3 pLrF”? 
Al, 


82ZA sin 9 sin 26 


But the atomic structure factor F in Eq. (2) has a value F’/2 for a simple 
cubic crystal consisting of atoms of one kind. Substituting 2F for F’ in the 
above expression, we then have for the energy in a powdered crystal line 
where the crystals are of the simple cubic type and contain similar atoms 


P, nrd* pLr F? ” 
pe ee  _ ine (5) 
Al, 2rA sin@sin 26 Z 

If the ionization chamber window is of such width as to subtend several 
degrees of angle at the scatterer, there will be several diffraction lines entering 
the chamber at large angles of scattering where the lines are close together, 
and probably more than one line at each of the smaller scattering angles if 
the scattering crystals have a grating constant that is not too small. The 
total diffracted energy entering the chamber is then equal to the sum of the 
energies in the lines and may be expressed as 


nd*pLr F -) 


S4' = (~ ——-— (6) 
2rA sin@sin 26 Z 


The total scattered intensity entering the ionization chamber is the sum 
of the diffuse and the diffracted intensities. Using previously measured values 
of F, Eq. (6) may be used to calculate the diffracted portion of the scattering 
and this portion may then be added to measured values of the diffuse scatter- 
ing to get values of the total scattered intensitites. Values of S,’ are in terms 
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of energy so the corresponding quantity representing the diffuse scattering 
will be represented by the area under the curve of diffuse intensity between 
the limits of the (sin ¢/2)/A ordinate corresponding to the angle subtended 
by the ionization chamber window. Hence we may divide each value of Sv’ by 
a value dx, which is the length of the ordinate corresponding to the range of 
scattering angle subtended by the ionization chamber window and add the 
quotient to the intensity of diffuse scattering at the value of the ordinate 
corresponding to the position of the middle of the window. The tetal scat- 
tered intensity is thus represented by 


S, = S.+ Sq'/dx = S.+Sa (7) 


where Sq is equivalent to S,'/dx. If the total scattering from a powdered 
crystal aggregation is measured with an ionization chamber whose window 
subtends several degrees of angle at the scatterer, a distribution of intensity 
equivalent to that given by Eq. (1) is observed. 


ul] 





Fig. 1. Arrangement of apparatus. 


The author has observed the scattered intensities from powdered crystals 
of KCl, NaF, and MgO, using an ionization chamber whose window sub- 
tended approximately 5 degrees of angle at the scattering plates. The ob- 
served intensities from KCI are in close agreement with the values obtained 
by Wollan® for argon and the intensities from NaF and MgO agree with 
Wollan’s values for neon.’ Wollan’s values are in agreement with values cal- 
culated from Eq. (1), using structure factors calculated from the Hartree 
charge distribution. 


Il. APPARATUS AND PROCEDURE 


Fig. 1 shows the arrangement of the apparatus. A molybdenum target 
x-ray tube was immersed in oil in a lead box placed close to the axis of an 
ionization spectrometer. The scatterers were mounted directly over the axis 
on the spectrometer table. The x-ray tube was operated at a potential of ap- 
proximately 40 k.v. and a current of 35 m.a. A Compton electrometer indi- 
cated the ionization currents. A Soller slit, S, collimated the primary beam. 
The use of this slit limited the divergence of the primary rays to less than a 


6 E. O. Wollan, Phys. Rev. 37, 862 (1931). 
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half-degree of angle. The scattered rays passed through a bronze tube 2.7 cm 
in diameter, 14 cm long, with 1.5 mm walls. This tube prevented stray radia- 
tion from entering the ionization chamber with the scattered rays. The 
window of the chamber was 1.67 cm high and 1.64 cm wide. The distance 
from the window to the axis of the spectrometer was 18.05 cm. The chamber 
contained methyl! bromide vapor at atmospheric pressure. 

Balanced filters of ZrO, and SrO were placed at 7 in the primary beam. 
These separated the Ka lines from the general radiation as shown in Fig. 2, 
where the difference in the transmission by the ZrO, and the SrO is re- 
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Fig. 2. Transmission through balanced filters. 


presented by the shaded portion. This is the portion that was used as the 
primary radiation in all of the measurements of scattered intensities. 

The scattering from the crystalline powders and from paraffin was ob- 
served throughout an angular range of from 20° to 100°. The scattering from 
each of the powdered crystal plates was compared with the scattering from a 
similar plate of paraffin. Since the scattering from paraffin can be calculated 
to a close approximation, the values of the scattered intensities for each of 
the other materials could be put on an absolute scale. For the Ka radiation 
of molybdenum Woo’ has found that the ratio of the incoherent to coherent 
radiation scattered from paraffin at 90° is 4.69. If it is assumed that for 


7 Y. H. Woo, Phys. Rev. 27, 119 (1926). 
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parathn at 90° the effect of constructive interference is negligible, the total 
scattering can be obtained from classical theory, if the incoherent portion is 
corrected by the Breit-Dirac factor. Thus, at 90°, the scattering from paraffin 
is given by 


S, = 0.825(R-*) + 0.175 = 0.93. 


From this value and the experimental measurements of relative intensities 
at various angles, S, for paraffin was found to be equal to unity in the range 
from 30° to 70°. The comparisons of the powdered crystal intensities with 
the paraffin intensity were all made at 60°. 

The pure, finely ground, powdered crystals were pressed into brass rings 
which fitted into holders that could be conveniently mounted on the spec- 
trometer table. The paraffin was melted into a similar ring so that it could 
be quickly substituted for any one of the powdered crystal plates. 


Potassium chloride 


The scattering from a thin plate of the KCI powder 2.5 cm in diameter and 
having a mass per unit area of 0.141 g/cm? was compared with the scattering 
from the paraffin which was also 2.5 cm in diameter and had a mass per unit 
area of 0.275 g/cm?*. The arrangement of the two was as shown in Fig. 1, 
where one scatterer is in the primary beam at A and the other in the primary 
beam at B, each plate being set so that the normal to its surface made an 
angle of 30° with the primary beam. The ionization chamber was set at the 
scattering angle of 60°. The intensity with the KCI at A and the paraffin at B 
was first observed. Then the intensity with the paraffin at A and the KCI 
at B was observed. The ratio of the first intensity to the second multiplied 
into the ratio of the masses per unit area and the ratio of the values of Z/ W 
for the samples then gave the value of S, for KCl at 60°. Z is atomic number 
and W is atomic weight. The paraffin was assumed to be equivalent to a 
molecule of Co;Hs0, for which Z/W has a value of 0.5. No correction for 
absorption of the x-rays in the samples was necessary since by this method of 
comparison the portion of the primary rays transmitted by the scatterer at 
A does not change when the scatterers are interchanged in their positions at 
A and B. The intensity of the scattered rays entering the ionization chamber 
from each substance was of the same order of magnitude by this method, even 
though the absorption in the KCI was much greater than that in the paraffin. 

Relative values of the scattering from KCI were observed by reflections 
from a plate of the powdered crystals 4.5 cm in diameter and thick enough to 
absorb completely the Ka-rays of molybdenum. The scatterer was always set 
so that the face of the crystal plate was at an angle equal to half the scattering 
angle. Thus the absorption of the primary and scattered rays was approxi- 
mately the same in the thickness of the plate penetrated, and proportional 
to the reciprocal of the absorption coefficient. For the relative values, there- 
fore, no absorption correction was necessary. The relative values of S,; were 
given by the relative intensities multiplied into the ratios of the values of 
the polarization factor, 1+cos*@. 
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Sodium fluoride 


The intensities from paraffin and the NaF plate, which had a mass per 
unit area of 0.470 g/cm?*, were of the same order of magnitude and the_ 
comparison between them at 60° was made without resorting to the method 
used for KCI of plecing both-samples in the primary beam. When a scatterer 
is at A and none at B (Fig. 1), S; is proportional to (W/mZ)(1,/is,)(1+ cos? 
)~', where Z is atomic number, W is atomic weight, m is mass per unit area, 
/, is the scattered intensity, and 7, is the portion of the primary rays pene- 
trating the scatterer. 7, is equal to Je~™/*°s‘¢/) where J is the incident ray 
intensity and w/p is the mass absorption coefficient of the incident rays in 
the scatterer. The absorption of the Mo Ka rays in the paraffin and in the 
NaF plate was measured with the Compton double crystal x-ray spectrom- 
eter.’ u/p for paraffin was found to have a value of 0.52 and for NaF 2.63. 
Relative intensities for NaF at various scattering angles were measured 
from the same scatterer that was compared with paraffin. 


Magnesium oxide 


The MgO plate that was compared with paraffin at 60° had a mass per 
unit area of 0.331 g/cm.* The intensity comparison was made exactly as it 
was for NaF and paraffin. u/p for MgO was calculated to be 3.18 from the 
values of the atomic absorption coefficients for Mg and O given on page 182 
of Compton's X-Rays and Electrons. The relative intensities were measured 
from reflections from a thick plate of the powdered MgO crystals as was done 
for KCI. 

III. EXPERIMENTAL RESULTS 


Froman® and Harvey'® have respectively measured the intensities of the 
lines from powdered KCI and the diffuse scattering from a single crystal of 
KCl. Froman’s values of F for KCl have been used to calculate values of 
S,for values of ¢ between 15° and 90°, each value extending over a range 
of 5° and differing from adjacent values by 5°. Thus for each value of @ the 
energy in all the diffraction lines included in the 5° range made up the value 
of Sy. The sum of the S, values would represent the sum of the energy in 
the diffraction lines between the minimum and maximum values of ¢, 15° 
and 90°. At each value of @ the value of S, was added to the value of S, 
as shown in Fig. 3. S, is the theoretical curve for argon as calculated by 
Hartree by his method of self-consistent fields. The data for the curve were 
taken from a paper by Herzog.'' His values for the incoherent radiation have 
been reduced by the Breit-Dirac factor before being added to the values of 
the coherent radiation to get the values of S, from which the curve was drawn. 
The F?/Z curve was drawn from Froman’s data. The circles are values of S, 
and the crosses are the values of S.+S,. The author’s experimental values of 
S,, as given in column 3 of Table I, are plotted as dots. 


8 A. H. Compton, Rev. Sci. Inst. 2, 365 (1931). 
® D. K. Froman, Phys. Rev. 36, 1330 (1930). 
0G. G. Harvey, Phys. Rev. 38, 593 (1931). 
1G. Herzog, Zeits. f. Physik 70, 583 (1931). 
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Jauncey and Williams'? have recently measured the diffuse scattering 
from a single crystal of NaF. Their values were kindly sent to the author by 
Professor G. E. M. Jauncey and have been used together with the F values 
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Fig. 3. Scattering from KCl. Dots S; experimental, circles Sy, crosses S.+5S,4, S, from argon. 
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Fig. 4. Scattering from NaF. Dots S; experimental, circles Sa, crosses S.+S., S, from neon. 


of Havighurst" to calculate values of S.+.S, as shown by the crosses in Fig. 
4. The F*/Z curve is also from the data of Havighurst. The S, curve for neon 
is again from Herzog’s values.'' The author’s experimental values of S, for 


2G, E. M. Jauncey and P. S. Williams, Phys. Rev. 41, 127 (1932). 
18 R. J. Havighurst, Phys. Rev. 29, 1 (1927). 
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NaF are shown by the dots. These values appear in column 4 of Table I. 
Values of S,; are shown by the circles. 

Fig. 5 shows the results for MgO. The author is grateful to Dr. Wollan 
for the values of diffuse scattering from a single crystal of MgO which he has 
recently measured but not yet published. These have been used to draw the 
S. curve. Values of Saas shown by the circles were calculated from Wollan’s™ 
values of F for powdered MgO crystals. The F?/Z curve is also drawn from 
Wollan’s values. Values of S.+.S,4 are shown by the crosses and the author’s 
experimental values of S; by the dots. These are the values appearing in 
column 5 of Table I. The S, curve for neon is the same as the one in Fig. 4. 
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05 
(SIN G/2)/A 
Fig. 5. Scattering from MgO. Dots S; experimental, circles Sz, crosses S.+Sa, S, from neon. 


The values of the total scattering, S,, observed by the author for pow- 
dered crystals of KCl, NaF, and MgO and for paraffin are all shown in Fig. 6 
together with the theoretical curves for neon and argon. The crosses indicate 
the values for KCl, the circles are the NaF values, the large dots are for MgO, 
and the small dots for paraffin. 





TABLE I. Values of S; for paraffin and powdered crystals. 











Angle (sin ¢/2)/ KCI NaF MgO Paraffin 
20° 0.245 7.25 4.16 1.22 
30 0.365 4.77 2.59 3.13 1.00 
35 0.423 1.96 
40 0.482 3.83 2.15 1.30 1.00 
45 0.539 1.78 1.95 
50 0.596 3.14 1.13 1.55 1.00 
60 0.705 2.65 1.28 1.31 1.00 
70 0.808 2.28 1.20 1.10 1.00 
80 0.906 2.01 1.10 1.09 0.93 
90 0.996 1.85 1.00 1.06 0.93 
100 1.082 1 


.80 0.89 








4 E, O. Wollan, Phys. Rev. 35, 1019 (1930). 








430 ALLEN W. COVEN 


IV. Discussion 

Fig. 6 shows that the measured values of the total scattering from KCI are 
in excellent agreement with the theoretical argon curve. Wollan’s experi- 
mental values of the scattering from argon are slightly higher than the 
theoretical curve at values of (sin ¢/ 2)/X below 0.7 and above 0.9. Hence the 
agreement between the author's values of the scattering from KCI and the 
scattering from argon is better for the theoretical values than for the experi- 
mental values. The experimental values of the total scattering from NaF and 
Mg@ are in good agreement with the theoretical neon curve for the higher 
values of (sin @/2)/X. At the lower values of (sin @/2)/X the points do not lie 

















0 5 1.0 
sin(/z)/2 
Fig. 6. Values of S;. Crosses KCI, circles NaF, large dots MgO, small dots paraffin. Curve | 
neon (theoretical), curve II argon (theoretical). 


on a smooth curve. This is to be expected, because crystals of NaF and MgO 
have a small grating constant. The diffraction lines are consequently widely 
separated and a good average of the diffracted energy is not obtained. How- 
ever, the average distribution of the points appears to lie along the gas scat- 
tering curve. The larger grating constant of KCl insures a good average of the 
diffracted portion of the energy throughout the whole range of scattering 
angles at which observations were made. This explains the excellent agree- 
ment of the KCI values with the gas curve at all points. 

Jauncey and Williams’ have found that the values of S.+ F*/Z for NaF 
are lower than the values of the scattered intensities from neon. The author’s 
calculated values of S.+5S., as shown by the crosses in Fig. 4, are also lower 
than the gas scattering curve. The author’s experimental values of the total 
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scattering from NaF are on the average higher than the values of S,+.S.,, but 
slightly less than Wollan’s neon values. The values of S.+5S, are probably 
more reliable than the experimental values of the total scattering because of 
the difficulty of averaging the diffracted energy in the ionization chamber. 
The values of S.+S. for MgO lie slightly above the theoretical curve for neon 
and are in good agreement with Wollan’s experimental measurements of the 
neon intensities. The S.+S, values for KCI are slightly lower than the argon 
theoretical curve, and the author’s values of the total scattering lie between 
the S.+.S, values and the argon curve. The difference between the values is, 
however, less than the probable uncertainty in the determination of the ab- 
solute value of the total scattering by comparison with paraffin. Using the 
F values of Froman® for powdered KCl, the values of S.+ F?/Z fall exactly on 
the theoretical argon curve. James and Brindley“ have measured the F 
values for a single crystal of KCl and have found them to be slightly greater 
than Froman’s values. Using these F values of James and Brindley, the S, 
+ F*,Z values for KCl agree better with Wollan’s values of the argon inten- 
sities than with the theoretical curve. Jauncey and Harvey" have used. 
Wollan’s and James and Brindley’s values in testing the S.+ F?/Z relation. 

In a former paper" the author has published the results of measurements 
of the total scattering from aluminum, copper, and lead sheets. Since a sheet 
of metal consists of a collection of small crystals, we may consider it to be a 
plate of powdered crystals of the metal. The total scattering from a metal 
sheet measured with a wide-window ionization chamber should thus be ap- 
proximately the same as the scattering from the same metal in the vapor 
state, for values of (sin @/2)/A greater than 0.5. This method of measurement 
vields values of the atomic structure factor for the atom at rest, since from 
Eq. (1) we have 


The author wishes to thank Professor A. H. Compton for suggesting 
measurements of the total scattering from powdered crystals and Dr. E. O. 
Wollan for his generous and helpful advice during the course of the experi- 
ment. 


 R. W. James and G. W. Brindley, Proc. Roy. Soc. Al21, 155 (1928). 
1 G. E. M. Jauncey and G. G. Harvey, Phys. Rev. 38, 1071 (1931). 
17 A. W. Coven, Phys. Rev. 38, 1424 (1931). 
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Sum Rules for Atomic Transition Probabilities 


By J. P. VINtTI 
Massachusetts Institute of Technology 
(Received June 27, 1932) 

Some new sum rules for atomic transition probabilities are developed; they are 
similar to the well-known rule of Kuhn, Reiche, and Thomas. They give, if unpolarized 
light is incident on any gaseous element in any given discrete energy level, the total 
absorption, excitation plus ionization, with various types of frequency weighting. 
The results are similar to those of Dirac and Harding, but are more general in that 
they apply to any discrete state of any atom, rather than to the ground state of hydro- 
gen only. The method differs from that of Dirac and Harding. 


INTRODUCTION 


HERE are well-known relations connecting the coefficients of spon- 

taneous and induced transitions in atomic systems (the Einstein A’s and 
B's) with the matrix components of electric moment between the initial and 
final states and the corresponding frequency differences. A»,, the spontaneous 
transition probability from a state 7 to a lower state 0 is equal (in the dipole 
approximation) to a universal constant multiplied by the cube of the fre- 
quency difference and the square of the absolute value of the matrix com- 


ponent of electric moment, e7() Xon) 7+ | You) 7+) Son) 7), Where ex, ey, es, are 
the rectangular components of electric moment, so that e.g., x is equal to the 
sum of the x coordinates of all the electrons. Bo,, the induced transition 


“4° ° . . o | Il» | lo | lo ° 
probability, is proportional directly to e7() Xon)?+ | You) ?7+ | Son,*), Where in 


this case, either 0 or 7 may be the lower state. 

Xon is given by an integral of the form /¥o(Z,-1"x,)Yodr, where Yo and y, 
are the wave functions of the atom for the states 0 and 7, N is the number of 
extranuclear electrons in the atom, and the integration is extended over the 
coordinates of all the electrons. We see, then, that in general one must know 
the wave functions of the atom for the states concerned in order to calculate 
A», and Bo,. One, can, however, give some general relations satisfied by the 
A's or B’s without knowing any wave functions, or knowing at the most the 
wave function of one discrete state only. 

The sum rule of Kuhn, Reiche, and Thomas! is an example of the type of 
thing one can say about the matrix components, and thus about the A’s and 
B’s, without knowing any wave functions. Using atomic units, and letting 
x denote »;"x,, the sum of the N cartesian coordinates, W, the energy of the 
state n with reversed sign, this rule can be expressed in the form Y,(Wo— W,) 
|xon |2= NV, where the state 0 denotes an arbitrary state of the discrete spec- 
trum, and the summation is understood to include a summation over all the 
discrete states that have dipole combinations with the state 0, and an integral 


1 W. Kuhn, Zeits. f. Physik 33, 408 (1925). F. Reiche and W. Thomas, Zeits. f. Physik 34, 
510 (1925). 
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over the continuous spectrum. (Similar relations hold for the y and z com- 
ponents of electric moment.) This rule is directly deducible from the commu- 
tation rules and the relation connecting the matrix components of momentum 
and coordinate. In terms of the A’s it enables us to evaluate at once *&,, 
Aon/Von?, Where vo, is the frequency of the spectral line associated with the 
transition from a running state n to the base state 0. Furthermore, Wigner? 
has shown that this rule may be sharpened by dividing transitions into three 
classes: +, —, and 0, + transitions being those in which the total azimuthal 
quantum number L increases by 1—those in which it diminishes by 1, and 0 
those in which it does not change. Wigner’s rule then gives the sum g,(L)=* 
+¢)(L)=° in terms of the L-value of the state 0, where g, and gy are simple 
functions of L, and the summands are the f values Ao,/¥o,?. We wish to men- 
tion at this point that all the sum rules developed in this paper can be 
sharpened in a similar way. 

Matrix multiplication of the simple type xx = x? furnishes another general 
relation connecting the matrix components: 


| | 9 
> 2 Xon | * = (X*)oo. 


n 


Here, as above, 0 means a discrete state and the summation is to be extended 
by means of integration over the continuous spectrum. This relation can also 
be looked on as a sum rule, although of a less general type, since the evalua- 
tion of the sum depends on a knowledge of the wave function of the state 0. 
It gives directly Y,Bon, if one knows Yo. 

The above rules tell us how to evaluate S,,(Wo—W,)’, xon|? when j has 
the value 1 or 0. We develop in this paper other rules of the same general type, 
letting 7 take on the values 2, 3 and 4. They constitute rules for the summa- 
tion of the A's and B’s (of a given base state) over all the states of the atom, 
with various kinds of frequency weighting. 


SECTION I 


We begin with the case j=2, and attempt to evaluate >>,.(Wo—W,)? 
(| xon| 2+ | yon! 2+ | Zo 2). (We abbreviate expressions of this type hereafter 
to Don.zy:(Wo— W,.)?| Xon|?). Let xy, Yu, 3, denote the coordinates of the uth 
electron, N the total number of extranuclear electrons in the atom, pi, Pyu, 
p:, the components of momentum of the uth electron, and p,,(0n) the matrix 
component of p,, between the stationary states 0 and n. Also, we shall use the 
contractions: 

N 


N 
> x, = x, Dope = pz, and Wo — Wr = Won. 
1 pol 


= 


Now 


_ @ 
Pxu(On) = (h/ 2ri) fas -W,dr, 
OX, 


2 E. Wigner, Phys. Zeits. 11, 450 (1931). 
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and it is well known that when at least one of the states 0 and » is a discrete 
state (to insure the vanishing of surface integrals at infinity which occur in 
the derivation), this can be thrown into the form m(277/h) Wo,x,(0n), where 
x,(On) is the matrix component of x, between the states 0 and m, and m is the 
rest mass of the electron. On summing over the index » from 1 to N, we 
obtain the relation: 


pAOn) = m(2ri/h)WorzXxon. (1) 

Similarly 
p.(n0) = p(On) = m(27i/h)W ox, (2) 
(since p, is hermitian, and Wo = — Wo,). On multiplying (1) by (2) we obtain: 
p(On)p.(n0) = m2(2r/h)2W 5,7) Kon \?. (3) 


Summing (3) over all states m (complete set), and making use of the simple 
rule for a matrix product, v/s: 


> p.(On) p.(n0) = (pz)oo, 


we obtain the relation: 


m?(2x/h)? >> Won? | xon |? = (p2?)o0- (4) 
Similar relations hold on replacing x by y or s. Adding, 
m(2x/h)? So Won? | xon| 2 = (p22 + py? + peo. (5) 
n,ryz 


We next attempt to give the right hand side of (5) a physical meaning. 
Since 


a 
p2= 2 Dpsb:, = Loz? + DL Ldz,dz, 
im v w=1 uv 

N 
(pz")oo = D (Pz,*)00 + > > (PzuP 2) 00. 
w=1 pw 


The right side of (5) thus decomposes into two parts: 


N 


Do (Pau + Py? + peu®doo and D> do (Prue + PyuPyw + PeuPs)oo- 


w= 1 u*y 


The former term is recognizable as 2m7 oo, where J is the diagonal compo- 
nent of kinetic energy in the state 0. By the virial theorem, making use of the 
fact that the interaction potentials in the atom are all Coulombian, we can 
equate 7) to — Wo. The cross-product term is expressible as 


—(h 2n)* f ¥ Vu Vaodr, 
or 


Z > (Pu P»)o0; where Pu - (P zu; Pw, pu). 


uy 
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Then 
m?(2x/h)? > Won? | Xon | 2 = — 2mWy + ) a > (Pu Pr) 00 (6) 
uy 


n.zys 
where | W,| is the energy necessary to strip the atom of all its electrons when 
in the state 0. 

Eq. (6), then, is the general sum rule for the case j= 2. If the wave func- 
tion of the state 0 is simple enough so that >> >-,.,(p,-B,)oo can be evaluated, 
it is of use. For states like the ground state of helium, where one can assume 
separability to a good approximation, 


(P suP x) 00 = P uP x» 


where the tilde indicates an average with respect to the wave function of that 
electron alone. But 


Pm _ Pw =(). 
In such a case, and for one-electron atoms, the rule becomes 


> Won? | Yon | 3 = (h? 232m) Wo. (7) 


n,@ryz 


For those cases above mentioned in which the cross term vanishes, it is 
of interest to express the rule in terms of the A and B coefficients: 





Aon 3277¢? 4 ; 
> = gE my | HW 0 | (8) 
- ie sme h? 
or 
4re? 
ini , 
y Von Bon = Pegi i 0 (9) 
n smh? 


the Kuhn-Reiche rule (j= 1) being expressible as 


Aon S8z7e? 
oe ~ = — 
n VOn mie 
or 
9 
Te - 
> von Bo», == ag A e 
n om 


Alternative method 


The sum >>, Wo,?|xon|? can be expressed alternatively as }>.)~,.Won?x, 
(On)>_,x,(n0). (In the rest of the paper matrix components of coordinates and 
momenta are understood to have time factors of the form exp|(27i/h) Wont |.) 
As before, p:, (On) =m(27i/h) Wonx,(0n), and making use of the above con- 
vention as to the time factor of a matrix component, we have, letting the dot 
denote differentiation with respect to time: 


Pxu(On) = m(2ri/h)*Wo,2x,(On) 
so that 


HW o,.2x,.(On) 


(1/m)(h/2ri)?p.,(On). 
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The sum then becomes (1,/m)(i/ 277)? >. >, So ypeu(On)x,(n0). We now make 
use of the fact that p,,(0n) = —(OV, Oxy)on, Where V is the potential energy 
of the atom,’ and of the rule for a matrix product, obtaining for the sum the 


mae) EEC, *) 


Summing on x, vy, 5, we get 


1/hk\ 
> Won?) Xo. 12 = (~) >, »( Val’: Rodeo 


Mi 


expression 


where \, is the gradient operator in the coordinates of the uth electron, and 
R, is the radius vector to the vth electron. This expression can be split into the 


1/h\? & 1/ hy? 
— ( ) > ( uV- R,)oo and - ( ) >» > ul- Ry )oo. 
m Ir sen'l Ml Jr) Fy 


parts 


=f 


The first part can be transformed by Euler’s theorem (since V is a homo- 
geneous function of the electron coordinates of degree —1), and becomes 
(—1/m) (h/27)?Voo. Application of the virial theorem, according to which 
Voo =2Wo, gives finally for the first part — (A? 272m) Wo, agreeing with the 
non-cross term of the first method. The cross term (1/m7)(i/27)°3> ¥oy., 
(V, V +R) ‘oe must thus be equal to that obtained by the first method, 
(1/m?)(h, 27)? > Do unol Pu P»)oo. It does not appear that this equivalence 
could be shown in any other way. The new form for the cross term is more 
useful for calculation, since it does not require us to differentiate a wave 
function which may be unwieldy. For an atom of atomic number Z with V 
electrons it becomes: 


m lr = Os 


r x d > e (ri ili r* - os ~ ro?) | 


od 
i-6 


For the special case of a two-electron atom the cross term becomes 


1 1 e 
=(- “) fv & rile COS (), {- . + ) + utr 
m 3 ry® ae 


where ©). is the angle between the radius vectors. If this is calculated out with 
a separated wave function, the term e?/7. leads to a nonvanishing result, in 
contrast to the first method, which gave zero. This discrepancy is un- 
doubtedly due to the inaccuracy of the assumption of separability. 


’ The proof is like that given by Sommerfeld, Wellenmechanischer Erginzungsband, pp. 
288, 289 for a diagonal component of momentum. The operators div, grad, and A, however, 
are now operators in 3 N-space, rather than 3-space, and with us the time component of a wave 
function of a stationary state k is exp | —(27i/h)1V,¢] rather than exp [+(2xi/h) W;t], so that 
Zon = fyopndr, rather than Jory »dr, as with Sommerfeld. 
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SECTION II 
y e ° ° Pa . , | o 
We continue with the case 7 = 3. The problem is to evaluate >_,. Vv On’ | Xon| 2. 
First, Wo,°| xon|*, because of the hermitian character of our matrices, can be 
expressed as —(Won?Xon)( Wrox). The factor Wrox,0 can be then placed 
equal, as before, to (i/2zim) p.(n0). Also, as before, 


, 2ri\? 
p.AOn) =m Woon? on 
h 


Thus 
; wart. 
Won?Xon = -) — p.(On) 
2Jrif m 
and 
h he, 
Wo.*| Xo. | ? = — p (On) p.(nO), 
2ri 4x°m? 
Now 
. N . 
p.(On) = >> P xu(On) 
p= 1 
and 
: oV 
P xu(On) = — ¥ 
OX, On 
Thus p,(0n) = — >°,(0V/Ox,)on, so that 
h h? oV 
Won? Yo, 27=— Sa pe = ( ™ -) Pr(nO). 
dri 4r7m? s OX,/ on 


Upon summing on », making use of the rule for a matrix product, we obtain 
h h? av 
ee ee Ce a, ee px) (10) 
. 2ri 4x*m? , OX, 00 


Eq. (10), then, is the rule for 7=3 in the general case of an N-electron 
atom. (The 7, of course, is cancelled on performing the indicated operations. ) 
We next consider the special cases V =1 and N =2. 


Case l: N=1 


Letting 0 denote any discrete state of an H-like atom 


a 1/h\i/oaVv a 
wi ou | Xu) ~ = 1 a 
k m*\ dr Ox IAx/ oo 


1/h\3 ; 
> Wo? | Xor| ? = —(—]/J, where J = fv 7V- VwWodr. 
k,xryz m~ ar 


But V=V (r only). So that 


and 


dV Apo Ze? do 
dr Or r Or 
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. 1 Ay» 
J= ze | Wo dr 
r- or 


and since Po = R(r) (8) B(¢), 


; »* dR ‘sites : 
J = Ze? | R dr } )? sin 0d9 
) dr v7 


(since V1 = —Ze? r). Thus 


en 
} Pddd. 


Letting R, O, and ® be separately normalized, 
; . dR Ze? Ze? 
J= ze f R dr = | R2\ x)— R20) | =— - R*(0), 
0 dr 2 Z 


since the radial wave function always vanishes at infinity. Then, making use 
of the properties of the Laguerre polynomials involved in the radial wave 
function, we find: 

J=0,140 


= — (Je?/a_)Zi/ n®, 1 = O 


where do is the Bohr radius, » the principal quantum number, and / the 
azimuthal quantum number of the state 0. Thus the cube sum vanishes for 


140, and equals 
Lf b\* 2e* Z 
—_ ~— 308 j= QQ, 
m*\2r/ ay*® n° 


It is interesting to express these relations in terms of the A's. 


Since 
O4rte? 
Aor = e I or | Voz ° + Nos : + =0 \ hs 
3h 
Dirlor = 0,140 
Se a" 


= ————— — , J/=(0). 
3mcao> 3 

Expressed in words, in any hydrogen-like atom, the sum of the sponta- 
neous transition probabilities from any state not an S-state to all lower states 
is equal to the sum of the spontaneous transition probabilities from all higher 
states down to the given state. For an S-state they are unequal, the difference 
varying inversely with the cube of the principal quantum number and directly 
with the fourth power of atomic number; for S-states near the ionization 
limit they approach equality. (By a sfa/e is here meant a level with given / 
and m; if one wishes to treat state as meaning a given / only, that is, to 
average over different orbital orientations, one must take into account such 
orbital degeneracy in writing down the formulas for the A’s or the B’s in 
terms of matrix components. The necessary quantum weight factors can be 
found most explicitly in a paper by Y. Sugiura.*) 


* Y. Sugiura, Phil. Mag. [7] 4, 495 (1927). 
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For example, the sum of the spontaneous transition probabilities from all 
higher states, continuous included, to the state 2P (with given m) in hydro- 
gen-like atoms is equal to the spontaneous transition probability from that 
state to the groundstate 1S. (m here means magnetic quantum number.) 

Applied to the case in which 0 denotes the ground state of an H-like atom, 
this cube rule furnishes an independent proof of a theorem familiar from 
direct examination of the hydrogenic functions. We have shown that this 
sum is of the order Z‘/n* for S-states of H-like atoms, and zero for other states. 
Since for the ground state all the Ao,’s have the same sign, it follows that 
their sum cannot vanish, so that we have: the ground state of an H-like atom 
must be an S-state. 


Case 2: \ =2, helium-like atoms 


On writing out the double sum, we obtain: 


a , h h? fav oV oV oV 

Ds Won®| ton |? = — — omy wane Bag a eo Dag a Baap = Pa « 

* 2ri 4r°*m*\ 0x, ONXe Ox, ON 00 
Using V = —Ze?/r,—Ze?/re+e?/ rig, and carrying through the indicated opera- 


tions in a straightforward way, we obtain: 


h\* Ze 1 d 1 d 
> Won? | Xon | Sos (<-) sa E ET . re — Siac 


oe Qe} mtr? gy: OX, 7 on OX 


ri > ayz ONe eee! 


, 1 ‘ 
+ - pT —- Tt - . > ¥2 | ° 
2° 00 


Letting 1,,, le,, 15, be unit vectors along the radius, meridian, and pa- 
rallel of latitude through the wth electron, 


> 0 0 
X,—— =r, — 
ryz OX, Or, 
Also 
0 
DX - rulry: Vv 
syz OX, 


where \, is the gradient in the coordinates of the vth electron. Using these 
abbreviations, we can write the cube sum for N = 2 in the form: 


h\+Ze? 1 @ 1 @ 1 1 
J Stolp ten Va than 0s [dade 
2x} m? r,? On, ro” Org ry? re” 


Knowing the wave function of any discrete state of a helium-like atom, 
we have here a rule for the evaluation of }>,Ao, for the atom. We might men- 
tion that for S-states (e.g., (1s)* or (1s)(2s)), an approximate function of the 
form Yo = Yo(", 72) (no angular dependence) is often used. In such a case 
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0 

1yi- ' Wo = COS @ i Vo 
Ore 
0 

] °° Wo = COs 0 Vo, 
Or; 


where ©. is the angle between the radius vectors r; and rs. But focos Oh. 
sin 6,d0\d¢, sin O2d62d@2=0. Thus if Yo= Wo(n, 72), the terms involving both 
electrons drop out of the integrand and the cube sum reduces to the integral 


h\+Ze ¢ 1 0a 1 @ 
pie fv. 4 _—— Wodr. 
2x} m? ry? Or, ro” Ore 


In applying this to the ground state of helium, it is convenient to express 
the cube sum in atomic units, in which case the factor in front of the integral 
becomes 8Z. Using a simple screening constant wave function Y= — (a?/7)!/? 


r 
4 


e “tm where a=Z—5,16, the cube sum becomes — 32Z(Z —5/16)*. 


SECTION III 7 =4 


To evaluate )>y.2y:Won'| Xo.|?. Taking the general case of an N-electron 


atom, we begin as before: 
(—) 
OX, On 


_ Peu(On) 


II 


m(27/h)*1V'9,2x, (On) . 


Thus 


OV , 
bo - = m(2r/h)* Won? von 
On 


OX, 


a 


oV oV 
m*( lr h)*W ont Von 2 = » ® x(— ) (< ) e 
m= v NOX on \OX/ no 


Summing over all the states n, 


and 


IV aV 
>> (27 h)*m? Wont | Xon| 2 = >» >(— —] 
v v 00 


n oe 
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For a hydrogen-like atom, V =1, we obtain for the fourth-power sum: 
( h )=( 1 ) 
2rJ m?*\r*/ 
3 (i + 1) 
Zi1-— 
an? 


; 
(".). aon + 3/2)(0 + 1)(U + BU — 3) 


where 


This fourth-power rule gives a rule for summing >> ,9,A ox. The presence of / 
in the denominator of the above expression leads to the result that a sum of 
transition probabilities with this sort of frequency weighting diverges for 
S-states of hydrogen-like atoms. 

For a helium-like atom, this rule takes a simple form. Using 


we find 


Summing over x, ¥, 3, 


oV\? 1 1 2 
5 (x)= nts ees) 
. of rit orgt ory? re? 


ryz rm 





where ©. is the angle between the radius vectors rf, f2, from which is obtained 
the fourth-power sum: 


h\*Z2e4 1 1 2 
a wasien fv —-+ + - cos io Wodr. 
2x] m? rit ore*t oye” 


Going back to H-like atoms, we can use the cube and fourth-power rules 
to deduce the asymptotic behavior of the A’s between the ground state and 
the continuous spectrum. >, .7%A ox diverges for S-states, and therefore for the 
ground state. For this state the vo,A 0; all have the same sign. For the discrete 
spectrum (Lyman series), vo,Ao, is asymptotically of the order 1/m,°, and 
since >.:“1/n;3 is known to converge, it follows that the sum > vA ox over 
the discrete spectrum only must converge. Thus /v(dAo,/dv)dv over the con- 
tinuous spectrum must diverge. In a similar manner, using the cube rule, it 
follows that /(dAo,/dv)dvy over the continuous spectrum must converge. 
Writing dAo,/dv asymptotically (large v) as const./v'** we see that /,,*dv/v'** 
must converge and /, “dv/v* must diverge, (where vp =frequency of Lyman 
series limit). The first requirement gives ¢€>0, the second €<1, so that 
0<e<1. Thus the f-value dfo,/dv(=v~? dAo,/dv) is of the order 1/v’, where 
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3 <j <4. This example serves to illustrate how the sum rules can be used to 
obtain information about the intensity distribution in the continuous spec- 
trum without knowing anything about the continuous wave functions. 

It should be mentioned that Dirac and Harding® have already obtained 
the sum rules of this paper for the special case when the state 0 is the ground 
state of hydrogen. Besides this restriction, their method encounters an 
ambiguity of a peculiar sort, in that part of an operator may operate back- 
wards and part forwards; this ambiguity has to be removed by a special 
investigation. 

The author is indebted to Professors J. C. Slater and M.S. Vallarta,who 
have read the manuscript and offered suggestions, and especially to Professor 
P. M. Morse, who has taken a very helpful interest in the paper and offered 
many valuable suggestions. 





® Dirac and Harding, Proc. Camb. Phil. Soc. 28, 209 (1932). 
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Multiplet Splitting and Intensities of Intercombination Lines 
Part I 


By Hucu C. Wo.Fe* 
Utrecht University 


(Received June 13, 1932) 


Kramers’ symbolic representation method for the treatment of those properties 
of free atoms which can be derived on the basis of irreducible representations of ro- 
tations in R; is outlined. Wave functions and operators are represented by functions 
of spinor quantities (£, ») whose properties are known from the theory of invariants. 
Integrations yielding matrix elements are only symbolic and the results, being de- 
termined by the rotation invariance properties, contain therefore undetermined con- 
stant factors of the nature of radial integrals. The method is particularly adapted to 
the problem of multiplets (interaction between two or more “vectors”—e.g., orbital 
angular momenta, spins—treated as a perturbation). The method is applied to con- 
figurations involving two valence electrons of which one is in an s state. Such a con- 
figuration gives rise to a singlet and a triplet with the / value of the second electron. 
The deviations from the normal / to /+1 interval ratio in the triplet are due to a re- 
pulsion bet ween the singlet level and the center triplet level, their 7 values being the 
same. The mixing up of the wave functions of these two levels gives rise to singlet -trip- 
let intercombinations. The constants which represent the interaction energies (1) in- 
terchange, (2) spin-orbit, (3) spin of one electron—orbit of the other (assuming that 
other types of interaction may be neglected) are found in terms of the three intervals 
of the multiplet. Corrections to the Kronig-H6énl intensity formulas are found in 
terms of the intervals of the two multiplets between which the transitions occur. 
These formulas then give the relative intensities of combination and intercombination 
lines. The sum rules for intensity hold for the complete multiplet and not for the 
singlet and triplet separately. Part II will contain application of the theory to ob- 
served spectra. 


N HIS book, Quantenmechanik und Gruppentheorie, \Wey| has shown that 

many of the properties of free atoms can be derived on the basis of irre- 
ducible representations of the rotations in three dimensional space. Kramers! 
has obtained a very elegant formulation of Weyl’s methods through the use 
of symbolic representations of wave functions and operators with properties 
well known from the theory of invariants. We propose here to outline this 
method and to show its application to a problem in complex spectra. 


THE SYMBOLIC METHOD 
Be 
Let £ and 7» be the components of a spinor or half-vector, i.e., a two- 
dimensional complex vector. The transformations of the monomials &'*"n'~" 
(m= —/, —1/+1, ---/), when £ and 7» suffer a unitary transformation, con- 
stitute an irreducible representation of degree 2/+1(/=1/2, 1, 3/2, -- - ) of 
the space rotation group. Except the identity, there are no other irreducible 


* Fellow of the Lorentz Foundation. 
1H. A. Kramers, Proc. Amst. Acad. 33, 953 (1930) and 34, 965 (1931). 
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representations. Then a series of 27+1 quantities G;, which transform into 
each other in an irreducible way under space-rotation, may be gathered up 
in the one formula 


> — 


‘ 


in 
>( Jar (— b)t¥G' = (an — DE)”, (1) 
/ 


where a and } are the components of an arbitrary constant spinor. The G; 
can always be so chosen that they transform like &**/y’ 
lant. 

2. 


f 


!. (an—bé&) is invar- 


In the study of multiplets, a free atom may be treated as if it consisted 
of two or more vector-frames (Vektorgeriiste), to each of which belongs a 
quantum number of the total angular momentum and whose mutual orienta- 
tion has, in first approximation, no influence on the energy of the stationary 
states. In the complete energy expression there appear small interaction 
terms which depend, in a not purely additive fashion, on the elements 
(Bestimmungsstiicke) of a pair of vectors and which denote an energetic 
coupling of these vectors.* As a consequence of this coupling every sfate, 
defined by the quantum numbers of the vectors, splits up into a multiplet. 

Consider the coupling of just two vectors. In analogy with the notation 
for Russell-Saunders coupling but without assuming anything about the 
physical nature of the vectors, we may designate their quantum numbers of 
angular momentum by / and s. 

The interaction terms in the complete energy expression may be written 


down in the form 
Q = doa rno, . (2) 


n 


Here X,, and o,, are functions, respectively, only of the elements of the /-frame 
and only of the elements of the s-frame in such fashion that A, and o, are not 
invariant under an arbitrary space-rotation but transform according to some 
linear law. The coefficients a, may contain any combination of rotation- 
invariant functions of the elements of either frame. In practice, 2 is usually 
given immediately in this form or else, when the /-frame and the s-frame are 
associated with different particles, Q contains terms involving 1/72. The 


expansion 
1 1 Va wi (fi °f2) ? 
~o— £{—) of). (3) 
rie rp m rh r\le 


where r, and r, are, respectively, the smaller and the larger of the distances 
r, and re (the distances of the two particles from the center of gravity), 
brings these terms into the required form. 

In an arbitrarily rotated coordinate system K’ a selected X,, call it A, 
assumes the form 


aq 
NO me Diver, 
p=l 


2 There may also be energy terms involving more than two vectors. To avoid complication, 
we shall not discuss these for the present. They are treated by an extension of the method. 
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where \’“ is the same quantity with respect to K’ as A“ with respect to the 
original coordinate system K. Thus there are g—1 other /-frame functions 
A™ (p=2,--+-, g) associated with A‘ and these quantities transform in a 
way which gives a representation of the space-rotation group of degree g. 
If this representation is irreducible we shall say, for brevity: A transforms 
irreducibly. If it is not, one can always reduce the transformation with the 
result that each \ is expressed by the finite sum >; Ay where the Ay“ 
transform irreducibly. Thus it is always possible to express Q in the form (2) 
in such a way that every A, and also every o, transforms irreducibly. 

Since we are dealing with a free atomic system, the series (2) is rotation- 
invariant and therefore can be divided up into parts such that the \ and ¢ 
in each part transform irreducibly with the same degree g.* Thus we may write 


Ge 2 T. (4) 
a=1 
- == Vary, Og), (5) 


where the summation indices have been omitted in the last sum. 

In the quantum theory, the elements which determine an atomic system 
are uniquely defined according to their significance; the elements of the elec- 
tron spin (its angular momentum components) as well as the ordinary co- 
ordinates and momenta have this property. That we consider them as oper- 
ators makes no difference. But, for even values of g, the transforming quan- 
tities are only defined up to a factor +1. Actually, in the series (4), only odd 
values of g occur and, if we set g=2r+1 and 7, = Q,, (4) assumes the form 


x 
Q= )2,. (6) 
r=0 


We are now ready to consider the symbolic representation of the energy 
operators. We have already remarked that 27+1 quantities which transform 
irreducibly under space-rotation may be represented by &*/n"-’. If we repre- 
sent the Ay) by x:7*/1 y:7“1 and the oy) by x27*/2 ye" “2, then Q, must contain 
these quantities in the form 


w, = (VeVi — Yor)". (7) 


For Q, is invariant under rotation and the only simultaneous invariant that 
can be formed of two spinors is (x2y; — 2x1)", where m is a positive integer. In 
general, 2, may be a sum of expressions like (7).‘ 


3. 
Let us now consider the symbolic representation of the wave functions. 
Each state of a multiplet is uniquely defined by the quantum number /j of the 
3 Tt is well known that Yad;,,)\,,') is never invariant for gg’. 
4 In case of an energy term involving three vector-frames, we shall have symbolic terms 


of the type 
(x2¥1 — y2%1)* (x39 — 3X1) B(x 3¥2—¥sx2)7, 


where a, 8, and y must be all even or all odd integers. 
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resultant angular momentum, where 
I—s'Sjslcts. 


Wave functions satisfying the unperturbed (i.e., interaction terms neg- 
lected) wave equation may be composed of products of two functions, one of 
coordinates associated only with the /-frame and one of coordinates associated 
only with the s-frame. There are 2/+1 functions of the first type which 
transform irreducibly under space-rotation and may be represented symbolic- 
ally by &!*" m'-" and 2s+1 functions of the second type, which may be 
represented by &°* "2 no*~ "2. For the state with quantum number j of the mul- 
tiplet there must be 27+1 functions which transform irreducibly under rota- 
tion and are composed linearly of terms &!*" !~" &3**"2 no*~ "2. We intro- 
duce a constant spinor (a, 0) and form the invariants 


P = (Eom. — neti); Q = (am — b&1); R = (ane — dks). (8) 


Then the state with quantum number j is represented by the invariant 


Visi = P*OSR?, (0) 
where 
atp=22Wjaty= 25; 8+ 7 = 2. (10) 
The 2j7+1 functions, which would be required for the study of the Zeeman 
effect, are the coefficients in y,,; of a2-"" b’*”™" (m'’=—j, —j+1,---,7)2 
4. 


\We now proceed to show that every term of (6) makes a contribution to 
the energetic splitting of the multiplet which is, up to a numerical constant 
C, which does not depend on j, a completely determined function of /, s, j, 
and rf. 

According to the principles of the quantum mechanics, the energy split- 
ting E,, due to the interaction term (7), is given by 


E, [rue - [vutedus (11) 


where / means integration and summation over all space and spin coordi- 
nates. Both integrals are invariant under rotation and every term must 
contain (a, 6) and (a*, b*) each to the power 2j. The only invariant of (a, }) 
and (a*, b*) is (a*a+b6*b) to some integral power. Therefore, these quantities 
appear in both integrals only in the form of a factor (a*a+*b)*/ and FE, does 
not depend on (a, bd). 


5 In case there were more than two vector frames composing our atomic system, we should 
form our wave functions of products of all the invariants (:n,—it,) and (an,—béx) in such 
a way that in every term the sum of the powers of £; and 7; is 2/; and the sum of the powers of 
a and b is 2}. In general, it will be possible to do this in a number of linearly independent ways, 
say f, which means that there will be f states of the multiplet with the same j value. To find the 
stabilized linear combinations of these f functions is a perturbation problem involving the 
solution of an f order secular determinant. 
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In (11), w, isan operator which works on y,,;. w, is built up by multiplica- 
tion and addition of operators which involve only the coordinates of the 
l-frame or only those of the s-frame. Accordingly, when we have written out 
all the terms in (11), every term may be put in the form 


Clils 


where J, and J, involve, respectively, the coordinates of the /-frame and of the 
s-frame and are integrals of the type: 


I; 


I 


& *kl+ m.,’ a na’ a £4 t l+n " 
fer “1'e,* my" y 7 my! mie mn! st (12a) 


wf 


The j value appears only in the factor C which contains, in a known manner, 
the expansion coefficients of the various factors in (11). We will show that, on 
account of the irreducibility of the transformations which the factors in /, 
and in Js suffer under space-rotation, both integrals are completely deter- 
mined in their dependence on the various m’s. There remains over only a 
common constant C,, in principle undeterminable by the symbolic method, 
which depends in general on /, s, 7, and the particular atomic system with 
which we have to do and which is of the nature of a radial integral. 

All integrals of the type (12a) are obtained by comparing coefficients of 
like powers of the components of the constant spinors on the two sides of 


fo** + 2’ no** ia! Yoh FH? yr? HiEy® + mono" me ( 12b) 


f (orn — b*E,*)" Cty, — Buy)**(an, — b&,)*! 
e 


= const (a*a + b*b)?"-"(a*¥A + b*B)(aB — dA)". (13) 


The evaluation of this integral follows directly from the invariance under 
rotation of all the factors in the integrand and of the domain of integration. 
The result must be built up of the three invariants of (a*, b*), (a, 6) and 
(4, B) and these quantities appear to the proper powers only in the form 
given in (13). 


5. 

Having seen that the results of the integrations in (11) are determined by 
the transformation properties of the factors and therefore do not depend on 
the special form of the operators which occur in w,, we may simplify the in- 
tegration in (11) by choosing special expressions for the operators x,"** 
y,”-* which transform in the same way. We have only to take care that in so 
doing we do not make the integral vanish identically. 

Replacing x,7** y,7* by operators acting on functions of the symbols &, 
and n,, we shall obtain, in place of (12a) and (12b), integrals of the type 


ee ae 
fer — 


Keren? e (14) 


Str 
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All integrals of this type are found by comparing coefficients on the two sides 


of 


[comm — b*E*)*0'(anm — bE)" = const. (a*a + b*b)7%6,.,, (15) 


where the results of integration are again determined immediately by the in- 
variance. We find for (14) the result 


2p I 
const ( Diteeres (16) 
pt+p 


Thus we find that our choice of operators to replace x,7**, y,7~* is limited to 
such as make y,,;* w, W)., homogeneous in (&;*, ,.*) to the same degree as in 
(E,, n-). The most satisfactory operators fulfilling this requirement are ob- 
tained from the formula 


0 0\' 
(Av, — Bu,)* = (Am — Be(. —+ B— ) , (17) 
OF; On) 


Both sides are rotation-invariant and the x,7** y,7 * are uniquely defined by 
comparison of terms with the same powers of and B. We may remark that 
the order of the factors on the right hand side of (17) is arbitrary since we 
have the same thing still after exchanging them. (17) is the special case p =0 
of the more general formula 


ra) fa] fp 
(dv, — Buy)? = (An, — By,)' te (.1 +B ) , (18) 
0 


gn On) 


THE CONFIGURATION: TWO ELECTRONS, ONE IN AN s-STATE 

6. 

We shall now see how this symbolic method works as applied to spectra 
involving two valence electrons, of which one is in an s-state.® 

The two electrons have angular momentum quantum numbers respec- 
tively / for the orbit, 1/2 for the spin and 0 for the orbit, 1/2 for the spin. The 
2/+-1 wave functions of the orbital coordinates of the first particle, which 
transform irreducibly under rotation, may be represented symbolically by 
the 2/+1 monomials &'*”" n'~™. Likewise we have two monomials S, and S_ 
(the two components of a spinor; hence the name) to represent the spin func- 
tions of the first particle and two monomials S, and S_ to represent the spin 
functions of the second. From the spinors (£, 7), (S,, S_), (S,, S_), and (a, b) 
we form the six invariants: 


(an — dé), (aS_ — bS,), (aS_ — bS,), (Sn — S_E), (Sin — S_8), (S,S_ — S_S,). 


6 Houston (Phys. Rev. 33, 297 (1929)) has treated this case, using the explicit wave func- 
tions and energy operators. The present treatment takes account of the interaction of the spin 
of each electron with the orbit of the other, which Houston did not consider. As we shall see 
in Part II, this interaction is not negligible. 
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These must be multiplied together to form polynomials homogeneous in 
(£, m) to the power 2/, in (S,, S_) and in (S,, S_) to the power one, and in 
(a, 6) to the power 2). 

We find (as we should expect from the vector model) that the only pos- 
sible values of j7 are /+1, 7, /—1. The symbolic functions for j7=/+1 and 
j=l—1 are uniquely determined but there are three functions for j=/, of 
which two are independent. Thus our multiplet will consist of four levels (a 
triplet and a singlet in case of pure Russell-Saunders coupling). We form the 
symbolic functions: 


f=lti1 I =Aj(an — d§)*"(S,b — S_a)(S,b — S_a)} 
IT, = A} (an — b&)?-\(S,n — S_t)(S,b — S_a)} 

j=l IT; = A} (an — bE)?-"(S,m — S_t)(S.b — S_a)} (19) 
IT, = A}(an — b)?(S,S_ — S_S,)} 

j=l—1 711 = A} (an — be)?-2(S,n — S_2)(S,n — S_é)}. 


The symbol A in front of a bracket means to subtract from the expression 
in the bracket the same expression with (.S,, S_) and (S,, S_) interchanged 
and (£, 7) replaced by (£, 7). Expressions in (&, 7) refer to functions of the 
coordinates of the second particle. Thus our wave functions are made to 
satisfy the Pauli principle requirement of antisymmetry. 


A 

Of all the terms in the interaction energy of two electrons,’ we shall neg- 
lect, to begin with, all except Q), the electrostatic interaction,® e?/ri2; Qe, 
the leading terms, 


Ze j ([r1, v1] -s:) ({r2,v2]-s2)) 


| 
—- 


2moc? \ r,° re® ,* 


in the interaction of the spin of each electron with its own orbital motion; 
Q;, the leading terms, | | 
ce? (({ri,vil-Se)  ({re,ve}-s1) 
@ f(lrvid-ss) | (I \ 


moc? \ ry? r,° f 


in the interaction of each spin with the orbit of the other electron. In Q and 
Q; we have neglected all terms but the first in the expansion (3) of 1/712°. 
To find the symbolic representation of Q:, we must use the expansion (3). 
Q, then assumes the form >, @,(r1-r2)", where the a, are functions of the in- 
Variants, 7; and re. Then the symbolic representation is 
doa, (xi — vz)", (20) 
n 
7 See, for example, Heisenberg, Zeits. f. Physik 39, 499 (1926). 
§ Of course, the term e?/7;2 does not properly constitute a small perturbation, since it is of 
considerable importance in fixing what one should like to consider as the unperturbed orbit. 


However, the true perturbation may be considered to be the difference between some well- 
chosen central field Hamiltonian and the true Hamiltonian in which e?/r;. appears. The sym- 
bolic representation of this perturbation is the same as that of e*/rjo. 
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In the matrix integrals which we wish to calculate, the fac tors coming from 
*) and (&, ») [or (& 7) and ( 


4) | to the same power 2/ or (£*, *) and (g, 4) lor (&*, en ) and (é, n)| to the 


same power 2/. Consequently, only the terms 7 =0 and »=/ in (20) will vield 


the wave functions will give us either (£*, 7 


integrals different from zero. The term from n=!/ is the so-called resonance 
term. Thus, for our purpose, the operator may be represented by® 


C; j a) a\?! = 0\?! \ 
OQ = Ao + —4(e— 4+ n—) +(F—435 4 (21) 
(2/)! \ og On 0g On f 


The operator Q2 is in the form to be represented symbolically by 
Tay2t 


C2 (Ny — Yu)? + (Vs — T3)3} 


where (x, v) is associated with the orbital coordinates and (.Y, ¥) with the 
spin coordinates. We find the representation of Q. in terms of the spinors 
which occur in the wave functions by a double use of (17): 


0 0 
(Xv — Fx)? = (Ny — re( x +) ) 
aks On 


, 0 0 _ @ _ @ ) 
~88— +985: +3. 7 
0g On OSs Os 


where the symbol S in front of the bracket means to add to the expression in 
the bracket the same expression with barred and unbarred symbols inter- 
changed. 

The operator Q; differs from Q only in the interchange of (.S,, S_) and 


(S;, S_). 
= o* 
(Sin S_E&){ - 
Os De 0g 


dis 
iz 
ae - a 
ae Os OS 


9 The term for n=/ comes from (18) with r=/ and p= —tr. 


Q; = C38; 
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Applying the operator 
2 =%+2% +29; (24) 
to the wave functions (19), we find 


Ql = {495 +C, —UC2+C;) {I 


QM, = {Ao — e+ (UL + 1)C3} Ta + (Ci + Co} Is 

(25) 
Qs = }Ao + (+ 1)Co — ICs} IIs + [Ci + Cs} Ie 
QUT = {Ao + Ci + + 1)C2+Cs) {IM . 


The multiplet energy levels €; and ¢€;, defined by 


«fi = [ver « f uri = f mom, (26) 


are given immediately by (25) without integration: 
€, = lp + Ci — UC2 + Cs) 
€3 = Ag “hk C; “+ (1 + 1)(Ce _ Cs). 


To find the other two energy levels and to find the linear combinations 
of I], and II; to represent the stabilized zero-approximation wave functions, 
we must solve the second order secular perturbation problem. The stabilized 
functions may be written in the form 


dg: = Yolla + dolls 
(28) 
oy = villa a 64113 ° 
Then our Y2, 62 and y4, 5; are the two solutions of 
Y ss — €Aca) + 6(Qas = €Aa3) = () 
(29) 


¥(Qsa — eAga) + 6(Q33 — eAss) = 0 


corresponding to the two roots €, €, of the determinant of the coefficients in 
(29). Here we have defined 


Qas — fueron; Aas = fet en. (30) 


Making use of (15) and of the integration formula 


few — b*E*) o+t(aE* + bn*) "(an — dE) OT! (aE + D*y) -" 


) -" 
= const. il ) (a*a + b*b)*4, (31) 
il 


which is obtained by replacing a by a+/fb* and b by b—fa* in (15) and equat- 
ing the coefficients of like powers of ff* in the result, we proceed to carry out 
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the integrations (30). It is convenient to integrate first over the spin variables 
and then to make use of the relation 


(a*a + b*b)(E*E + n*n) = (a*n* — b*E*) (an — bE) + (aE* + bn*)(a*E + b*n), (32) 


which leaves us integrals of the type (31). We find 


Aaa = Asz3 = const. (2] + 1)(a*a + b*d)?'- (1/7) 
(3.3) 
Aas = Asa = const. (a*a + 6*b)?'- (1/7). 
The integrals Q., etc. are given in terms of these by (25). 
For the secular determinant, we find 
1 Cy ACe+ (274-1) flo — C2 + (14+ 1)C3—-€ | 
(272 +1)C, + (/+1)(Co+C3)+A0—-€ 1 
=(). (34) 
(27+1)C,;+(4+1)(CotC3)+Ao—e 
Cy +C 34+ (2/41) } Ao +(4+1)C2—-IC3-€} { 
The solution is 
e = lo + (3)(C2 + C3) + (UL + 1)(C2 — Ca)? + [C1 + (A) (Co + C3) J?) "2. 
If we define the quantities p and @ by 
C; oo (3)(Coe oo Ca) = pcos 20 
(35) 
[72 + 1)] (C2 — Cs) = p sin 28, 
we have, as the roots of (33), 
€o = Ao + (3)(C2 + C3) + p 
; (36) 
‘= Ao + (3)(Ce + C 3) = Pa 
The complete multiplet, given by (27) and (36), is 
€: = Ag + C, -—- (Co + C3) 
€2 = Ag + Ci + (C2 + C3) + 2p sin? 
(37) 


€; = Ag+ C; + (/ + 1)(Ce + C3) 
és = Ag a C; —_ 2p sin? @. 
We observe that, in the limit @=0, €,, €, €; form a Russell-Saunders 
triplet with the intervals in the ratio /+1 to/, while €, appears as a singlet at 
a distance proportional to the resonance energy, C;. Using these limit values 


of € and & in (29) and substituting in (28), we find for the pure triplet and 
pure singlet linear combinations of II, and II; 


N (Ua + Is) 
N,(Ila — Is) = N,IIy. 


gr 


(38) 
ds 
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From (33), we find for the normalizing constants 
N, = const. /'/*;N, = const. (J + 1)!/?. (39) 


¢:, like I and III, is antisymmetric in the coordinates, symmetric in the spins; 
¢, is symmetric in the coordinates, antisymmetric in the spins. 

Now, substituting for II, and II3 in terms of ¢, and ¢, in (28) and using 
the values of y and 6 obtained by putting (36) in (29), we find for our stabil- 
ized wave functions 


i sin ] Ps + cos 7] oi 
cos 8, + sin 0 @;. 


de 
os 


(40) 


Here the normalizing factors have already been introduced. 





€ jrlel 

















es j 
Fig. 1. 
Introducing the notation dj. = €;— €, etc. and using the symbols A and D 


to represent the distances from the normal position of the center triplet level 
to €2 and & respectively (see Fig. 1) we find from (37) 


L 
II 


| [+ 1)de3 — Idig] /(21 + 1) 
A= |dy|—A 
[( + 1)de3 — ldye|/| (20 + 1)de4] 
C, = (3)(dos + dig — dos) 
Co+Cs3 = — (dig + do3)/(21 + 1) 
"9 "3 = (AD/I(I + 1))'”?. 


tS 
—s 
~ 
— 
te 
SS 
II 


(41) 


oO 
| 
Oo 
II 


Since sin® @ must be positive, the sign of (/+1)d2;—/diz must be the same 
as the sign of da;. This means that, if the singlet, €;, lies below the triplet, the 
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center triplet level, €, must lie above its normal position and vice versa. A 
convenient way of stating this rule is that: the center triplet level is repelled 
by the singlet level.*° 

In (41), the quantities C. and C; are not completely determined since 
we do not know the sign of C2—(C;. They become determined when we make 
the physical assumption that the interaction of the spin of each electron with 
its own orbital motion is larger than with that of the other electron, i.e., 
'C,'!> |C;|. This is equivalent to making the sign of C:—(C; the same as the 
sign of C.+ C3. This determines also the sign in (40) of the angle, @, which was 
defined by (35). The sign of C.+C; (and so likewise the sign of C2—(C3) is 
always the same. Since we can express C2— C; in the form 


Co — C3 = doy sin 6 cos @ uC + 1)}} ;. 


we must change the sign of @ according as dz; is positive or negative, i.e., 
according as the singlet lies above or below the triplet. 


8. 

While € is regarded as a triplet level in spectral classification and ¢€; 
as a singlet, we observe in (40) that the wave functions for these states are 
partly singlet and partly triplet in character so that we may expect transi- 
tions between both of these levels and other singlet and triplet levels—i.e., 
we shall expect forbidden intercombinations between singlet and triplet 
systems as well as combinations within these systems—with relative inten- 
sities determined by @. 

The intensities of transitions between levels are proportional to the 
squared absolute values of the corresponding matrix elements of the polariza- 
tion. Kramers has shown! how the symbolic method is adapted to calculating 
the matrix elements. The operators for the components of the polarization 
must have the transformation properties of the components of a space-vector. 
These operators may all be represented symbolically by the spinor invariant 


Q = (AY — BX)’, (42) 
where we make the association, according to the transformation properties, 
P,+ iP,— X*; — P,+ iP, Y?; P:. > XY. 


Brinkman! has shown in detail that calculation with the wave functions 
corresponding to Russell-Saunders multiplets leads to the well-known in- 
tensity rules as given by Kronig and by Hénl.'* Hence, for our purpose, we 
may regard the matrix elements associated with the functions I, ¢,, III, ¢, as 


10 We must not expect this rule to hold in case of very small displacements of the center 
triplet from its normal position, since this may arise from cancellation of the opposite effects 
of the perturbations { and Q, in which case the effect of the other perturbation terms we have 
left out may not be negligible in comparison. 

1 Brinkman, Dissertation Utrecht 1932. 

® Kronig, Zeits. f. Physik 31, 885 (1925); 33, 261 (1925). Hénl, Ann. d. Physik 79, 273 
(1926). 
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known. We may indicate the matrix elements for the group of transitions 
I+1-—/ by 

J 1980 = Q),, etc. (43) 
where quantities referring to the final state are distinguished by a dash. We 


shall make explicit use of the fact that all matrix elements Q,,, Q,,, etc., and 
. ! - | . « 
all matrix elements for Aj | >1 vanish and that* 


I+1 
nme ERg (44) 
[72 + 2) ]1/2 


Then we obtain, for the actual multiplet, the matrix elements 


fae =O 
fear = cos6éQ,, 
f mer’ = Os) 


| 6,21’ = sind Q,, 


=-ss 


| $222’ = J oneteiin — sin @ sin 6’ + cos 6 cos Pat Qe, 
SUPA + 2) 2 f 


| IIIQ¢2’ = cos &’ Qs, (45) 
_—— +1 hci a \ 
{ O40. = < — — — cos @sin 6’ + sin@ cos 0 + Qy, 
[m+ 2) pee 
, j l + l ° ; ; ° a 
go2b, = 4 — — — sin@cos# + cos@sin & -Q:, 
(fa + 2)] 2 f 


[ 11126.' = sin 0’ Qs, 


f ?+1 — 
$;2¢,' = <—— —— ¢os 6 cos 6’ + sin 6 sin 0’ > Qs, 
Ufa + 2)? 


a 


| LIQUit’ = Qs. 
The sum rules for intensity hold for the complete multiplet and not for the 
singlet and the triplet separately.“ In applying the formulas (45), we must 
use the rule at the end of section 7 for determining the signs of @ and @’. 

'8 This relation is given by carrying out the symbolic integrations. Apart from the sign, 


it follows from the sum rules. 
4 Compare with Ornstein and Burger, Zeits. f. Physik 40, 403 (1926), 
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The formulas for the relative intensities of spectral lines, which are ob- 
tained by squaring the expressions in (45) and substituting in the Kronig- 
Hénl values for Q) 7, Qn *, ete., hold only in a sort of first approximation, 
i.e., With neglect of quantities of the order of magnitude Avy, where Av is the 
interval between the frequencies of the lines compared. In many cases this 
ratio is not small and the ambiguity may be correspondingly large. Thus we 
should not attempt to compare our results with experimental intensity ratios 
except for groups of lines of about the same frequency. The origin of this am- 
biguity lies in the fact that the matrix elements calculated, while assumed to 
be those of the polarization, might equally well be those of, say, the second 
time derivative of the polarization; for this quantity has the same rotation 
properties. But the squared matrix elements for these two quantities differ by 
a factor v*. 

As remarked at the beginning of section 7, we have taken account only of 
certain selected terms in the interaction energy of two electrons as given by 
Heisenberg (reference 7). With the exception of the leading term in the spin- 
spin interaction, all of the neglected terms involve more than two vector- 
frames. To show how the method is applied in such cases, we treat now the 
leading term (i.e., we take only the first term in the expansion (3)) in the in- 
teraction 

e* j([r1, v2] s))+ ({r2,vi]-s2)) 
, . en (40) 
moc? ] ro" f 
The expression ([{1r, v2|-s:) involves the orbital motions of both particles and 
the spin of one of them. To obtain a symbolic representation, we associate 
each of the vectors with a spinor: 
ri — (x, ¥); ve > (XY, 1); si (r,s). 


We may make the association in this way. If (a, 6, c) are the components of 
a vector and (£, 7) are the components of a spinor, we may set 


l 
cer =f Fe ~-—04+ 93; Co — &. 


This gives us, as the symbolic representation of ( [ro, v»| So), 
(rY¥ — sX)(ry — sx)(Xy — Fx). (47) 


We proceed to replace each of these spinors by operators involving the 
spinor representing the corresponding vector-frame in the wave functions. 
Thus we substitute 

for (x, y) expressions involving &, n, 0/0 
for (Y, Y) expressions involving &, 7, 0/0&, 0 
for (r, s) expressions involving S,, S_, 0S, 0/0S_. 


On; 


We observe that, in general, it will not be possible to find an operator repre- 
senting this interaction term whose matrix elements do not vanish identically. 
Only in case /=1 do we get integrals different from zero. The operator may be 


1 £, » are the coordinates by stereographic projection of a point a, b, c on the zero sphere. 
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put into a form where it contains first derivatives with respect to S, and S 
and either second derivatives with respect to £ and y and no derivatives with 
respect to £ and 7 or second derivatives with respect to & and 9 and no deriva- 
tives with respect to € and 7. We then obtain nonvanishing resonance in- 
tegrals when / = 1. 

To get (47) into a form where we can use (18), we introduce a symbol 7; 
which means to pick out the coefficient of f in the expression following. 


(ry — sx)(Xy — Vx) = 37;)(N + fry — (V + fsx}? 
a __¢ _@ \ 3 0 aw) _@ _°¢ 
— 3754 (X + fr)— + (V + fs)—p -(, +s ) a= T- ), 
dé On f d& On 0g On 


Applying (18) a second time™ 


; ; 0 0 j 0 Ai te 
(7 — sNX\r + 5 = 37;,(V +f )y —-{[X—/ Ss? 
0g On ] ~ OF On f 
iT (1 + ; *)s (x “\s \ (x =) 0 +(1+ f =) a) | 
— al; y* -~~ta =~] . ct 
— \ ~ OE ~ On | ~ On/ dS, © OE/AS 
j a o* 0 fa] 0 ra] } 
= sfour sf) (52 452\xt ort) 
] d&aS Onas , dg On, Os, as f 


In the final step, we have simply Y¥—£, Y— 7. This gives us the symbolic 
operator for (47) in the form 


{ . 0? a 
4 (Sif — S_&) — -_- — — 
T 0&aS OndS , 
ra) ra] ra] ra] | ra] ; 0 
+(s,—+s—\i— +4 —)}(i—+4-). 
0g On Os, As dé On 


The other usable form of the operator differs from this only in the interchange 
of (&, 9) and (§, n). The term ([r2, v: | - se) vields operators which differ from 
these two by a complete interchange of barred and unbarred quantities. Ob- 
taining symmetrical operators by lumping together those operators that 
differ from each other by a complete interchange and putting them in a 


f 


Str 


somewhat different form, we find as the operator representing the leading 


term in (46), 
( 0 oe ) 
OS, On OS_OE 


 ffe8 , 8 
w=C,-S4[E—+ 79 (S49 — 5 
( eke On 


Ser' 


, 0 0 _@ _ @ | 
~-25— ¥ 8-3 —— + 3..— ( 
0g On OS , Os 
«all 0 0 ; 0° 0° 
+ C,:S: (: + n- Sa + 3. - 
rks On OS ,O7 OS_OE 


0 fa) ra) ; 0 ) 
-{ _ tT S4— 3 -* 
0k On Os, Os 


16 9/0& transforms like » and 0/dy like —£. 


(48) 


tr 
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In the special case ] = 1, we apply the operator (48) to the wave functions (19) 
and find 
el = (Cy a C;)1 


— Cyll, + (Cy — 2C;5) 13 
wills = — C;II3 + (C's — 2C,)11, 
wlll = (Cy +C;)111. 


wll, 


In addition to neglecting what we believe to be the unimportant terms in 
the interaction energy, we have made another omission which may in some 
cases completely invalidate our results. \We have assumed in our treatment 
of the interaction of two electrons as a perturbation that each electron con- 
figuration (defined by the # and 7 values for the two electrons) may be 
treated quite by itself. This treatment is valid only so long as there are not 
similar terms of nearly equal energy arising from different contigurations.™ 
As we follow through a series of terms from configurations differing in the » 
value of one of the electrons, irregularities in the values of 7* reveal the pres- 
ence of perturbing terms arising from other types of configuration. Shen- 
stone and Russell have shown that perturbations of this sort are common. 
We cannot expect our formulas (41) for the energies of the different inter- 
action terms as functions of the multiplet intervals to have much meaning 
for a configuration which is much perturbed by a neighboring one. 

Part II of this paper, containing the application of the theory developed 
in Part I to observed spectra, will appear later. The relative energy values 
of the different types of interaction will be calculated for many spectra. The 
calculated intensities of certain multiplet lines in the spectrum of mercury 
will be compared with the experimental values. 

The writer wishes to express his thanks to the Lorentz Foundation for 
the grant of a fellowship and to Professor Kramers for suggesting this prob- 
lem and for many discussions while the work was in progress. 


'7 This was first emphasized by Condon, Phys. Rev. 36, 1121 (1930). 
18 Shenstone and Russell, Phys. Rev. 39, 415 (1932). 
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The Isotope Shift in Hyperfine Structure 


By Jenny E. ROSENTHAL AND G. BREIT 
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The effect of deviations of the nuclear field from the inverse square law such as 
are suggested by the theory of radioactive disintegration is calculated. The isotope 
shift due to a change of the nuclear radius by 1 percent is found to be larger than those 
observed. The spectrum of Tl appears not to fit any simple modification of the theory. 
The isotope shifts in Pb I and Pb II can be explained roughly on the supposition that 
the deviations of the nuclear field from the inverse square law are local but that for 


the valence electrons these deviations should be taken to be smaller than is usually 
supposed in the theory of radioactive disintegration. 


$1. INTRODUCTION 


N THE analysis of the hfs of T1,)? Pb**5 and Hg*’? there have been ob- 

served isotope shifts of the order of magnitude of 0.05 cm in the arc 
spectra to 0.5 cm~! in the spark spectra. It has been pointed out! that this dis- 
placement is greater than that to be expected on account of the mass correc- 
tion to the Rydberg constant and that it might possibly be due to deviations 
of the nuclear electric field from the Coulomb law. The effect of such devia- 
tions has been investigated by Bartlett® with the apparent result that they 
may easily account for the observed isotope displacement. In view of the fact 
that Bartlett’s estimate does not take into account the relativisitic nature of 
the electron and the screening of the nucleus by the other electrons it has 
been thought of interest to make a more detailed calculation. Since this work 
was completed there appeared a letter of Racah® in which the same effect 
is estimated for S terms of Hg with the result that it is hard to explain why 
the isotope shifts are as small as they are. This is in agreement with our cal- 
culations. In the course of the work the effect has been estimated for other 
terms as well with the same result for some and quite different results for 
others. Also the effect of the deviations of the nuclear field from the Coulomb 
law on the interaction of the electron with the nuclear spin has been estimated. 


1H. Schuler and J. E. Keyston, Zeits. f. Physik 70, 1 (1931). 

2 D. A. Jackson, Zeits. f. Physik 75, 233 (1932). 

3H. Kopfermann, Zeits. f. Physik 75, 363 (1932); Naturwiss. 19, 400 (1931); 19, 675 (1931). 

4 J. L. Rose and L. P. Granath, Phys. Rev. 40, 760 (1932). 

*> H. Schuler and E. G. Jones, Zeits. f. Physik 75, 563 (1932). 

6 H. Schuler and J. E. Keyston, Zeits. f. Physik 72, 423 (1931). 

7H. Schuler and E. G. Jones, Zeits. f. Physik 74, 631 (1932). 

8 J]. H. Bartlett, Nature 128, 408 (1931). 

®° G. Racah, Nature 129, 723 (1932). In this letter the potential inside the nucleus is sup- 
posed to be constant and at the boundary of the nucelus it is taken to be continuous. Below the 
potential is taken to be discontinuous at the nuclear boundary so as to represent the potential 
wall of nuclear models for a-particle emission. In comparing the formulas in this paper with 
those of Racah one should start with (11) since in (11a) Av is supposed to be zero. The resultant 
isotope shift becomes then 1/(29+1) of (11a) with v=0 in the latter formula. 
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$2. SOLUTIONS OF THE WAVE EQUATION 


The central field equations for the Darwin-Gordon radial functions ¢, 
and @» for Dirac’s equation are: 


(dg, dv) — jo. v= (da)[1 — (E— V)/me2|de 
(1) 
(dg. dy) + j’o: v = (}a)[1 + (E — V)/me? Jo, 


vy = 8r2*mZe*r hh? = 2WZr/an; a = Za = laZe*, he 


Z is the atomic number, ay, the radius of the first Bohr orbit for hydrogen, 
FE is the energy and V the potential energy and j’ is the quantum number 
(analogous to Sommerfeld’s &) taking the values —1, +1, —2, +2, —3, - 
for S)2, Pie, Ps2, Ds2, Ds 2 > > + terms respectively. 

We consider first the solutions of Eqs. (1) for y<yo, where vo is the value 
of y corresponding to the nuclear radius 7». If we assume a Gamow field inside 
the nucleus, V is constant for y< yo. A satisfactory approximation is obtained 
by letting E/ mc? =1 which means that the term value energy is negligible in 
comparison with mc*. The solutions can then be written: 

dy = sy! *J) pia (v'*(e — 2)'?y/ (2a)) 

de = F sv'2(y — 2) Ay, py ye (o'/2(o — 2)'/2y/(2a)) (’ z 0). 
The positive index of the Bessel functions is required by the condition of in- 
tegrability; v= V/mc? and s are constants, y<yo is small and if |v is not 
very large the Bessel functions are sufficiently accurately represented by their 
first term. The values of v for which this is valid may be taken to be given by 
—20<v7<20. In this approximation: 


od, = ky’; be = R22 — v) y't/[2a(2j’ + 1)] Gj’ > 0) (3’) 


by, = k® yt li'lp/[Ja(1 + 2) 7’| )]5 be = R@ yl’! Gj’ < 0) (3”") 


and 
1 2 = 2a(2j’ + 1)/[(2 —a)y] VG’ > 0) 


o:/o2 = vy/|2a(1 + 2) j’\)| Gi’ < 0) 


k and k® are constants which will be determined later. 
In the region y>yo, V is the Coulomb potential and the equations for 
go, and ge are: 
(do;/dy) — j’¢:1/y = (1/2a)[1 — (E/me?) — 2a*/ yoo 


(5) 
(do2/dy) + j'¢2/y = (1/2a)[1 + (E/me?) + 2a*/y]91. 
Letting as before E/ mc? = 1, we have as the solution of Eqs. (5) 
o:/a = CJo,(2y!/*) + C_J_2,(2y'!*) 
(0) 


de = CAs,(2y"!2) + C_A_2,(2y!/) 
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with the abbreviations: 
lo, Ci’ — p)J2p + v" Joost 


Aa, = (j" — p)I-2 — YS oy 


where C and C_ are constants. 


The condition of continuity across the boundary is: 
Pil Vo —-eée = oi(vo + €), Pol Vo —e) = o2( Vo + €); e—(). 


Substitution of the actual solutions for @; and @e into (8) gives: 


Cl(o, dade) Ao, o“ Jo, | = C [J 2p — (¢) dds) A 2,1. (9) 


The argument of the functions J, etc. is 2vo'/*; o1/ade is given by (4’) or 
(4’’) depending on whether j’>0 or j’<0. Condition (8) also gives Rk 
(respectively &’) in terms of Cand C_. 

If the field were Coulomb at the nucleus, C_ would be zero since the con- 
dition of integrability for @: and ¢2 must be satisfied. The presence of the 
Gamow field results therefore in a modification of the eigenfunctions repre- 
sented by the term in C_. These modified or perturbed eigenfunctions will be 
denoted by @ and @2; the unperturbed functions (i.e., functions with C_=0) 
will be written as ¢," and @,»". 

The ratio C_/C as given by (9) was computed for a wide range of values 
of v in order to see to what extent the wave function is affected by deviations 
of the potential from the inverse square law. The results of this calculation 
for 7/=+1 are tabulated below. They are accurate within 5 percent for 
—20<v7<20. The formulas also may be used to compute C_/C for values of 
v lying outside these limits but the results will be less accurate due to the 
fact that for large values of |v! (3’) respectively (3’’) are no longer satisfac- 
tory approximations to (2). For j’=—2, C_/C was estimated to be quite 
negligible in comparison with those for j’= +1. 


TABLE I. Z=81 
ro = 107" cm; vo = 0.03, a = 0.590, p = 0.807 
Jo, = 0.0403, Jo,41 = 0.0027, J_2, = — 4.736, J_x-1 = 43.0 


> 2p 
As, = 0.0082, A_», 8.356, o:/a¢2 = 200/(2 — 2) 
C(1.58 + 0.04032) (4.7367 + 1662) 
Ae, = —0.0724, A_2, .116, o)/dee = 
C(0.07240 + 2.81) "(1.1167 + 330.0 
— 20 = 0 4 20 
1); 0.0005 0.0009 0.0009 0.0010 0.0014 


— 1); 0.0044 0.0077 0.0085 0.0093 0.0120 
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Z= 91 
ro = 10°" cm; vo = 0.035, q = 713/4, p= 
Jo, = 0.0000, Jo4,1 = 0.0044, J_2, = — 3.722, J_2,_, = 28.62 
j =1, lo, = 0.0158, A_s, = — 6.286 
C_(5.440 + 0.12597) = C(2245 + 7.8067) 
fe <-' Aw = — 0.1042 A_o, = 1.158 


C_(7.806 + 0.03242) = C(0.1259 + 0.0029) 


v: — 20 —4 0) 4 20 
C_/C(y’ = 1); 0.0014 0.0022 0.0024 0.0026 0.0033 
C_/C(j’ = — 1); 0.0095 0.0149 0.0161 0.0173 0.0218 


$3. Tue IsoTrorE DIspLACEMENT 


Corresponding to a change 6V of the potential inside the nucleus there is 
a change 6W in the energy. The usual way of evaluating it is by means of the 
perturbation method using the unperturbed functions ¢," and ¢2" correspond- 
ing to C_=0 





OW = | (1 + (Ze" r) |(o. 2+ go? )dr (10) 


where the normalization is such that 


x 


{ (pi? + bo *)dr = 1. 


To evaluate (10) ¢:°° and ¢o are expanded in terms of y. For small values 
of y such as 0<y<yo these expansions converge very rapidly, and for a first 
approximation it is sufficient to use only the first term of each expansion. 


G12 + po? = 2C27’"(7’ — p) y#/T?2(2p + 1). 
Performing the integration we get: 

6 = 2Ze%XC2j"(j’ — p)[(yo%/2p) + (v/2a2) vo2*t?/(2p + 1)]/T2(2p + 1). (11) 
Supposing that yo changes by Ayo and that v remains unchanged 6W changes 
by: 

ASW = 2Ze*C2j’(j’ — p)| v0 + (v/a?) yo2*t!|(Ayo/vo)/P2(2p + 1). (11a) 


For S terms the magnitude of the isotope shift can be conveniently ex- 
pressed in terms of the value of ¥?(0) which the square of the nonrelativistic 
Schroedinger function Y would have at r=0. The normalization of the 
function yw is for S terms such that 


inf y?-r2dr = 1. 
0 
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The function ¢» for all but very small values of r satisfies the same differential 
equation as ry. It may be equated to the latter provided ¢: and ¢ are nor- 
malized so as to have 


a 


inf (o:> + o2*)dr = 1. 
0 


The expression for AéW.given by (11) must then be also multiplied by 47. 

The leading term in the expression for ¢2 at not too small y is Cy"? Jop41 
(2y'?). Nonrelativistically the function ry must have the form —C’y"?J, 
(2y''*). The constants C and C’ should be the same because the asymptotic 
expansions of the Bessel functions differ only by a constant phase in the trigo- 
nometric term. \WWe have thus 


y(0) = —- C(v/r) = — C(2Z/an). 
Substituting into (11) and noting that 
2Ze*C? = (e? 2anZ)ay*p*(0) = (R, Z)an*V7(0) 


where R,, is the energy corresponding to the Rydberg constant we have 


4rR, ; l+p v Avo 
M6 = —— an*y2(0) ———— } yo? + — yor! | — (12) 
Z [r(2p + 1)]? Yo 


which is in agreement with Racah’s® formula. 
Supposing 7=0 and using Racah's'® value ¥?(0)=0.49X10” cm-* for 
the 7s term of Tl I we obtain for it 


ASW & 360(Ayo/ yo) cm™!. 


If Ayo/ vo £0.01 the isotope shift for this term should be 3.6 cm~! which is 
much too large according to Schiiler and Keyston! and to Jackson.* This 
value may conceivably have to be multiplied by a factor of about 1/2 (see 
$4 below) in order to take into account the change in the wave functions due 
to the finite size of the nucleus and it may have to be multiplied by 1/3 to 
take into account the fact that the nuclear radius would be expected to vary 
with the cube root of the atomic number rather than with the first power. 
But even 0.6 cm~! appears much too large for this term particularly when it 
is considered that according to Schiiler and Keyston the isotope shift for it 
is considerably less than that for 7 *Pi,2 which itself is about 0.05 cm~*. It 
should be noted that the magnetic moment of the Tl nucleus as calculated!!! 
using the splitting of the 62S term and the value of ¥?(0) used above is 
smaller than that calculated using the splitting of 7 *P\,2 by a factor of about 
4. It appears that the interaction of the 6s electron of Tl I with the nuclear 
spin as well as with the Gamow trough is smaller in reality than would be 
expected by comparison with other terms of the same spectrum. 


10 Racah, Il Nuovo Cimento 8, 178 (1931); Zeits f. Physik 71, 431 (1931). 
1G, Breit, Phys. Rev. 38, 463 (1931), see p. 471. 
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For P, D,-- - terms of single electron spectra the magnitude of C can be 
obtained from the doublet splitting of these terms as has been already 
done in the discussion of the hyperfine structure by Racah.'" Using Sommer- 
feld’s doublet splitting formula one obtains for the energy difference between 
two parts of the doublet 


5 = (+ 1)(2l + 1)Zpo?} | + 1)? — a? ]2? — 1 — (2 — a®)8 2h a 29) see 
(uo = eh (4amc)) 


where the suffix “Schr” indicates that the average of r~’ is to be taken using 
Schroedinger's nonrelativistic wave functions. If the relativistic effects were 
absent we would have 

o2 = rw =Clim As, = + Cy Vayi(2y'?) (7 Z 0). (13) 

10 

For sufficiently large values of y there is no difference between the Dirac and 
the Schroedinger wave equations. The amplitude of the oscillations of the 
wave function as given by (7) is the same as that given by the nonrelativistic 
Eq. (13) provided the same values are given to C in the two expressions. The 
value of C is thus determined in terms of (7~*) by 


x 


(r-8)schr = C2(32Z2ay yf Jory i2(s)s-8ds = 2C%ZPay-?/ | + 1)(27 + 1)] 


where ay is the Bohr radius. Using now the relation between (7~*)s.4, and 6 
and the expression (11) we have 


6W = 2)’ — p) [F204 1) )?{(@+ 1% - a@}?-1 
— (2 — a®)¥2t-1[(ye2#/2p) + (2/2a2) vo2et/(2p + 1) |. 


Keeping v constant the change in 6W due to a change in the nuclear radius 
corresponding to a change Ayo in yo is 


AdW = 2)’ — p) [P20 + 1) 2h [+ 12 - a@]!2 - 1 
— (12 — a2)'/2}—1[ yg? + (2/2a2) vot! |(Avo/ vod. (14) 


For the P;2 normal level of Tl ! substitution of j7’=1, p=0.807, yo =0.03, 
6=7792 cm~! gives 
ASW = 49(Ayo/ vo) cm~!. 

The difference in the term values of P;,2 for the isotopes 203 and 205 would be 
estimated therefore to be of the order of 0.5 cm~! if it is supposed that the 
nuclear radius differs for these two isotopes by 1 percent. The isotope shift for 
P32 should be negligible in comparison. According to Schiiler and Keyston! 
the difference in A(6W) for P12 and P3,2 is not greater than about 0.05 cm7! 
and according to Jackson? it is probably of the order of 0.01 cm~'. The cal- 
culated value is, therefore, considerably greater than that observed. 

Essentially the same result is obtained using direct numerical calcula- 
tions'! for the wave functions ¢:, de. One finds that at the nucleus ¢);7+¢2" 


=9.82X10-* y?*/ay which leads to 


AbW = 60(Ayo/ yo) cem=. 
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$4. VALIDITY OF THE PERTURBATION METHOD 


Although 61 is small the fields responsible for it are large so that it is 
not clear without an additional calculation to what extent the first order of 
the perturbation method is a valid approximation. It will be seen that a cor- 
rection factor of the order of 1/2 must in fact be applied for the calculation 
of A(é61V). This correction is clearly not large enough to bring about agree- 
ment with experiment for the S term of TI. 

As long as one is interested in changes of the nuclear radius the only thing 
which matters for A(6W) is the value of ¢,;°7+ 2? at y= yo because A(61V’) may 
be considered as the change in the energy due to the perturbing potential 
(Ze*/r)+V applied in the region ry9<r<ry+Ary. It is, furthermore, clear 
from the numerical calculations given above that the wave functions are 
practically unaffected by this perturbation. It is, therefore, sufficient to mul- 
tiply the values so far obtained for A(6W) by the ratio of the (perturbed) 
value of $:°+ 2" at y= yo to the approximate (unperturbed) value which has 
been used so far. This ratio is found by observing that the constant C in (6) 
must be practically independent of yo and v for the range of values considered. 
In fact for large y the asymptotic expansions of the Bessell functions make 
(6) equivalent to 


do = vite 21C2 + C2 — 2CC_ cos (2p + 1)r]!/2 cos (2y"!? — &) 


where € is a constant phase (independent of y) which is immaterial for the 
present purpose. The function @2 is therefore essentially an oscillating func- 
tion the amplitude of which depends on C, C_. The dependence on C_ is seen 
however to be very slight in view of the small values of C_/ C. The amplitude 
of the oscillations for large y must be the same whether the nucleus has a 
finite size or not because it is determined by the continuation of the wave 
function from still larger values in towards the nucleus. We will thus suppose 
that C has the same value whether the nucleus gives a purely Coulomb field 
or whether it has a finite size. The factor by which the previously calculated 
values should be multiplied is then by (6) 


[a2Jo,? + o,2 + 2(C_/C)(aV oJ 2, + Asp-A_2,) 


+ (C_ C)*(a*J +" + A 20°) | (u*J 9," + Ae,*). 


For Tl P; 2, S; 2 terms this factor has the values: 


Pi 2(j’ 1): 1 — 422.8(C_/C) + 1.22 & 10°(C_/C)* 
Si. (j’ = — 1): 1 — 50.76(C_/C) + 1561(C_/C)?. 
Using the values of C_/C tabulated above it is seen that the correction factor 


is generally greater than 1/2 and lies between 1/2 and 1. For any given v 
it is readily computed by the above formulas. 


$5. EFFECT OF NUCLEAR FIELD ON HyPERFINE STRUCTURE 


The coupling of the electron to the nuclear spin axis is known to depend 
on the radial integrals'’"! 
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I -{ didov “dy. 
0 


It has been seen in $4 that ¢@:?+¢,? in the neighborhood of the nucleus is 
smaller than its unperturbed value ¢,“°?+ ¢2""? corresponding to an inverse 
square field of force. It may be expected that @:¢2 is also smaller numerically 
on account of the changes brought about in the functions @1, ¢2 by deviations 
of the field from the Coulomb law. The dependence of the integrand of J on y 
has been examined graphically for Z = 91, 81, y= —4. The main modification 
due to the finiteness of the size of the nucleus consists in reducing @:¢2\~? 
to a relatively negligible amount for ySvo. For vy>yo the integrand ap- 
proaches rapidly the value which it would have if the nucleus had negligible 
dimensions. An inspection of the figures indicates that the area under the 
perturbed ¢@i¢2y~* curve is approximately the same as the area under the 
unperturbed ¢:'"¢2'" y-? curve from y=2y9 to y=. We have for the un- 


perturbed functions 


x x 


[ $1 G2 yA"dy = [ SaC?| (j’ — p)Jo,(s) + (2/2)F294.1(5) |Jo,(s)s ds 
7 7 0 

= aC?(2j’ — 1)/[2p(4p? — 1)] 
f $1 ga y2dy SB aC2(j’ — p)(pyo)?-!/ [(2p — 1)(P (29 + 1))?]. 


Thus the integral over the perturbed functions is 


, [f° — i. 2(j’ — p)p(2p + 1) (bye)? il “ 
= idey “dy = cide er eeneetx opumemeera Th! allie a 
et) OG ~ ees OP ‘ s 


) 
oo 
{ $1 go yAd y 
0 
f 


where p=2. The quantity in , j is a correction factor which must be ap- 
plied to J as evaluated on the hypothesis of the inverse square law. It is seen 
that the correction factor depends in addition to yo also on j’ and p. It is 
largest for S terms and is appreciable for P:,2 but is negligible for P32 terms. 
Taking p=2 for Z=91 the correction factor is 0.67 for Si,2 and 0.86 for Ps) 
terms. Schiiler and Keyston® and Schiiler and Jones’ find apparently that 
the hyperfine structure splittings of the two odd isotopes of Hg are not al- 
ways proportional to each other as one would expect on a simple theory, the 
splitting of 199 being sometimes greater and sometimes less than that of 201. 
One might suppose that this is due to the difference between ¢:¢2 and 
$1°ds. The above estimates show however that the effect is not accounted 
for by the theory because the nuclear radii of the two isotopes must differ 
from each other by much less than 1 percent as shown by the smallness of the 
isotopes displacements. The resultant difference in I is then expected to be 
surely not greater than 0.33 percent. Nevertheless for the absolute calcula- 
tion of magnetic moments the above correction factors may have to be con- 
sidered. 
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$6. SUMMARY AND DIscUSSION 


1. The observed isotope displacements in the spectra of Hg, Tl, Pb are 
smaller than those calculated on the supposition that the nuclear radius 
changes by about 1 percent for isotopes differing from each other by 2 in 
their atomic number. 

2. The validity of the perturbation method has been examined and it has 
been found that it introduces at the most a decrease in the theoretically ex- 
pected value by a factor of 1/2. 

3. The change in the radial functions due to deviations from the inverse 
square law is found to cause an appreciable effect on the coupling of the 
valence electron to the nuclear spin for S,2 and P;» terms of single electron 
spectra. This change is not large enough however to account for the difference 
between the ratio of the splittings of the isotopes 199, 201 of Hg observed for 
different terms by Schiiler and Keyston and by Schiiler and Jones. 

The disagreement between theory and experiment presented in §3 ap- 
pears to be universal. Thus in the spectrum of Pb II Schiiler and Jones’ have 
found a displacement between the 206, 208 isotopes of 0.052 cm! for the 
lines 6660, 5609 (7 2S,2—7 ?Pi2, 7 ?Si2—7 *P3/2). By formulas (13), (14) it is 
estimated that if this shift should be due to a change in the nuclear radius it 
would have to be attributed to the 7 °S;,2 term. The value of y?(0) to be used 
for 7S Hg II should be approximately the same as that for 7S Tl II. The 
latter is given by Racah"® as 0.4410" cm~*. The expected displacement is 
thus approximately the same as for 6S of Tl I and is obviously too large. All 
of the shifts could of course be decreased by supposing that yo and v change 
just in the correct proportion. Such an explanation appears to be a forced 
one at present, however. 

The situation is quite different with the 6s76d configuration of Pb II. The 
°—D term is inverted. The /fs of lines terminating on it has been observed by 
Kopfermann,’ Schiiler and Jones,> and in especial detail by Rose and Gran- 
ath? who subsequently to the publication of their results have measured the 
structure of the higher series members of 4245 v7s., AA2948, 2526 (6 7D: 2—5, 
6, 7°Fz/2) also \A5372, 3451, 2887 (6s6p? *D5,2—5, 6, 7 ?Fz2) as well as AA3786, 
2719 (6s6p? ?D3,2—5, 6, ?F 52), 45367 (6s6p? ?Ds,2—5 *F5,2) and 44242 (6°D;,2 
—5 °F;.2). According to their data as well as that of Kopfermann and 
Schiiler and Jones the shift between Pb*®* and Pb’ is practically the same 
for all lines ending on the same D term so that the shifts must be attributed 
to the D rather than the F terms. The displacement for 


6s°6d 2Ds)2 is 0.20 cm™! 6s6p??Ds,2 «0.28 cm" 


6s°6d 2Ds 2 0.09 cm! or (0.08) 6s6p? "D322 0.43 cm™!. 


An application of Eq. (14) shows that for the 6s*6d 7D; terms one expects 
an entirely negligible shift (10-cm™') as a result of a 1 percent change in the 
nuclear radius. Still smaller displacements are expected for 7); 2. The op- 
posite is true experimentally. It is also remarkable that 6s6p* gives larger 
isotope shifts than 6s°6d. It may be however that the 6s*6d terms are per- 
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turbed'? by 6s6p?. If this is the case the isotope displacements of all four 
terms could be attributed to the /ack of one 6s electron in the 6s6p? configura- 
tion. This view is in agreement with the fact that for the 6660, 5609 lines the 
Pb? component is displaced towards shorter wave-lengths with respect to 
Pb*"S while for 4245 and the other lines ending on the 7) terms the Pb*’ com- 
ponent lies towards longer wave-lengths. The direction of the displacement 
is correctly explained by an increase in the radius for the heavier isotope since 
the heavier isotope has the looser binding to the nucleus for the 6s as well 
as the 7s electron. Furthermore the sum of the isotope displacements of 
6s*6d7D; 2 and 6s6p? 7D; 2 is 0.20+0.28 =0.48 cm~! while for the two 7D, » 
the sum is 0.52 cm~'. This sum may be taken to be the isotope shift of a single 
6s electron. It is approximately 10 times as great as the shift for the 7s elec- 
tron (0.052 cm~'). It appears reasonable to suppose that the ratio of ¥*(0) 
for these two electrons is also approximately 10 since Racah’s numerical 
calculations gave a ratio 6.4 for of the 6s and 7s electrons of TI IT. 

Thus even though the absolute value of the isotope shift is smaller than 
would be expected on the hypothesis that it is due to a change of the nuclear 
radius by (2/3) percent (for 206 and 208) the comparison of the 6s and 7s 
displacements indicates that the shift is roughly proportional to the value of 
¥*(0) at the nucleus that the shifts for p, d, electrons are small and that 
therefore the shift may be due to deviations of the field from the inverse 
square law. 

In the hfs of the spectrum of Tl I as it has been analyzed by Schiiler and 
Keyston! and by Jackson,” even this appears untrue. The main argument 
against such an explanation lies in the fact that according to Schiiler and 
Keyston! AA6550, 5528, 5110 (7 27S; 2—8, 9, 10 *P3 2) and ANG714, 5584 (7 7S; 2 
—8,9?P, 2) show no isotope shift. For this reason the 7 *.S; 2 term is supposed 
to have no shift and therefore the shift of 0.05 cm~! present in the 5351 line 
is attributed to a displacement of ?P3.2. It should be noted, however, that 
while Schiiler and Keyston' have found no isotope shift for 3776 (6*P; » 
—7 *S,2) Jackson? reports a displacement of roughly the same amount as for 
5351. If, therefore, the experimental difficulties present in the observation of 
6550, 5528, 5110, 6714, 5584 have masked the shift in these lines one could 
suppose that the 6 *P terms have a negligible isotope displacement while 
7 *S of TP has an energy higher than that of 7 ?S of TP by 0.05 cm~'. This 
assignment would not fit in however with the sense of the displacement in 
TI II as observed by Schiiler and Keyston! because the energy difference be- 
tween the 6s7p and 6s7s configurations is higher for 205 than 203. In the 
TI II spectrum the 7s electron appears to be more tightly bound to the nucleus 
in the 205 than in the 203 isotope while in Tl I even forcing the P terms to 
have no displacements the opposite is the case. Thus this spectrum appears 
not to fit in with a theory of isotopic shifts due to deviations of the nuclear 
field from the inverse square law even in a roughly qualitative manner. 

There is another possibility for accounting for the isotope displacement 


® Shenstone and Russell, Phys. Rev. 39, 415 (1932). Note remark on influence of pertur- 
bations on isotope shifts in Hg. 
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as a direct effect of the mass." It is a striking feature of the observed isotope 
displacements, however, that the odd isotopes of Hg and Pb fall not exactly 
between the even isotopes of closest atomic weight but somewhat closer to 
the even isotope of smaller atomic weight. This displacement cannot have 
anything to do with the mass and it is not due to observational error inas- 
much as in the case of Pb it has been found by three independent sets of 
observers.* +” 

The fact that the observed shifts for S terms are much smaller than those 
expected on the supposition of a progressive change in the nuclear radius 
shows that either the field used in the theory of radioactive disintegration 
cannot be applied in the discussion of the interaction of extranuclear elec- 
trons with the nucleus or that the change in the radius of nucleii of the same 
element is smaller than has been supposed. However, even with such modifica- 
tions the shifts in Tl remain quite unexplained. 

Note in Proof: Since the paper was written we learned of a new nuclear 
theory of Heisenberg according to which the nucleus is supposed to consist 
of protons and neutrons. The Gamow mountain may then be understood as 
the result of exchange forces. We have also learned of the new considerations 
of Bohr according to which the collision between neutrons and electrons is 
very improbable. According to these views the electron may be regarded as 
subject only to the influence of protons within the nucleus. The density of 
protons may be imagined either to be approximately uniform or else to vary 
as a simple function of the distance from the center. It is clear that the term 
displacement is smaller on this view than according to the model used above. 
A calculation for the isotope shift was made using a potential energy 
V= | —(n+1)/n+(1/n)(r ro)" |\(Ze? ro) inside the nuclear radius 7. The dis- 
placement of the energy for s terms is 

dr 1+ p (a + 1) yo? 
é6wW = —a,*y¥?(0)—_—_—_—__ ———- —_ ' 
Z (2p + 1) 2p(2p + 2+ 1)(2p + 1) 


= oI 


and for/]/>0 


2j'(i'- +1 yr? 
sv= Ae Se GO te —} [(/+1)?—a? }i—1-—(P?—a®)!} 8. 
I?(2p+1) 2p(2p+n+1)(2p+1) 








The expected differences in term displacements for isotopes are given there- 
fore by the same expressions as before except for the factor 


A(n) = (n + 1)/(2p + 1)(2p + n + 1). 


For n=2 (uniform charge density in the nucleus) and p=0.81, A =0.25 and 
for n=1/2, A =0.18. Using A =1/5 the isotope shift of Pb*"* and Pb*"* would 
be estimated to be 0.46 cm™ for the 7s and 2.9 cm for the 6s electrons while 
the observed values are 0.052 cm! and 0.50 cm~ respectively. In order to 
bring about agreement between theory and experiment one would have to 
use nm =—0.84, i.e., the nuclear charge would have to be concentrated practi- 


1 TD. S. Hughes and C. Eckart, Phys. Rev. 36, 694 (1930). 
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cally entirely at the center. It is probably incorrect however to give such a 
literal interpretation to the nuclear charge because forces between neutrons 
and electrons may not be entirely negligible and because a-particles may form 
groups having a volume which is affected only slightly by increasing the 
atomic weight. It must be recalled that in the cases considered here one iso- 
tope can be obtained from another one simply by the addition of neutrons. 
The small magnitude of the isotope shift suggests that the volume occupied 
by the positive charges is not affected greatly by the addition of neutrons to 
the nucleus as though the repulsive forces between the positive charges were 
balanced by the cementing influence of neutrons in their immediate neigh- 
borhood. 

Our thanks are due to Professor W. Heisenberg for informing us of the 
new developments before publication and for a helpful discussion of the 
interaction of electrons with the nucleus. 
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Isotope Separations in the Infrared Absorption Spectrum of HCl 
and the Possible Existence of a Chlorine Isotope of Mass 39 


By J. D. Harpy* anp G. B. B. M. SUTHERLAND** 


Physics Department, University of Michigan 


(Received June 30, 1932) 


The absorption band of HCI at 1.74 has been reéxamined under high dispersion 
in an attempt to verify the work of Hettner and Béhme on the existence of a third 
chlorine isotope, Cl’. No trace of the isotope was observed although a path length 
of three meters, some ten times the path used by Hettner and Béhme, was used. The 
isotope separations due to Cl* and Cl*? have been measured very accurately and do 
not agree with the separations calculated by the simple theory. When the anhar- 
monic terms in the energy are taken into account the agreement between the calcu- 
lated and observed isotope separations is within the experimental error. The isotope 
separation for HCI** molecule has been calculated and found to be slightly different 
from that observed by Hettner and Béhme. Our results indicate that the relative 
abundance of Cl** to Cl? cannot be greater than 1 to 500. 


ECENTLY Becker! has predicted the existence of a chlorine isotope, 
Cl’, from a consideration of the form of the infrared absorption curves 

of HCi. This prediction was made from the curves published by Meyer and 
Levin? and by Becker' himself. The absorption maxima in these curves which 
Becker attributes to HCI*® are so small as to be probably within the error of 
the experiment. In fact, Meyer and Levin state specifically that the uneven- 
ness of their base line was due to errors in galvanometer readings, and if 
one examines the original curves of Meyer and Levin it is seen that the max- 
ima mentioned by Becker do not fall at equal distances from the maxima due 
to Cl and in many cases do not show up at all. However, the existence of 
Cl8*° was later apparently confirmed by Hettner and Béhme* working on the 
R branch of the HCI harmonic. In this work the maxima due to HCI** stand 
out in all clearness with the expected constant frequency difference between 
these maxima and those due to Cl**. The present work was undertaken with 
two objects in view, to verify the work of Hettner and Béhme and to measure 
more accurately, under high dispersion, the separation of the maxima due to 
HCI and HCl’. The interest in the third isotope arose naturally from the 
apparent disagreement of observers about its existence. Aston‘ at one time 
considered the possibility of a Cl** isotope but finally decided against its 
existence. Ashley and Jenkins® studied the absorption of AgCl in a region 


* National Research Fellow. 

** Commonwealth Fund Fellow. 

'H. Becker, Zeits. f. Physik 59, 583 and 601 (1930). 

2 C. F. Meyer and A. Levin, Phys. Rev. 34, 44 (1929). 

* G. Hettner and J. Béhme, Zeits. f. Physik 72, 95 (1931). 
‘F. W. Aston, Isotopes, p. 68, 2nd ed., 1924. 

5’ M. Ashley and J. A. Jenkins, Phys. Rev. 37, 1712 (1931). 
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favorable for showing the isotope but found no trace of an AgCl® molecule. 
EK. IX. Plyler (private communication) made an attempt to find the absorption 
lines due to the third isotope by a close examination of the fundamental band 
at 3.4u. He used the same instrument that Meyer and Levin used in their 
work and his results were negative within the error of his observations. The 
interest in the isotope separation arose from a discrepancy in the observed 
values and those calculated on the simple theory. 

With a new type spectrometer® whose resolution and sensitivity are 
greater than any of those yet applied to the problem, the work of Hettner 
and Bihme, Becker, and Meyer and Levin has been repeated. That is, the 
short wave-length side of the harmonic band at 1.74 has been carefully 
examined. A description of the spectrometer has already appeared and it 
need only be mentioned here that by altering the dimensions of the thermo- 
couple a 50 percent gain in sensitivity has been realized over previous ar- 
rangements. 

The region to be investigated, between 1.714 and 1.75, is a particularly 
favorable one in which to work. It is practically free from water vapor ab- 
sorption and in it lies the energy peak of the Nernst heater. Further, the 
impurity of the spectrum due to overlapping orders is very slight even with- 
out the use of any filtering device. The first order energy of 0.85 is small and 
in the second order the energy will be only 10 percent of that in the first order. 
The difference between curves taken with a filter and without a filter is 
negligible. The filter used was opaque to visible light, transmitted 9 percent 
from 0.8u to 0.94, and had 70 percent transmission at 1.73y. 

The grating used in this experiment was ruled by Wood with 15,000 lines 
to the inch and had a ruled surface of 6X 3.5 inches. Previously the emission 
lines of helium had been examined with this grating and the definition was 
found to be good enough to resolve lines at 10829.1A and 10830.3A in the 
second order. The spectrometer was calibrated by means of the helium lines 
at 10830.3A and 20582.1A. The width of the slits for the most detailed exami- 
nation was 0.075 mm and for the rougher curve 0.15 mm. The wave-length 
range across the slit in the first instance was 0.97A or 0.32 cm~'!, and in the 
rougher curve 0.64 cm~!. Readings were taken one slit-width apart, and each 
point represents usually a single observation. The steadiness of the instru- 
ment was such that an observation could be repeated to within 2 percent in 
the detailed curves, where the deflections were about 30 cm between lines, 
and to within about 1 percent in the rougher curve where the deflections were 
about 60 cm. The deflections at places of maximum absorption were practi- 
cally zero for the strong lines, the residual deflection probably being due to a 
slight impurity of the spectrum coming from overlapping orders. No radical 
departures from the base line were observed anywhere except in regions of 
the absorption lines of HCI* and HCI”. 

The present work was done with two cell lengths. The first work, with 
very narrow slits, was done with a cell 31 cm long. This is the length of cell 
used by previous workers and it was found that the strong lines due to 


6 J. D. Hardy, Phys. Rev. 38, 2162 (1931). 
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HCI gave nearly 100 percent absorption with this cell length. When this 
length of cell failed to show the lines due to HCI*® a cell of three meters 
length was used, and it was felt that had the isotope been present in the HCI 
gas which was examined, it would not have escaped detection, assuming that 
its relative abundance to Cl*® is about that shown by Hettner and Béhme's 
curves. 

The 31 cm cell was a brass tube provided with strong windows of plate 
glass. The cell was evacuated and filled with HCI gas from a liquid air trap 
so as to ensure dryness of the gas. The HCI was prepared in the usual manner 
by allowing sulfuric acid to combine with pure NaCl in an evacuated Kipps 
apparatus. The HCI was frozen out by liquid air and then introduced into the 
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Fig. 1. Two sets of observations: R branch of HCl harmonic with 30 cm cell. 


evacuated cell to a pressure of 90 cm of mercury. The trap from which the 
HCl was evaporated was never allowed to come up to the temperature of 
0°C during this process. The three meter cell was of galvanized iron thickly 
coated inside with paraffin. The cell was eight inches in diameter and pro- 
vided with stout mica windows. The HCI was introduced into the cell from a 
liquid air trap and the air gradually displaced by continuous flowing of the 
gas for several hours. Tests at the end of the filling showed that the gas inside 
was about 92 percent HCI. 

Readings were made by taking a complete curve with the empty cell and 
then a run with the cell full. The background with the empty cell was found 
to be very smooth and since lines and line positions were desired, percentage 
absorption was plotted for only one curve, and deflections with the cell full 
for the detailed curves. 
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The results of observations with the 31 cm cell are shown in Fig. 1. Lines 
are drawn to indicate the places where the third isotope might be expected 
from theory, account being taken of the anharmonic terms in the energy. In 
Table I are given the wave-lengths and frequency separations measured in 
this experiment as compared with the values obtained by Meyer and Levin. 


TABLE I. Wave-length and frequency separations in absorption spectrum of HCl. 


Wave-length (A) Isotope separation 
Line Hardy and Meyer and Hardy and Sutherland 
Sutherland Levin (A) (em~') 
2 17519.8 17520.2 [2.3 4+.04+0.08 
3 17467 .7 17467 .9 12.2 4.03 
+ 17418 .3 17419.2 12.2 4.03 
5 17373 .8 17374.9 12.7 4.21 
6 17334.6 17334.0 12.7 4.21 
7 17298 .3 17298 .1 12.4 4.14 
8 17266.3 17265 .2 12.7 4.24 
9 17238 .1 12.6 4.23+0.09 
10 17213.3 13.0 4.34 
11 17193 .2 12.8 4.29 
12 17176.4 
13 17163 .4 


! 
| 
| 
} 





The lines 9, 10, 11, 12, and 13 were observed with the three meter cell. The 
isotope separation for the last two lines is not given since the convergence of 
the band is so strong that the HCl? component of one line overlaps the 
HCIl* component of its neighbor. In Fig. 2 are shown the results obtained 
with the three meter cell. Again no trace of the third isotope is visible. The 
strong lines toward the center of the band are very wide because of the very 
intense absorption in this region. 

Hettner and Béhme state that the theoretical separation between the 
lines due to HCI and HCl? is AX=13.54A, whereas the value found by 
Meyer and Levin is AN=12.0+1A. Hettner and Béhme attribute the dis- 
crepancy to the lack of good separation in the curves of Meyer and Levin, 
and thus to a resulting “attraction” of the stronger line for the weaker. In the 
present experiment the lines are well resolved, not overlapping at all, and the 
resulting separation is in good agreement with the value found by Meyer and 
Levin. Upon closer examination it is found that the separation expressed in 
wave numbers is not constant but steadily decreases as one goes toward 
longer wave-lengths. This variation in isotope separation from line to line 
illustrates the well-known rotational isotope effect. The agreement between 
the accurate theory and the most precise measurements is worth examining 
in some detail. 

We start from the Fues’ expression for the energy of a rotating vibrating 
dipole. 


E =A + hy(n + 3)[1 + BY + 3)?] 
+ (h?/8r*J)(7 + 4)2[1 — K2(j + 4)2] — (A2/8r2J) (mn + 3)°C. (1) 


7 E. Fues, Ann. d. Physik 80, 367 (1926). 











wy 


~ 
= 
S 
i 
w 
x 
aR, 
Rk) 
ny 
<3) 
S 
i 
> 
A 
sy 


“sel OH 40J Suorztsod pajsadxa Zurmoys d1uow sey [DP] JO youRIq y uO 


[99 Ja}9U 9dIY} YUM SUOI}BAIISGGY “7 “B14 

















jQuceqe yweway 








476 J.D. WARDY AND G. B. B. M. SUTHERLAND 


Where w is the vibrational quantum number; j is the rotational quantum 
number; 


A = h dey J; B= (3 2)AR1 + 203); C = 3 + 15Ca + (15/2)C3? + 3C, 


and C; and C; are constants in the potential energy expression which measure 
the degree of anharmonicity of the motion. From this we can immediately 
write down the expression for the frequencies in a vibration-rotation band 
corresponding to the jump of ” from 0 to An. 
vy = vy An — fvoBAn — (hk 8r7J)C(An? + An) 
+ j(— mB — voBAn + h 4n°J — hK* 82*J) (2) 
— j*yyBAn — Ph &SwJ)4AR*. 
Here we have arranged that j is the number of the line from the center of the 
band, positive for the positive branch, and negative for the negative branch. 
Colby’ has shown that this formula is in excellent agreement with the experi- 
mental work of Meyer and Levin on the fundamental and overtone bands of 
HCI. 
As usual let us take p as the ratio of the square roots of the reduced masses 
of the two molecules, then 
Vv, — vo = vyAn(1 — p) — fropBAn(1 — p*) — (hh 8x7J)C(An? + An)(1 — p?*) 
+ 7) —mBA+An)(1—p*)+h 42 (1—p2)—(h 82°) K2(1—p*)} (3) 
— j*vyBAn(l — p*) — jh / 277) KA7(1 — p*) 
Where »; corresponds to the HCl molecule, v2 corresponds to the HCI? 
molecule, 


po = |u( HCI) ‘u( HCI) ]!2 
and the dependence of the various factors on p can easily be verified. The 
values of the constants taken from Colby’s analysis of the bands show that 
the terms with the coefficient )?/ z°J are quite negligible. 
Taking out (1—p) asa factor this may be written 
v, — vo = (1—p) } voAn — lypBAn(1 + p+ pp?) — (h/ 8r2J)C(An? + An)(1 + p) 
+ j[— mB + An)(1 + p + p?) + (h/4x°J)(1 + p)| (4) 
— jBAn(1 + p+ p’)}. 
We may simplify this further by putting p=1 everywhere except in the first 


factor. So the frequency separation for the general band when » changes from 
0 to An is 


v) — vo = (1 — p)yn — (A/8rV)C(An? + An) — 3noBAn 
+ j[h/2x2J — 3Byo(1 + An)| — 37?BroAn}. 


8 W.F. Colby, Phys. Rev. 34, 53 (1929). 
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Where v.=vopAn— jyoBAn—(h/8r?J) C(An?+An) denotes the center of the 
band. The values of the constants deduced in Colby’s® paper for the funda- 
mental band of HCI*® are C=4.885, Byy=0.303, h/42°J = Kv) =21.17. Sub- 


stituting those we obtain for the separations: 


Fundamental 


vy — ve = 0.000751} 2885.88 — 103.3 + 40.57 — 0.91)?}. (6) 


llarmonic 


V} os Vo 


0).000751 } 5667.96 — 309.9 + 39.67 — 1.82j*}. (7) 


Tables II and III show the agreement between the values of the separations 


as determined experimentally and as calculated from these expressions. 


TABLE II. Jsotope separations in the fundamental band in cm. 


Positive branch 


j Calculated Meyer and 
value Levin's value 
1 2.32 2.09 
2 2.14 2.09 
3 247 2.10 
4 2.20 2.16 
5 2.23 2.22 
6 2.25 2.32 
7 2.20 2.13 
8 2.29 2.08 
9 2.31 2.26 
10 2.33 3.43 
11 2.34 2.25 
2 2.36 2.34 


| 


Negative branch 


Calculated Meyer and 
value Levin's value 
2.06 2.10 
2.03 1.97 
1.99 1.98 
1.96 1.90 
1.92 2.02 
1.88 1.72 
1.84 1.74 
1.80 Lae 
1.76 1.83 
1.72 1.61 
1.67 Sues 
1 1.57 


.63 


TABLE III. Jsotepe separations in the harmonic band—positive branch. 


Calculated Observed 


Hardy and Hettner and Meyer and 
Sutherland Béhme Levin 
I (em~) (A) (cm~) (A) (A) (A) (cm~) 
2 4.08 7.5 4.04 2.2 12.5 12.5 4.20 
3 4.10 2.5 4.03 12.2 12.0 ob 3.87 
4 4.12 iZ.o 4.03 12.2 12.5 ‘2.2 4.03 
5 4.14 iZ.s 4.21 2.7 12.0 12.0 4.19 
6 4.15 12.5 4.21 12.7 12.0 4.22 
7 4.16 12.45 4.14 12.4 12.0 4.14 
8 4.17 12.4 4.24 i2.7 4.26 
Mean 4.13 12.5 4.12 2.3 12.3 4.13 


The discrepancies between individual line 


separations as computed and 


as observed lie in each case within the experimental error. The agreement of 


the various values for the mean separation in the harmonic band is very 
satisfactory. The variation in the separation from one end of the funda- 
mental band to the other illustrates particularly well the rotational isotope 


etiect. The result of considering the anharmonic terms is to reduce the value 
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of the separation from that calculated on the simple theory which would be 
given by Eq. (5) with B=C=0. The anharmonic term involving C is the 
most effective in reducing the separation. The interaction terms involving B 
begin to contribute appreciably only when j is very large. 

If now we use the above formulae (with an appropriate value of p) to 
compute the isotope separations for an HCl*® molecule we obtain a set of 
values which disagree slightly with those found experimentally by Hettner 
and Béhme. The comparison is shown in Table IV. 


TABLE IV. Jsotope separations in the harmonic band of HCI”. 


j Calculated value (A) Hettner and Béhme's value (A) 
2 23.8 24.5 24.5 
3 23.8 24.5 25.5 
4 23.8 24.5 24.5 
5 23.8 23.5 23.3 
6 23.8 2369 
7 23.7 23.5 
Mean 23.8 24.5 


The difference is small but the close agreement of calculated and observed 
separations in the case of the Cl’ isotope make it doubtful whether these 
lines are due to HCI. 

To obtain an estimate of the smallest relative abundance of HCI®® to 
HCl} which we could detect, we may compare the absorption coefficient of 
the weakest line of HCI** which we could detect with that of the strongest, 
vis., line 3. This can be done theoretically with formulae due to Dunham.?® 
Now when a is the absorption coefficient so computed for a certain line it is 
well known that!? 


a= KN 


where N is the number of molecules in the initial energy level and K is an 
abbreviation for a number of factors which would not change appreciably 
for an HCI or an HCI*® molecule. If then we know the least ratio of a’s 
which we can detect, we know the least ratio of N’s which we can detect. We 
employ an equation of Dunham’s, modified for use in the positive branch. 
a3 Ps v3 { 2 ‘ ; a 

— {1 + 2a(vm — v3) }e(Bma- Be) 47 
Am Pm Vm 


where m denotes the number of the line in the positive branch; p is the mean 
of the statistical weights of the initial and final states of the molecule; 2a is 
a small correction factor which we can calculate from Dunham's paper to be 
very close to 0.0154; E denotes the energy of the molecule in the mth rota- 
tional level of the zero vibrational level. When we apply this equation to cal- 
culate a;/ai3 using energy and frequency values from Colby’s paper® we ob- 


§ J. L. Dunham, Phys. Rev. 34, 438 (1929). 
10D. M. Dennison, Phys. Rev. 31, 503 (1928). 
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tain a ratio of 460 to 1. The calculated ratio of a3/ai, was 1600 to 1. Using the 
three meter cell we could easily detect line 13 and there is an indication on 
the curve of what may be line 14. The latter is, however, uncertain, but we 
can safely state that had the Cl** isotope been present in a ratio to Cl® of 1 
to 500 we could not have failed to detect it." 

There is still another viewpoint from which our inability to substantiate 
the work of Hettner and Béhme seems remarkable. The instrument we used 
was of considerably higher resolving power than that employed by them. Our 
lines are consequently deeper and have higher maxima than theirs and so our 
instrument is better adapted for the detection of faint companion lines. When 
we remember that the path length used in the present work was ten times that 
employed by Hettner and Béhme our failure to detect any sensible change in 
absorption at the expected places is quite irreconcilable with their con- 
clusions. 

It is a pleasure to express our thanks to Professor D. M. Dennison for the 
suggestion that the anharmonic terms might be important in the calculation 
of isotope separations and for his interest in the problem. 


11 Note added in proof: Since the above article was written the experiment has been re- 
peated using an absorption cell seven meters long. There were still no indications of lines which 
might be interpreted as coming from an HCI** molecule. Line 14 of the positive branch was 
easily detected giving quite a distinct absorption peak on the curve. This enables us to place 
the abundance ratio of Cl** to Cl** lower than 1 to 1600, 
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The Spectra of Na II, III, and IV in the Extreme Ultraviolet 
By Bruce B. VANCE 
Ryerson Physical Laboratory, University of Chicago 
(Received July 7, 1932) 


New lines of Na II, III and IV have been found and classified in the region between 
434,000 and 242,000 cm 7. 


INTRODUCTION 


HE spectrum of Na II in the near ultraviolet has been classified by 

Frisch! in a recent article, who also gives a list of references to previous 
work on this spectrum. In the extreme ultraviolet region Bowen? and Edlen 
and Ericson* have observed the resonance lines 372.065A and 376.350A2 The 
spectrum of Na III (fluorine-like) has been studied by Edlen and Ericson* 
and Mack and Sawyer.® The lines corresponding to the 2s—2p transitions 
are the only lines recorded by these investigators for this ion. Both Edlen and 
Ericson* and Mack and Sawyer’ predict a PP® group of Na IV in the region 
of 244,000 cm~!. This group had been found and identified previously by Gale 
and Hoag of this laboratory, but had not been published. The observation of 
this group and of other lines in this region, which might be attributed to so- 
dium, caused Professor Gale to suggest the present problem to the author. 


(GENERAL PROCEDURE 


The source was a high potential vacuum spark, (hot spark).® This took 
place between hollow aluminum electrodes stuffed with metallic sodium 
which had first been rolled in nickel gauze. It was found that when the 
pressure in the discharge chamber was below 10~‘ cm of Hg and the electrodes 
were separated about 1.5 mm, the vacuum gap (source) would not break down 
below 40,000 volts. The maximum voltage of the high potential source was 
above 100,000 volts. The primary circuit of the transformer was automatic- 
ally made and broken by a very quick acting switch twenty-two times each 
minute, the duration of the sparks being from one-twentieth to one-tenth 
seconds each. The total time required for an exposure was from two to four 
hours. 

The vacuum spectrograph used in this investigation was the one designed 
by Hoag,’ with a few changes. The grating was glass, instead of speculum, 
with 1500 lines to the inch. The plate holder accommodated 2 X 18 inch plates 

1 R. Frisch, Zeits. f. Physik 70, 498 (1931). 

21. S. Bowen, Phys. Rev. 23, 1 (1924). 

3 B. Edlen and A. Ericson, Comptes Rendus 190, 116 (1930). 

4B. Edlen and A. Ericson, Comptes Rendus 190, 173 (1930). 

5 J. E. Mack and R. A. Sawyer, Phys. Rev. 35, 299 (1930). 

61. S. Bowen, J.0.S.A. 13, 89 (1927). 

7B. Hoag, Astrophys. J. 66, 225 (1927). 
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and was so arranged that it could be raised to allow three exposures on one 
plate without removing it from the vacuum. Very thin glass Schumann 
plates, made by Hilger and Company, were used. 

In reducing the plates a linear dispersion was assumed between the Al IV 
line of 160.073A% and the O III line at 599.598A,° error curves were then 
plotted by the aid of intermediate standards (O III 266.992," Na II 372.065,' 
Na III 380.0914 and O III 525.795A).° These error curves were found to be 
extremely smooth for all plates. 
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Fig. 1. Typical spectrogram. 


With a very narrow slit and the pressure in the discharge chamber below 
10-4 cm of Hg (as observed with an ionization gauge) plates were easily ob- 
tained which were free from fog and weak gas lines. Under these conditions 


8 J. Soderqvist and B. Edlen, Zeits. f. Physik 69, 356 (1931). 

9 From a list of standard oxygen wave-lengths furnished through the courtesy of Dr. 
Bengt Edlen of Uppsala. 

10 A. Ericson and B. Edlen, Zeits. f. Physik 59, 659 (1930). 
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only the prominent oxygen, carbon and aluminum lines appeared as impuri- 
ties, Which were used as standards in reducing the plates. 


Na II (NEON-LIKE) 


In the region investigated one would expect to find, besides the promi- 
nent resonence transitions s*p* ‘Sy—s*p*3s °P), ‘Pi, transitions from two or 
three of the three 3d levels with J=1 and possibly from higher s orbits than 
the 3s. 

It was found possible to form two series formulas of the type 


y= py —4tR (m+ u+ am)? where m = 1, 2,3,--- 


whose limits differed by the same amount (in cm~!) as the separation of the 
ground terms of Na III (s*p° *P; 2,32). In addition to the three pairs of tran- 
sitions used in forming these series formulas a fourth pair and possibly a 
fifth pair were found. (The latter was later found to fit better into the term 
scheme of Na III.) The levels thus found correspond to the ms: and nis, 
levels of neon as given by Paschen,'! whose limits are 7P°;» and *P%s », re- 
spectively. From the value of the latter limit 47.068 volts has been obtained 
as the ionization potential of Na Il. The lines resulting from transitions from 
these levels to the ground term s*p* Sy) can be seen easily on the plates (Fig. 1 
typical). They are given in Table I and II and are shown diagrammatically 
in Fig. 2a. 


TABLE I. Observed and calculated transitions of Na IT. 


I \ (obs.) vy (em~) Config. and Terms v (cale.t) 
6 376.350° 265710 stp 'Sy—stp'3s 8P,°(s,) 265710 
10 372 .065° 268770 Spi IS)—sp3s 'P°(s2) 268770 
0 302.28 330830 stp Sy —s*pts (sy) 330945 
3 301.61 331540 s*p* 1S) —s*p'3d Zi) 331502! 
4s (s2) 332811 
2 300.38 332910 Spt 1Sy—s*p* 
3d 2. 333050! 
1 282 .96 353200 Spi IS,—s*p’5s (s,) 353283 
1 281.81 354850 SpPIS,—sp 5s (s2) 354854 
0 275.10 363500 Sp IS)—s*pOs (s,) 363510 
0 273.99 364960 s*p* 1S)—s*p's (sy) 304888 
(0 271.01 369000 stp 'Sy—stp’7s (sy) 369028)* 
(0 269.98 370400 sp’ 'S,)—s*p7s (s) 370459)* 


* Probably Na III. 

t Calculated v obtained from the following formulas: 

msy + + + p=382800 —4R/(m +0 .963845 —0 .001872/m)? 

ms, + + + v=381430—4R/(m +0 .949650 —0 .002057 /m)* 
in which m=1, 2, 3, °° > 


Of the three 3d levels with J =1 one should have the limit s*p° ?P°, » and 
the other two the limit s*p*° *P°; 2, of the latter two only the transition to the 
ground term from one of them is observed. This level is calculated to be 
331,540 cm~! above the ground term, and corresponds, within experimental 
error, to the Z;, term of Frisch.' The other 3d (J =1, ?P°; 2 limit) falls too close 
to the 4s(s2) level to be resolved, it would correspond to the Zz term of Frisch. 


11 Paschen-Goetz, Seriengesetze (1922). 
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The transition from this level to the ground term is evident only from the 
fact that the line classified as the \Sj)—4s(s2) transition is unusually wide for 
a line of that intensity. The 4s(sz) and 4s(sy) levels correspond to the Z; and 
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TABLE II. Term values of Na II. 
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(c) 


Fig. 2. Term schemes for Na II, Na III, and Na IV for transitions in the far ultraviolet. 


Z\4 levels of Frisch. The other 4s and 3d terms, which have been calculated 
from Frisch’s values, are indicated by dotted lines in Fig. 2(a) and are given 
the following J values; Z1.5.8.10.15 J=2, Zs.9.12 J=3, Zi3 J=4 and Z; J=0. 
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The levels in Fig. 2(a) between 335,000 and 330,000 cm~! are Z; to Zig reading 
down, Zs.9 are indicated as one level. 
Na III (FLUORINE-LIKE) 

In addition to the previously observed transitions 2s—2p reported by 
Edlen and Ericson* and Mack and Sawyer,’ lines have been found which can 
be classified as s*p*—s*p'3s transitions. These are given in Table III, shown 
diagrammatically in Fig. 2(b) and the term values are given in Table IV. 


TABLE III. Observed transitions of Na IIT. 


I N vy (cm~) Config. and Terms 

6 380.0914 263095 Ss? p® *P°ij2—sp® *Sy2 

8 378.1203 264466 S2p? 2P°s.2—sp® 2Si2 

0 272.59 366840 Sp 2*P°s2.—s*pt3s *P52 

0* 271.01 369000 sp 2P°s—S ps *P 3/2 

Q* 269.98 370400 Sp 2P®s.—Spi3s *Pr2 

2 208.81 372060 Sp 2*Po2—s*pi3s 2 P32 

6 207.78 373430 s*p° 2P oso —S2pi3s 2 P 3/2 
374600* Sp 2P®sj.—Spi3s 2P ij 

2 252.26 396410 Sp 2P® yy. —Spi3s 2Dsi2 

5 251.25 398010 Sp 2P°s.—s2p'3s 2Dsi2.3/2 

1 231.03 432840 Sp 2P®eij.—S*pt3s 2Syys 

2 230.25 434210 Sp 2P°s,,—s*pi3s 28 


° Possibly Na Il. 

t Masked by O III line 266.992A. 

The configuration s*p'ns theoretically yields inverted *P)2.325 2, 7Pi2, 
terms with limit at s*p* *P; *D3 05,2 terms with limit at s*p* 'D2; and a 


Si term with limit at s*p* ‘So. Lines are found which are considered to 


so > 


TABLE IV. Term values of Na ITI. 


Config. Term Term Value Limit 
(cm!) 

s*p? 2P 3/0 0 

s*p° 2P 1/2 1370 

s pS 2Si/2 264400 

Spi3s $P 5/9 366840 “Pp 

Spi3s §P 2 369000 7. 

Spi3s 1P yj. 370400 ‘Pp 

Spi3s *P 3/2 373430 <a 

Spi3s *Pre *374000? a 

Spi3s *Ds/2 398010 1p 

Sp'3s 2Ds/2 . 398230 1p 

Spi3s 2Siy2 


434210 1s 


t Masked by O III line 2606.992A. 


come from all these levels, for 7=3, except from the *Pj,2 level, and these 
levels fit well into a Moseley diagram with the same terms of other members 
of this isoelectronic sequence (F I'? and Ne II"). 

The transitions from the ?P/2 level to the ?P°3,2,1;2 levels are assumed to 
be masked by the prominent O III line at 266.992A and the line correspond- 
ing to the transition ?P°;,2—?P3)2 of Na III respectively. The relative position 

2 H. Dingle, Proc. Roy. Soc. A117, 411 (1927). 

8 T. L. deBruin and C. J. Bakker, Zeits. f. Physik 69, 19 (1931). 
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of this *P1/2 level on this assumption is indicated by a dotted line in Fig. 2(b). 

The lines here identified as corresponding to the transitions between the 
*P1/2,3,2 and the *P°s/2 levels possibly may be the fifth members of the ms2 and 
ms, series of Na II, but it seems more likely that they belong to Na III. 


Na IV (OXYGEN-LIKE) 


The PP® group at 244,000 cm~! that had been found and identified as 
that of Na IV on photographs taken by Gale and Hoag,“ during their in- 
vestigation of multiply ionized lithium, was obtained on nearly all the plates 
taken during this investigation. The central line was not resolved, which 
made it difficult to determine exactly the position of the *P; and *P»° levels. 

In addition to the triplet levels there should be a ‘Dz and a 'Sy level for the 
configuration s*p', and a 'P, level for the configuration sp*. Lines have been 
found which probably correspond to the allowed transitions between these 
levels, and these transitions fit well into a Moseley diagram with the same 
transitions of Ne III given by Boyce and Compton." The actual position 
of these levels can not be given from the results obtained, but it is possible 
to place them roughly from the information obtained from the positions of 
the 'P, *P, *D, and 2S levels of Na III. This has been done in Fig. 2(c). 
TABLE V. Observed transitions of Na IV. 


Config. and Terms 








I d vy (cm~) 

3 412.46 242448 stp! 3P, — sp 8P°, 
P 411.25 243161 spt §P,—sp® *P°, 
6 410.41 243654 stpt 3P, ,—sp>®P% , 
2 409 .46 244224 Spt §P, —sp® §P% 
3 408 .56 244762 Sp! §P,—sp *P°, 
0 351.00 284900 spt 1S)—sp> 'P%, 
3 


319.61 312880 stp' 'D.—sp 'Py 





TABLE VI. Term values of Na IV. 





Config. Terms Term Values (cm!) 
s*p* 3P, 0 
stp 3p, 1110 
stp 3P, 1600 
S*p* 1D» A 
Spt 1S A +27980 
sp 3p°, 244050 
sp’ 3P°, 244760 
sp> 3P°, 245330 


sp 1P®, A+312880 


In conclusion the author wishes to express his appreciation to Dean Gale 
for suggesting the problem and for his helpful encouragement throughout the 
investigation, and to Professor R. S. Mulliken for his advice in the analysis. 

Note added in proof: J. Soderqvest, Zeits. f. Physik 76, 316 (1932), gives 
an analysis of the Na III spectrum which agrees well with the results given 
above. 


4H. G. Gale and B. Hoag, Phys. Rev. 37, 1703 (1931). 
% J. C. Boyce and K. T. Compton, Proc. Nat. Acad. Sci. 15, 656 (1929). 
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The Absorption Spectrum of Benzene in the Region of 6.75 


By S. SILVERMAN 
The Johns Hopkins University 


(Received July 5, 1932) 


The absorption band of benzene lying at 6.754 has been investigated under high 
resolution for both liquid and gaseous states. The spectrum for the vapor shows three 
distinct branches with evidence of unresolved fine structure. The liquid shows a less 
regular spectrum with a definite shift towards longer wave-lengths. In addition, sonte 
new water vapor lines in the 6.264 band are reported. 


N RECENT years an increasing amount of work has been done investigat- 

ing the infrared absorption and Raman spectra of organic substances. The 
success which followed the unravelling of the fine structure of some of the 
more simple compounds led to attempts to resolve the spectra of molecules 
of increasing complexity. Among the group offering promise due to its com- 
plete symmetry is benzene, which has been investigated by a number of 
workers.!? 4.4.5.6 Particularly in the last two or three years, considerable work 
has been done in the near infrared in this laboratory by Barnes® and Leber- 
knight.’* Continued improvements in apparatus and resolving power, con- 
sisting chiefly of the employment of the Pfund anastigmatic spectrometer 
with echelette gratings of high resolution, have led to discovery of new bands, 
and increased complexity in some of the older ones. These studies have thus 
far been restricted primarily to the region in the neighborhood of 3u. Bands 
which previously were so broad as to permit, at best, mere designation as 
combinations of various vibrational frequencies began to show some evidence 
of vibration-rotation structure. As a result of several discussions with Dr. 
C. E. Leberknight, the writer felt that it might be profitable to examine in 
great detail one of the bands somewhat farther out. The band at 6.75u, which 
is quite intense and hence probably represents a fundamental vibration, of- 
fered excellent possibilities. 

For the purpose of this investigation, the grating spectrometer employing 
the resonance radiometer, which was used in a recent mapping of the ab- 
sorption and reflection spectrum of calcite,* was found readily adaptable. 
The grating had 1364 lines per inch; the slit-width was 20A. As in this pre- 
vious experiment, a magnesium-oxide filter was used to remove the higher 


* In this work, which is soon to be published, the spectrum of benzene was investigated 
in solid, liquid and gaseous states in the region of 34 under high resolution. 

1 Dougherty, Phys. Rev. 34, 1549 (1929). 

2 J. Ellis, Phys. Rev. 27, 298 (1926). 

3K. W. F. Kohlrausch, Der Smekal-Raman Effekt (Springer 1931). 

4 Schaefer and Matossi, Das Ultrarote Spektrum. P. 270 ff. 

®’ R. B. Barnes, Phys. Rev. 35, 1524 (1930). 

6 R. B. Barnes, Phys. Rev. 36, 296 (1930). 

7C. E. Leberknight, unpublished results. 

8S. Silverman, Phys. Rev. 39, 72 (1932). 
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orders of shorter wave-lengths. The instrument was calibrated by means of 
visual observation of higher orders of the 5461 mercury line, and also by 
prominent atmospheric absorption lines.’ The vapor cell was prepared by 
sealing polished rock-salt windows to the ends of a brass tube. The cement 
used was a clear lacquer which was found to be almost completely unattacked 
by the benzene. The optimum cell length was found to be about 22 mm, which 
length of vapor absorbed 50 percent of the incident energy in the neighbor- 
hood of the absorption maximum. A free surface of liquid benzene was kept 
continually in the bottom of the cell in order to maintain constantly satura- 
tion pressure. No difficulty due to condensation on the windows was exper- 
ienced. A dummy cell was prepared to give a comparison set of readings, as 
all measurements were taken cell-in and cell-out. The liquid cell consisted of 
two flat plates of rock salt, separated by thin strips of silver-foil 0.01 mm 
thick. This thickness of liquid gives the equivalent molecular path length of 
22 mm of gas; this is desirable, as it is essential for direct comparison of gas 
and liquid to have very much the same molecular thickness to traverse. The 
cell was made tight on 3 sides with lacquer, and the 4th, open, side was fed 
occasionally with liquid from a pipette; an identical dummy cell was also 
prepared. The cells maintained their performance unaltered throughout the 
duration of the experiment. 


RESULTS 


The absorption in the vapor cell shows the 6.754 band to be clearly 
divided into three distinct bands lying at 6.683; 6.743; 6.790u respectively ; 
the band at 6.743 being somewhat broader and slightly more intense than 


TABLE I. Absorption maxima for benzene. 














Vapor v Liquid y 
r | r 
6.683 +0 .0015. 1496.3+0.4 6.6434 1505 .3+0.4 
6.683 1496.3 
6.743 1483.0 6.753 1480.8 
6.790 1472.8 6.830 1468.1 








its companions. All three of the branches show clear evidence of unresolved 
fine structure, with small but distinct maxima separated by intervals as 
small as 2 wave numbers. It is entirely probable, however, that the true fine 
structure of the band is such as to defy resolution by the instrument at our 
disposal; in which case these small maxima would have no meaning other 
than to indicate the presence of other fine structure. Inasmuch as the benzene 
is so complex and has such a large moment of inertia (in the neighborhood of 
11X10-*°g X cm?) one would expect the structure to be very complicated and 
the lines to be densely packed. 

The absorption of the liquid presented a less attractive map. The central 
band from the vapor absorption is sharpened slightly with a shift from 6.743 
to 6.753u. The band at 6.683 has been so altered as to bear little resemblence 


® Plyler and Sleator, Phys. Rev. 37, 1493 (1931). 
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to the gaseous state. It is much weaker, being nearly washed out, and shows 
two maxima: one at 6.643 and one at 6.683u. The 3rd band has also changed 
so as to present an entirely new picture. It is much broader than its counter- 
part in the vapor state, and it is quite irregular in outline, with its maximum 
at 6.830u. Its intensity is about that of the corresponding vapor band. 
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It is perhaps worthy of note that in the plotting of this spectrum many 
runs of the atmospheric absorption bands were obtained. An admirable check 
of many lines found recently by Plyler and Sleator® was obtained. In addition 
some new ones, marked by the above authors as missing, were observed. 
These are listed in Table IT. 


TABLE II. New absorption lines in water vapor. 


DN v r v 


6.284 +0.0015.z 1591.3+0.4 6.367 +0.00154 1570.6+0.4 
6.294 1588.8 6.480 1543.2 
6.299 1587.6 6.485 1542.0 


1519.1 





6.313 1584.0 6.583 





In conclusion, it might be remarked that at present it appears impossible 
to derive a detailed description of the quantum states of the benzene molecule 
from the data at hand. The present investigation has shown the existence of 
hitherto undiscovered fine structure, and it is planned further to probe the 
spectrum with instruments of still higher resolution in an attempt to unravel 
the fine detail. 

The writer takes great pleasure in expressing his thanks to Dr. C. E. 
Leberknight for many fruitful discussions concerning this experiment. He 
also wishes to acknowledge a Grant-in-aid of the National Research council 
which has helped make this investigation possible. 
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The Theory of the Faraday Effect in Molecules 


By ROBERT SERBER 
University of Wisconsin 
(Received June 27, 1932) 
In treating the Faraday effect two cases may be distinguished, depending upon 


whether the frequency of the incident light is near resonance or well removed from 
resonance with absorption lines of the molecule. 


Frequency of incident light well removed from resonance with any absorption 
lines. In this case it is imperative to include the perturbation of the intensities by the 
magnetic field, as well as the perturbation of the energies. A general expression is ob- 
tained for the rotation by molecules (poly-, di-, or monatomic), of the form 


: | vA(nn’) v?B(nn’) v?C(nn’) 


(n'n)?—v?)?—v(n'n)?—-v? T(v(n'n)? 3) i : 








where V is the Verdet constant. This formula contains dia- and paramagnetic terms 
of the usual type, but augmented by terms arising from perturbation of the intensi- 
ties. It contains, in addition, other diamagnetic terms which have the same frequency 
dependence as the paramagnetic terms. For atoms this expression reduces to that 
given by Rosenfeld. However for diatomic molecules our results differ from Kronig’s, 
since we include the effects of the components of magnetic moment perpendicular to 
the axis of figure. The terms arising in this way were omitted by Kronig, although 
generally they are of the same order of magnitude as the contribution of the parallel 
component of the moment. 

Independence of spin. When the over-all spin-multiplet width is small compared 
to kT/h the rotation is completely independent of spin. As a consequence the para- 
magnetic terms vanish for nonlinear polyatomic molecules, and for linear polyatomic 
and diatomic molecules in & states. 

Magnitude of the rotation and comparison with experiment. The classical Becquerel 
formula for the Verdet constant is V=y(e/2mc?)vdn/dv, with y =1. It is shown that 
the rotations, in the visible and near ultraviolet, of the gases for which data are 
available should be approximately representable by a formula of this form, provided 
7 is given the proper value. The y value should lie between zero and one. This con- 
clusion, in all cases but one, agrees with the known facts. The exception is oxygen, 
but the data are probably in error because of polymerization effects which seriously 
alter the absorption at high pressures. 

Frequency of incident light near resonance with an absorption line. Here only 
the perturbation of the energies by the magnetic field need be considered. It is shown 
that the rotation in iodine vapour, observed by Wood, is due to rotational distortion 
of the excited “IIo level, which partially uncouples the spin moment from the figure 
axis. The magnetic rotation spectrum of the alkalis, observed by Wood and Loomis, 


is also explained, in particular the quenching of lines of large rotational quantum 
number. 


WO different experimental procedures have been employed in measuring 
the Faraday rotation in gaseous molecules. One method uses a contin- 
uous band of incident frequencies, and a low gas pressure; rotation is then 
observed only in the immediate neighborhood of an absorption line. The other 
uses monochromatic incident radiation, and high pressures; it serves to mea- 
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sure the rotation in “transparent” regions of the spectrum. The first method is 
exemplified by Loomis’s measurements of the magnetic rotation spectrum of 
the alkalis, the second by Siertsema’s work on the rotation of hydrogen in the 
visible. 

Correspondingly, the theory must distinguish two cases: 

(I) Incident light near resonance with an absorption line of the molecule. 

(II) Incident light well removed from resonance with any absorption line. 


Anticipating our later results, we may remark that (1) is the simpler of 
the two, as here the rotation is caused by the Zeeman effect alone. In case 
(II) it is essential to include, in addition, the effect of the perturbation of the 
amplitudes (i.e., intensities) by the magnetic field. 

$$1 to 5 of this paper deal with case (II). §1 discusses, in a non-mathe- 
matical way, the theoretical attack, the general conclusions, and the experi- 
mental data. $§2, 3, and 4 give the mathematical development of the theory. 
$5 is concerned with the magnitude of the rotation and with the application 
of the theory to a number of simple examples. 

$6 treats case (1), in particular, the rotation of iodine as observed by 
Wood, and the work of Loomis mentioned above. 


$1. Discussion OF THEORY AND EXPERIMENT 


The theory of the Faraday effect in atoms has been treated by Rosenfeld! 
and others. Rosenfeld found, for the case of multiplet widths small compared 
to k7T/h, that when the frequency of the incident light is not near any ab- 
sorption line the rotation can be expressed as a sum of terms of two types: 
the diamagnetic, characterized by independence of temperature, and the 


r 


paramagnetic, proportional to 1/7. The diamagnetic terms are given by 


dre Nv? v(n'n) | P(nn’) | 2 


‘ —— (1) 
3hme?  ~— (v(n’n)* — v?)? 


where V, the Verdet constant, is the rotation per gauss per cm, the P(nn’) 
are the matrix elements of the resultant electric moment, and vy is the fre- 
quency of the incident light. Eq. (1) can also be written in the Becquerel form 


: e On 
V =—py—> 
2mc? av 
in which x is the index of refraction, as given by the Kramers dispersion 
formula.” This value of V will be called the normal Verdet constant. 
The purpose of the present paper will be to extend the theory to molecules. 
A start in this direction has already been made by Kronig,’ who, under the 
same condition that v is not too near the frequency of any line, finds a Verdet 


'L. Rosenfeld, Zeits. f. Physik 57, 835 (1930). 

2 See, for example, J. H. Van Vleck, The Theory of Electric and Magnetic Susceptibilities, 
p. 362. 

5 R. de L. Kronig, Zeits. f. Physik 45, 458 and 508 (1927). 
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constant of only one-third the normal value.‘ This result seems unsatisfac- 
tory, for it does not agree with Rosenfeld’s calculation in the limiting case of 
a central field. It will be shown that the error lies in Kronig’s neglect of the 
components of magnetic moment perpendicular to the axis of figure of the 
diatomic molecule. Since in a central field the components of magnetic mo- 
ment along any three orthogonal axes contribute equally to the rotation, 
neglect of two of these components will naturally lead to Kronig’s result. 

We shall suppose throughout that v is well outside the Zeeman pattern of 
any line, i.e., that y—v(n'n) is large compared to the Zeeman displacements. 
When only this restriction is made the rotation of any molecule’ is found to 
be given by a formula of the form 





pee Lon'n)? — vp)? v(n’n)? — vy T(v(n'n)? — vy? 





j v°a(nn’) v*b(nn’) +. v?c(nn’) \. WalkT, (2) 
) 


with B=1/>_,e-"»/*?, The three kinds of terms appearing in Eq. (2) will be 
called “a terms,” “b terms,” and “c terms,” respectively. The a terms are the 
contributions of the Zeeman effect; the c terms are due to the altered distri- 
bution of molecules in the various normal states caused by the magnetic field. 
The 6 terms arise from the perturbation of the amplitude elements of the 
electric moment by the magnetic field. This perturbation is, of course, a 
consequence of the existence of matrix elements of the magnetic moment be- 
tween the various states. 


Definition of the term ‘‘magnetic interaction” 


For brevity, when an element of the magnetic moment between two states 
exists, we shall speak of the magnetic interaction between the states. 


Coalescence of b terms into a and c terms 


, 


The magnetic interaction between a state 7’ and a state n’ 
tion proportional to the difference quotient® 


gives a rota- 


1 { 1 1 \ 
ee ee (3) 
v(n'n’”) \ vy? — p(n'n)? sv? — v(n'’n)? 
If v(n’n’’)/(v—v(n'n)) is small the difference quotient can be expanded in 


powers of this ratio (since v(n’n’’) =v(n'n) —v(n''n)) giving 


2v(n'n) ! 2v(n'n)v(n'n"’) 4 \ 


(v? — v(n'n)?)? 








v? — v(n’n)? 


* Kronig has calculated the rotation explicitly only for diatomic molecules in '= states. 
In comparing Kronig’s result (p. 511, reference 3) with our Eq. (1), the factor 2 in his expression 
must be absorbed by summing over A= +1 rather than only over A=1. Making the substitu- 
tion | ao(lo)|?=4! Pe(lo)!2, and remembering that in a central field | Pe(JJ:)|?=4| P(loh)|?, 
one readily verifies that the two expressions differ by a factor 1/3. 

5 Throughout this paper the term molecule is to be understood as including atom. 

6 It must be cautioned that this picture is incomplete, as not all b terms appear in the 
form of difference quotients (see Eq. (8)). 
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Neglecting terms after the first we see that the two } terms contained in 
Eq. (3) have reduced to a single a term. Similarly, b terms due to magnetic 
interactions between states with slightly different Boltzmann factors give 
rise to c terms when expanded in powers of hy(n'n’’)/kT. Of course 6 terms 
caused by interactions between states of widely different energies cannot be 
simplified in this way. 

The net result of the operations just described is the enrichment of a and 
c terms at the expense of ) terms. These modified terms will be called A, 
B, and C terms to distinguish them from the a, 6, and c terms. It should be 
noted that the A and C terms can not be ascribed any such simple origin as 
was possible for the a and c terms. Hereafter, whenever we refer to Eq. (2) we 
shall mean the equation modified in this way, rather than the equation as 
originally given. 


Elimination of Boltzmann factors 


The exponential factors drop out of Eq. (2) when the following conditions 
are satisfied: 

a. The over-all spin-multiplet width is small compared to k7/h; or 

a’. The separation between the two lowest spin-multiplet components is 
large compared to k7/h. 

b. There are no electronic states with separations from the lowest state 
of the order of R7/h. 

c. The separations between the rotational levels of the normal state are 
small compared to k7/h. 

d. The molecule can be considered rigid (e.g., rotational distortion neg- 
lected). 

e. The states of the molecule can be divided into groups.’ A set of states 
M1, M2, - + * M, Will be said to constitute a group if v(,n,) is small compared 
to both y—v(n,p) and v(n,m), n, and n, representing any two members of the 
set, p one of the normal states, and m any state, normal or excited, not be- 
longing to the set. Thus in Fig. 1 the states g form a group. 

It is evident that condition (e) restricts us to case (II) discussed in the 
introduction, for if y—v(n,p) is too small (i.e., if the incident light is “near 
resonance with any absorption line”) division of states into groups will be im- 
possible. When vr is restricted to such a range that (e) is satisfied, the allow- 
able spread of a group will be widely different in different parts of the spec- 
trum. For example, groups lying near v or near the normal level may consist 
of states differing only in their rotational quantum numbers, while groups 
far from v and the normal level may include different vibrational states, or 
even states arising from different electronic configurations. 

When either set of conditions, (a), (b), (c), (d), (e), or (a’), (b), (c), (d), 
(e), is met, the rotation is given by a formula of the form 


y v?A(nn’) v?B(nn’) v?C(nn’) \ 


(v(n’n)? — v?)? — v(n’n)? — v®? = T(v(n’n)? — v*)} 





(4) 


7 In dividing the states into groups not all states need be considered, but only those having 
nonvanishing terms in Eq. (8). 
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with A(nn’), B(nn’), C(nn') constants independent of temperature. Both A 
and B terms are thus diamagnetic, the C terms paramagnetic. One must not 
be misled by the terminology into supposing that the signs of paramagnetic 
and diamagnetic terms are necessarily opposite, for this is not the case. The 
presence of diamagnetic B terms, having the same frequency dependence as 
the paramagnetic terms, has been generally overlooked in the previous 
theories (it was noted, however, by Van Vleck®), although, as will be shown 
later, the B terms are actually no less important than the more customary A 
terms. 
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Fig. 1. Fig. 2. 


In virtue of condition (e) and in view of the discussion of Eq. (3), we can 
now say that the A terms are due to magnetic interactions between states 
lying in the same group, while B terms are due to interactions between states 
in different groups. 

The A terms of Eq. (4) are of a more general type than the A terms ob- 
tained by Kronig and Rosenfeld, since ours include the interactions between 
various states in a group. These interactions were neglected by Kronig® and 
did not exist in the case considered by Rosenfeld. Rosenfeld’s case of atoms 
with multiplet widths small compared to k7/h is simplified by the fact that 
My,and M, can be considered constants of the motion,!® and so the z com- 


8 J. H. Van Vleck, reference 2, p. 370. 

® Kronig, to be sure, included the elements of u, between states of different j, but not be- 
tween states of different electronic quantum numbers. It is to the latter sort that we refer. 

10 Strictly speaking, this is true only if the magnetic field is strong enough to produce the 
Paschen-Back effect; however the distinction is trivial in view of the invariance of the spurs 
appearing in Eq. (9). 
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ponent of the magnetic moment has no nondiagonal elements. Thus there 
are no interaction terms like Eq. (3), with the consequence that the A terms 
reduce simply to the normal Verdet constant, while the B terms vanish en- 
tirely (see Eq. (9)).'! Kronig’s calculation of the A terms for diatomic mole- 
cules is valid in the limiting case of an axial field so far from central in charac- 
ter that the components of the magnetic moment perpendicular to the figure 
axis have elements only between states in different groups. Here obviously 
the A terms are due only to the parallel component of the moment, as Kronig 
assumed. The B terms will not vanish in this case, and, in fact, are of the same 
order as the A terms, although omitted by Kronig. Our A terms agree with 
Kronig’s and Rosenfeld’s in these limiting cases, but otherwise they are not 
the same. For example, for nonlinear polyatomic molecules the mean value 
of the magnetic moment along any avis is zero;'* hence Kronig’s method 
would give no rotation at all, although we obtain both A and B terms. 


Independence of spin 


An important characteristic of Eq. (2) (and consequently of Eq. (4)) is 
that, when (a), (c) and (e) are satistied, it is completely independent of spin. 
In any calculation of the rotation, subject of course to these conditions, the 
spin can be regarded as completely uncoupled from the orbital angular mo- 
mentum, and, in effect, forgotten completely. This result is analogous to the 
well-known theorem that the magnetic susceptibility is independent of the 
Paschen-Back effect, and like it is a consequence of the fact that we sum our 
expressions over all Zeeman components. In view of this property of the rota- 
tion the attempts to explain the anomalous behavior of oxygen as due to 
its spin are not justifiable. Another conclusion which can immediately be 
drawn is that the paramagnetic terms of Eq. (4) vanish for nonlinear poly- 
atomic molecules, and for linear polyatomic molecules and diatomic mole- 
cules in Y states. This is because the paramagnetic terms are due to the mag- 
netic interactions between the normal states. Since spin need not be consi- 
dered, these interactions are clearly zero in the cases indicated. 

The absence of a spin-paramegnetic rotation is in sharp contrast to the 
presence of a spin term in the magnetic susceptibility under the same sup- 
position of narrow multiplets. The difference arises because, under the cir- 
cumstances, the spin and orbital moments can be regarded as completely 
independent systems. They thus will contribute additively to the suscepti- 
bility; on the other hand, the spin appears in the rotation only in cross terms 
with the orbit, and such cross terms average to zero. 


Magnitude of the rotation 


It frequently happens that the rotation has not been measured much be- 
yond the visible, while all the v(m’n) are in the far ultraviolet. The A and B 


" Considering the phases proper to a central field, it is readily seen that, in this case, 
Eq. (9) is identical with Rosenfeld’s expression. 

2 This statement might not hold if the molecule had a very high degree of symmetry; 
however, whenever we consider nonlinear polyatomic molecules we shall assume that this is 
not the case. 
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terms are then practically indistinguishable, the net result of the B terms being 
to change the apparent values of the A(nn’). Thus the diamagnetic terms 
of Eq. (4) will have, effectively, the frequency dependence, v?/(v(n'n)? —v?)?, 
of the Becquerel formula. In consequence, as will be shown more fully in §5, 
the rotation can often be represented fairly well by the Becquerel formula 
with an anomalous value of e/m. 

In general our knowledge of the molecule is not sufficient to permit ex- 
plicit calculation of all the terms required by Eq. (4), but simple considera- 
tions lead to the conclusion that the Verdet constant should ordinarily be less 
than the normal value. The experimental results in all cases confirm this. 
The reason is roughly as follows: The rotation is largely due to the nearest in- 
tense group of lines in the ultraviolet, which are heavily weighted compared to 
those further out. If all the lines lay close to this nearest group a normal 
Verdet constant would result, as is shown in $5. The effect of removing some 
of the lines further into the ultraviolet, where in fact they belong, is to 
diminish greatly their contribution to the rotation. If the field is not central 
this diminution is greater than the corresponding decrease in 0n/dv, and so 
leads to a smaller Verdet constant. 


Discussion of experiment'* 


Siertsema’s measurements of the rotation in hydrogen in the visible are 
representable almost exactly by the normal Verdet constant. This result, 
remarkable in its agreement with the classical theory, will be discussed in §5. 
However Sirk’s measurements in the ultraviolet show some evidence of the 
presence of a B term. The magnitude of the observed rotation in nitrogen is 
63 percent of the normal Verdet constant, i.e., the Becquerel formula is 
valid if we insert a factor 0.63. As the measurements have not been carried 
bevond the visible no differentiation between A and B terms is possible. The 
case of N.O is similar. Carbon dioxide has been observed in the ultraviolet, 
but not far enough to give any conclusive evidence as to the presence or ab- 
sence of B terms (see $5). 

Oxygen, on the other hand, is definitely anomalous. The measurements of 
Siertsema and Sirks in the visible and ultraviolet cannot be fitted by Eq. (4) 
using the observed absorption frequencies, and, as we have seen, this cannot 
be blamed on the spin. However the experimental data are open to serious 
objection because of the high gas pressures used: Siertsema, working in the 
visible, used a pressure of 100 atmospheres; Sirks, 80 atmospheres in the 
visible, 40 in the ultraviolet below 2800A. At these pressures a considerable 
amount of oxygen would be in the form of the polymer Q,.'** The equilibrium 
constant is somewhat uncertain; at 80 atmospheres the estimates of the 


8 An excellent account of the experimental data on the Faraday effect has been given by 
C. G. Darwin and W. H. Watson, Proc. Roy. Soc. 114A, 474 (1927). The reader is referred to 
this paper for references to the literature. 

a The existence of this polymer was kindly called to my attention by Professor R. Laden- 
burg. 
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partial pressure of O, range from 1 to 10 percent of the total pressure." 
According to Wulf," O,; shows continuous absorption below 2400A, as well as 
a number of bands extending towards longer wave-lengths. Moreover, the 
source of light used by Sirks and Siertsema was a quartz mercury arc, and 
Warburg® has shown that at high pressures the 2537A line is efficacious in 
producing ozone. Ozone shows strong absorption in the infrared, and in the 
ultraviolet beginning at about 2800A. Sirks found, in his 230 cm column of 
oxygen at 80 atmospheres, that absorption began at 2805A; in fact it was 
for this reason that he was compelled to reduce the pressure for his shorter 
wave-length measurements. At 40 atmospheres the absorption had moved 
down to 2654A. QO, shows no absorption in this region, but both O,; and Os; 
bands have been reported at 2810A and 2654A;"° these might readily account 
for Sirk’s results. 

In view of these facts we believe the existing measurements of the rotation 
entirely unreliable; solution of the problem must await more unambiguous 
data. 

Monatomic gases should, according to Rosenfeld, have a normal Verdet 
constant. This has been verified for neon, although at only one wave-length, 
5460A. The same result has also been found for argon.’ 

In general the rotations in the large number of liquids and solids cited by 
Darwin are not accurately representable by the Becquerel formula, even with 
an anomalous value of e/m. There are two obvious reasons for this: first, the 
A(nn’) are not proportional to the f values of the Kramers dispersion formula; 
second, B terms are also present. 


$2. DERIVATION OF THE GENERAL FORMULA FOR THE ROTATION 
Let the magnetic field and the light ray be directed along the s axis. The 
rotation of the plane of polarization per centimeter is’ 


6 = CB Dor(n'n)[P(nn’) Py(n'n) Jem ¥" 7, (5) 


mr 


where 


B=WN doe Walk?) C = 4r°v7i/ch, 


n 


and P,(nn’), P,(nn’) are the elements of the x and y components of the 
electric moment in the presence of the magnetic field. Here, and throughout 
this paper, the square brackets have the meaning 


| P..(nn’) P,(n'n) | 


P(nn’)P,(n'n) — P,(nn')P(n'n); 
and 
t(n'n) = 1/(v? — v(n'n)?). 


"O.R. Wulf, Proc. Nat. Acad. 14, 356, 609, 614 (1928); Jour. Am. Chem. Soc. (II) 50, 
2596 (1928). 

s References to Warburg’s work are given by Wulf, reference 14. 

® Wulf, reference 14; D. Chalonge and M. Lambrey, Comptes Rendus 184, 1165 (1927). 

17 R. de Mallemann, L. Gabiano, and F. Suhner; Comptes Rendus 194, 861 (1932). 

's H. A. Kramers, Kon. Acad. Wet. Amsterdam 33, 959 (1930). Kramer’s Eq. (34) should 
be multiplied by 47, as can be seen from his Eq. (32). 
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We expand P,, P,, v(n'n), W,, B in power series of the form P, = P,°+ HP, 
+--+; then, under the supposition that »y—v(n'n) is large compared to the 
Zeeman splitting, the terms of Eq. (5) linear in 3 give 


V = 0/K = CB df °(n'n)([P.%(nn’) Py (n'n) | + [P.(nn’) > (n'n) |) 


+ 2v(n'n)v™ (n'n)79(n'n)?| P(nn’) P,(n'n) | (6) 
— (W,/kT)7°(n'n) [P2(nn’) P,(n'n) || eWa AT, 
P‘»(nn’) is expressed in terms of the unperturbed amplitudes by means of'® 


1 » (Han) P (n"'n’) P(nn")H® (n"'n’) 
P,(nn’) ~~ pe — Smee me <eceeesnarersimentttne \. (7) 


y"(nn"’) y*(n''n') 


"i 


The prime on the summation sign indicates that n” =n and n” =n’ are to be 
omitted from the summation. Denoting the z component of the magnetic 


moment by u., we have J] (nn’) = —yu.(nn’). Using these, and a similar ex- 
pression for P,(nn’), we obtain 
, + ‘ ‘ , a f f 
V = — (CB*/h) ; v(n'n")— | P(nn')p(n'n"”) Py (nn) | 
nn’n’’ 


f _ > 
r°(n'n) — 7°(n'’n) eo Wl kT 


— (2CB"/h) do v°(n'n)79(n'n)?[P"(nn’)y.(n'n'’) P,"(n'n) |e Walk? 


nn’ 


+ (CB°/h) ‘y(n’ n)—79(n'n) 1 [P.(nn’) P(n'n”) |u(n’n) 
' 


nn’n’’ 
+ p(nn”’)|P,(n''n’) Py (n'n) | erWe kT 
+ (2CB°/h) Sov%(n'n)r%(n'n)2y.(nn) [P2(nn’) Py(n'n) je AT 
po 7 9 , ) > 0 , P i} , ¥ Wa" kT 
+ (CB°VRT) r"(n'n)p-(nn)|P,"(nn') Py (nn) je . 
Here we have extended our definition of the square bracket to include 


7 


[P.°(nn’)y(n'n") Py(n''n) | = P22(nn’)p(n'n”’) PY (nn) 


— P,(nn’)p(n'n”")P (nn). 

In accordance with our division of states into groups, all matrix ele- 
ments will now be written (dropping henceforth the superscript") in the form 
P,(nk;n'k'), where n and n’ represent the quantum numbers which are the 
same for all states in a group, k and k’ the quantum numbers which differ- 
entiate the various states within a group. 

The terms of the first sum in Eq. (8) for which »” =n’ (although, of 
course, k” #k’) can be simplified by expanding the factor in curly brackets in 
powers of v(n’'k’;n'k"”)/(v?—v(n’k' ;nk)*?). To simplify the terms of the third 
sum for which 2” =n (k” #k) the method of pairing terms previously used by 
Van Vleck?" is employed. 


1” See, for example, Born and Jordan, Elementare Quantenmechantk, p. 198. 
20 J. H. Van Vleck, reference 2, p. 191; Phys. Rev. 29, 727 (1927). 
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The pair of terms 


can be written 


; | P.(nk; wk’) P (n'k'; nk'’) |u:(nk”; nk) + (nk; nk") [P(nk"; n'k')P,(n'k’; nk) |} 


Putting 


and expanding gives (the first bracket representing the first line of the above expression) 
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(1/h)v(nk"’; nk)~'r(n'k’; nk) } |P(nk; n’k’)P,(n'k’; nk’) Jus(nk’’; nk) 
+ p(nk; nk”) [P(nk”; n'k’)Py(n'k’; nk) |} eek? 

+ (1/h)v(nk; nk") 1(n'k'; nk’) | | P2(nk"; n'k’)P(n'k’; nk) |u(nk; nk’) 
+ u(nk”’; nk)[P.(nk; nk’) P(n'k’; nk’) |} eM nate? 


1 ( e Wik kT e “Wine? kT 


hv(nk”; nk) bv? — v(n’k’s nk)? v? — v(n’k’s nk’)2)- 
Ware = Wak + ho(nk"; nk), 
v(n’k’; nk’) = v(n’k’; nk) — v(nk"’; nk) 





ppp erat, (teint) 
4 Y¢2QkT vp? — v(n’k’; nk)? k7 


Qv(n'h’; nkyewW nx a hv(nk"’; nk} 


kT 


+ h(v? — v(n'k’; nk)?2)2 


plus higher order terms which are very small. Apart from them, this is the same as 
| | Ps(nk; n'k!)Py(n'k’; nk” |u(nk”; nk) + wlnk; nk”) |P.(nk"; n’k’)P(n’k’; nk) |} 
| (1/h)v(n’h’; nk)r(n'k’; nk)? + (1/2kT)1(n’h’; nk) | eM nel kT 
+ || P.(nk”; n’k’)P,(n'k’; nk) Jus(nk; nk’) + u:(nk”; nk)|P.(nk; n’k’)Py(n'k’; nk’) |} 
{(1/h)v(n’h’; nk!)r(n’h’; nk’)? + (1/2kT)r(n'k’; nk’) | eM nw kT, 


The purpose of these manipulations is to remove all small »’s from the denominators of 


Eq. (8). A primed sum is no longer required for the terms written in this form, since the diagonal 
terms (i.e., 


pearing. The expansion of the first sum of Eq. (8) leads to a similar simplification, the diagonal 
terms being contained in the second sum. 


It is now convenient to introduce the tibermatrix notation used, for ex- 
ample, in Born and Jordan’s Elementare Quantenmechanik. In this notation 
P,"" signifies the matrix of all the elements P,(nk;n’k’) having a given n 
and ’, Sp is an abbreviation for spur. 

One finds, on carrying out the operations described above that Eq. (8) 


becomes 


V = (CB/h) Yor(n'n)r(n'n)? Sp \(- 2[P 60m pn) Py’) | 


k” =k) are already present in Eq. (8): namely, the fourth and fifth sums there ap- 


+ [P(e Pr) | yr + pe"™ [P60 Pym) Jy e-W™ kT | 
+ (CB/h) Zz — v(n'n"’)-') r(n'n) — r(n'’n) | 


n’n!? (n!#n!") 

Sp {[Pet9 nina» Pye enw a 

+ v(n'’'n)—'7(n'n) Sp {([Prm Py’? | r™ 
“- pz (nn) | P (nn) P (n'n) |) e-W™ [47 | 


+ (CB/2kT) >o7(n'n) Sp i ([ P26 Py’™ | 


(9) 





+4 ge™™ [P (mn) Pyin'n) |) e-W™ 147 | ; 
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As all small v’s have been removed from the denominators of Eq. (8) we have 
been able, in Eq. (9), to write v(n’k’ ;nk) =v(n'n). In the first sum of Eq. (9) 
small terms involving v(nk’;nk) have been omitted. This corresponds to 
neglecting the pure rotation spectrum. 


$3. INDEPENDENCE OF SPIN 


The spurs appearing in Eq. (9) are invariant of a transformation to a 
system of representation appropriate to strong fields. If the over-all spin- 
multiplet width is small compared to k7/h the matrix e-”“’/"7 is independ- 
ent of spin. Then the “strong field” representation the spin appears only in 
a:, Which breaks up into two parts, w:.+p,z, of which only the second con- 
tains the spin. On summing over m, this term vanishes, since y,, is a diagonal 
matrix with elements ranging from —28s to 26s. Thus, when the over-all 
multiplet width is small and when + is not near resonance with any line, the 
rotation is independent of the spin.*! This result, of course, is valid for all 
field strengths (barring saturation effects). 

In $5, dealing with the magnitude of the rotation, we shall consider only 
narrow multiplets; 4, may accordingly be taken as representing only that 
part of the magnetic moment due to the orbital angular momentum. 


$4. ROTATION IN RIGID MOLECULES 


We will make the approximation of supposing the molecule rigid, that is, 
we suppose the “internal” wave functions, which are expressed in terms of a 
coordinate system rotating with the molecule, independent of 7. Then the 
summation over j, m, j’,m’, 7”, m” implied in Eq. (9) can be readily performed 
by a method given by Niessen.” Let x’, y’, z’ be the rotating system of co- 
ordinates. On applying Niessen’s method to Eq. (9) we find” 


V = (CN/3hg) >) Sov(n'n)r(n'n)? Sp | — 2[ Pe py OO Py | 


z’y’2’ n’ 
+ [Pp (mn’) Py (nn) Ty, (an) + peer (™ [Pym Py n'nd| 
+ (CN /3hg) > > — v(n'n!)-"} r(n'n) — r(n''n) } 


z’y’2! on’n’’(n’/Pin’) 


Sp [ Pepe (mn) (nn) Py (nn) | + v(n'’n)—'7(n'n) 
Sp [Pp mn) Py Cun) |, (nn) oa pee (nn) [ Pa (n'a) Py, (n'n) | | 
+ (CN /6kTg) > dor(n'n) Sp | [Pa Py ’™ | em 


z’y’2’ n 


(10) 


a gee") | P,(mn’) Py n'n) | ' 


The symbol >>,.,.: indicates a summation of the three terms obtained by 
a cyclic interchange of x’, vy’, z’. g is the statistical weight of the normal state 


21 This result has been given, for central fields, by Rosenfeld, reference 1, and in the old 
quantum theory, also for central fields, by Darwin, Proc. Roy. Soc. 112A, 814 (1926). 

2 K. F. Niessen, Phys. Rev. 34, 253 (1929). Niessen’s derivation is valid only for diatomic 
molecules, but the proof can be given for the general case. 

% It should be noted that Eq. (10) holds when either condition (a) or (a’) is satisfied. 
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of the stationary molecule. The elements of matrices such as P,-‘""” are, of 
course, independent of j, m, j’, m’. 

Making a simple rearrangement of the terms in the second sum of Eq. 
(10), we find the explicit expression for the A(mn’), B(nn’), C(nn’) of Eq. (4): 


A(nn’) = (CN /3hgv?) do v(n'n) 
z’y’2’ 


Sp i 2[ Pz’ pane Pp fn'nd | 
oo [Pp (mn’ Py | em pe [Pn Pca | 
(CN /3hgv?) > : 3 — v(n'n'’)-! Sp } [Pg (me Dy, (nn) Py (nn) 


r'y’s’ n’’(n’’#n’) 


B(nn’) 


+ [Pepe (mn’ yy, (nn) p(n’) | + v(n'’n) 1 Sp YP (me Py Cem) |, vm 
a pe (nn) | Pi (nn) Py (nn) | 
Clan’) = (CN /okgo4) Sp [Pe Py ae pe [Po Py] 


a* y's’ 
$5. MAGNITUDE OF THE ROTATION 


For nonlinear polyatomic molecules the matrix u“"” will be zero, there be- 
ing no “precession frequencies” which are not large compared to k7/h. Thus 
the C(mn’) will vanish. Diatomic molecules and linear polyatomic molecules 
in S states will also have no elements in u‘"" different from zero, so the 
C(nn') again vanish. In these cases there will be no paramagnetic rotation. 

The work of Rosenfeld and Kronig give limiting values for the A(nn’). 
We have shown in $1 that in a central field the A(mm’) have just the values 
which give the normal Verdet constant, and in the other limiting case (for 
diatomic molecules) of an extreme axial field the A(nm’) are reduced to one- 
third the previous value (the three terms of the >>,,,.. sum reducing to a 
single term), unless both » and n’ are © states, in which case A (nn’) vanishes 
entirely.** The A (nn’) of polyatomic molecules similarly vanish for sufficiently 
asymmetrical fields. 

The whole question of the magnitude of the rotation is greatly clarified 
by consideration of a few simple examples. The molecular wave function (in 
the x’, y’, 3’ system) will be supposed the product of the hydrogenic wave 
functions of the individual electrons. This makes A a constant of the motion. 

Let the normal state be (mso)?=, as in the hydrogen molecule. Take for 
excited states nson’prIl and nson'’po >, and let 4, v2, vs be respectively the 
nson' poX—(nso)*X, nson'pril—(nso)*X, nson' pox —nson'pril separations. 
This situation is illustrated in Fig. 2 (a). No assumption will be made as to 
the relative magnitudes of »; and ve, other than that condition (e) is satisfied ; 
they must accordingly be considered as belonging to two different groups. 
This procedure, of course, is valid even when »;—P (i.e., v3) is small. Using 
the well-known hydrogenic moment matrices we readily find from Eq. (10), 


*4 We are thinking here of the practically important case, in which the normal state is a = 
state. In general the one-third remaining after passage to the axial field limit will be divided 
among the states of different A into which the central field group breaks up. 
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_ tre Nv? { Pp Pp? veP* \ 
O} Semmetee, ceveminencneneen Say Snmmmeneiionn, ip enaien my" {> (11) 


shmec? Lvs(v2 — yp?) v3(v,? — v?) (vo? — pv? 


where P? =2 Pi (n ‘pr, nso) [? in the notation of footnote (25). 
The normal Verdet constant is 
4reNv? j vy, P? 2v2P? \ 


3hmc? lo — py’)? (v2? — vp): 





The asymptotic properties discussed above are easily verified. If vis 
made small compared to v,;—v, Eqs. (11) and (12) both reduce to 


V = V, = 4reNv?veP?/hmc*(v2? — v*)?; 


thus V is just the normal Verdet constant. If v3 is very large, so that the 
first two terms of Eq. (11) and the first term of Eq. (12) are negligible, the 
rotation, as given by Eq. (11), is reduced to one-third the above value, which 
is Kronig’s result. In this case V=3V, (the index of refraction being reduced 
to 2/3 its central field value, while the rotation is reduced to 1/3). 

If y is defined by y= V/V,, we have seen that y =1 for a central field, 
and y = 0.5 for an extreme axial field.*° In intermediate cases y will lie between 
these values. For example if we take vy; = 212, vy =v2/6 we find y =0.817. Tak- 
ing v=r2/4 gives y= 0.813. If v. corresponds to a line at 1000A these two 
values of v correspond to 6000A and 4000A respectively. This simple calcu- 
lation strikingly illustrates the point made in $1 that A and B terms are, in 
such a frequency range, practically indistinguishable. Further, it shows very 
clearly why a Becquerel formula with an anomalous value of e/m is approxi- 
mately valid for many substances. If we put y=vr2/2.5, corresponding to 
\ = 2500A, we obtain y = 0.79. Thus even measurements this far in the ultra- 
violet would be of little value in distinguishing between Eq. (11) and the 
modified Becquerel formula. 

Magnetic interactions between some of the excited states may be due to 
the “core,” or part of the molecule exclusive of the excited electron, for when 
we remove an electron from a closed shell the residuum may have a magnetic 
moment. Of course states which have different core-configurations are greatly 
different in energy, and so the magnetic interactions of the core may easily 
connect widely separated levels. On the other hand, a normal Verdet con- 
stant requires narrow, rather than wide, separations. The essential difference 


* For example, denoting the normal state by y,, the II state by ¥2, and the wave functions 
of the individual electrons by (| x;) we have 
yi = (nso | x1)(nso| x2), ¥2 = (nso | x;)(n' pr! x2), 
P,(12) _ ef fur* (x1 + x)~dridr, = — ef (nso | x2)*x2(n' pa | x2)dr2 
P,, (nso; n'pr). 


The elements such as P,,(nso; n'pr) are given in all the standard texts, e.g., H. Weyl, The 
Theory of Groups and Quantum Mechanics, p. 200. 

26 If the = state is taken as the lower of the excited states y would range from 1 to 0 (we 
are supposing vy <<», 2). For polyatomic molecules also, y would be expected to range from 1 to 
0. These figures are exclusive of negative terms, which will be discussed later. 
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between the rotations in atoms and molecules lies in the spread, in molecules, 
of levels which coincide in atoms; and the more complicated the core the 
larger this spread is likely to be. One would therefore expect the rough rule 
that molecules with complicated cores have lower values of y than simpler 
molecules. The data confirm this, the y values for Hs, Ne, COs, NeO being 
respectively 1.00, 0.63, 0.56, 0.34. 

The central field behaviour of hydrogen in the visible is now readily ex- 
plained, for the core consists of a single electron in ag state, and thus plays 
no part in determining the rotation. As one would anticipate, the potential 
energy curves”’ show that the excited lso2po¥ and 1sa2pzll states lie quite 
close to each other as compared to their separation from the normal state. 
Since there is little possibility of exciting the core electron to obtain other 
states which could contribute to the rotation, the rotation is evidently due 
almost entirely to the two excited states mentioned above; and, as we have 
seen, these closely spaced levels should give a normal Verdet constant. 

As an example of the case in which the core does have a magnetic mo- 
ment we shall take for the normal state (7 pz)*(npo)*X, the nitrogen configura- 
tion, and for excited states (mpz)*(nfo)*n'soll and (npr)* nton'’soy. Magnetic 
interactions between these states will evidently be due to the core of the II 
state. Proceeding exactly as in the previous example we again find the rota- 
tion to be given by Eq. (11). In a central field the energies of our two excited 
states will coincide, and a normal Verdet constant will result: in such a field 
no distinction need be made between the magnetic moment of the core and 
the magnetic moment of a valence electron. In the axial field of a diatomic 
molecule the energies of the excited states will be widely different, since these 
states have different core configurations; y will thus be considerably less than 
unity. This result, as contrasted with the result obtained for the hydrogen 
molecule, illustrates the effect of a complex core in diminishing the Verdet 
constant. 

When we deal with homonuclear diatomic molecules another factor may 
enter which we have not yet considered. For such molecules P(nn’) has ele- 
ments only between the normal state and states of opposite symmetry with 
respect to reflection of the electrons in the origin. States ’’ of the same sym- 
metry as the normal state will have no effect on the refraction, but Eq. (10) 
shows that they will affect the magnitude of the rotation if the elements 
u.(nn’’) are not negligible, although no terms containing 7t(’’n) will appear. 

A situation in which yu, has elements between the normal states and an ex- 
cited group of states can be obtained by inverting the energy levels of our 
first example, so that the (msa)*= state lies above the nson’prlIl and nson’ paL 
states, as is shown in Fig. 2(b). Eq. (11) then holds only if v3 is small compared 
to k7T/h, so that there is no difference in the concentrations of the molecules 
in the mson'pril and nson'po® states.** The second term of Eq. (11) is then 
due to molecules in the latter state, the first and third to those in the former. 


27 R. S. Mulliken, Rev. Mod. Phys. 4, 1 (1932). 
*8 A weight factor 1/3 should be included in Eq. (11) in this case. There will also be a 
paramagnetic term; however it is of no interest for our purposes. 
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However if the II state is the higher of the two, that is if v3 is positive and 
large compared to k7/h, the first and third terms vanish and we are left with 
a negative B term. A negative B term is also obtained if the = state is the 
higher. 

In our examples the effect of magnetic interactions between two excited 
levels is to give, when vy <7, ve or vy >, v2, two B terms which, although they 
are of opposite sign, give a net rotation which is positive; when ve<v <7, 
two negative B terms. The effect of interaction between the normal state and 
an excited state of frequency v; is to give a negative B term when vy<», a 
positive term when vy>»,;. Negative B terms may account for the negative 
rotation of TiCl,. 

It is of some interest to calculate the rotation under the assumption that 
all excited levels lie in the same group. Suppose the electric moments have 
elements only between the normal states, 7, k, and a group of states, ”’, R’, 
having the same n’. Under this hypothesis the second sum of Eq. (10) 
vanishes, and the teilmatrices appearing in the remaining terms can be re- 
placed by the complete matrices (the spur, of course, is still taken only over 
the normal states). The third sum of Eq. (10), and the second and third 
terms of the first sum, are zero in virtue of the commutability of the coordi- 
nates. On summing over all electrons we have 


P w:Py Pyu:P, 


— (e*/2mc) > xix iPy; — ViPri) Ve — Vel XjPyi — ViPei)Xi 


ijk 


— (e?/2mc) So xin (Py irs — VePyi) + vey (Pej Xi — XiPz;) 


ijk 


II 


— (e*h/4rmci) Yo (xixj + yi) = —(eh/4amci)(P + P,’). 


Adding the terms obtained by cyclic interchange of x, y, 2, we find 


dreNv? y(n’ n) 1 


—_—_—__— > | P(nk; n’k’) | 2, 


3hme? (vy? — v(n'n)*)? go ge 


' 


which is just the normal Verdet constant. 


$6. ROTATION NEAR AN ABSORPTION LINE 


Going back to Eq. (6) we see that when + is close to one of the v(n’n) the 
rotation is 
v(n'n)v(n'm’; nm)|P(nm;n'm’)P,(n'm'’;nm)| 
v= 8 rit, (13) 


bag (v? — pv(n’n)*)? 


for, since this term contains 1/(v?—v(n’n)?)? rather than 1/(v?—v(n'n)?), it 
is very much larger than all the others. The rotation thus depends only on 
the Zeeman effect, the perturbation of the amplitudes being negligible. 











504 ROBERT SERBER 


Rotation in iodine 


The iodine lines observed by Wood** come from a *IIp—'Z transition, and 
so exhibit no Zeeman effect if we use the “rigid molecule” approximation. 
Actually, the Zeeman effect, and consequently the optical rotation, is due to 
rotational distortion, that is, to the effect of molecular rotation in uncoupling 
the spin from the nuclear axis of figure. The suggestion that rotational dis- 
tortion is responsible for the behaviour of iodine was originally made by 
Kemble,*’ but an adequate mechanism was lacking, as at the time the tran- 
sition was mistakenly classified as 'X—'Y. Kemble supposed that the orbital 
moment, rather than the spin moment, was uncoupled by the molecular ro- 
tation, but it now seems much more probable that spin uncoupling is the 
preponderant effect. 

We shall use the rotating system of axes, x’, y’, 2’, previously employed by Van Vleck;*! 
the z component of the spin in fixed axes is given in terms of this rotating system by 

S: = S$, sin @ + s, cos @, 
@ denoting the angle between the axis of figure of the molecule and the fixed axis s. The term 
involving s.’ gives the spin part of the ordinary rigid molecule Zeeman effect. In our case this 
just cancels the effect of the orbital angular momentum, since 2=0 in the “Ilo state; we shall 
not be further interested in this term. The matrix s,’ sin @ can be written 
(s,7 sin 0)(a.jm; o,/j’m’) = 2s,(o.0,')L(o.jm; o.'j'm’). 


The elements L(c.jm; ¢,'j’m') are defined and evaluated in a paper by Rademacher and 
Reiche.® As the multiplet separations of the *Ilg state of iodine are very wide, only the first 
order perturbation of s, by the molecular iotation is required. This is given by a formula similar 
to Eq. (7), from which we find for the diagonal elements of s., when 2=0, o,=1, 


s2(— Lym; — Um) = 2H) (— 10)L(Ojm; — 1jym)s,-(0 — 1)/iv(— 10) 
+ 21(— 1jm; Ojm)s,(— 10)H©(0 — 1)/hv(— 10) 
= 4H‘)(— 10)L(Ojm; — 1jm)s,-(0 — 1)/hv(— 10), 


since all the quantities involved are real. 

It is convenient to write s:™ (o,jm; o,'j’m'), L(o.jm; o,'j'm'), and v(o.0,’) in terms of & 
rather than o,. The elements appearing above are, in this notation, s, (Ojm; Ojm), L(1jm; 
Ojm) and »(01). We then have® 


Sy'(o.5 0, + 1) = 3[S(S + 1) — o,(0, + 1)]"; 
H(o,; 0. + 1) = b[j(j + 1) — O(@ F 1)]"2[S(S + 1) — oo, + 1)"; b = H2/8x2MP?; 
L(1jm; Ojm) = — m/2|j(j + 1)}"?; 


and from these 
sz) (Ojm; Ojm) = — 2bm/hr(01). 


29 R. W. Wood and G. Ribaud, Phil. Mag. 27, 1009 (1914). 

8° E. C. Kemble, Bull. Nat. Res. Council 57: Molecular Spectra in Gases, Chap. VII. 

tJ. H. Van Vleck, Phys. Rev. 33, 467 (1929). The primed and unprimed systems of his 
Eq. (5) will be interchanged, to accord with our previous notation. We shall follow this paper 
in our choice of phases for H and s,’, and in notation: ¢, will denote the component of angular 
momentum along the figure axis, ¢, the component of spin along the figure axis, 2=¢,+e; 

82 H. Rademacher and F. Reiche, Zeits. f. Physik 41, 453 (1927). See Eq. (13) and the table 
on p. 475. 

* Van Vleck, reference 31, Eqs. (30) and (31). From Rademacher and Reiche’s table we 
find L(1jm; 0jm) =| m| /2[j(j+1) |"*. The minus sign must be taken when m>0 to match Van 
Vleck’s phases for H™ (see Van Vleck’s footnote (25)). 
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Using the fact that L(1jm; Ojm) =L(—1jm; Ojm), it is readily verified that the same result is 
obtained when o,;= —1. 
Thus for our “lly —'S transition 


v\OG'm’; jm) = — 48bm'/h*(01), 


8 being the Bohr magneton number. 
The elements of the electric moment appearing in Eq. (13) are the familiar amplitudes of 
the symmetric top, given, for example, by Kronig.*' The terms different from zero are 


iP? (j+m+2)(j7+m+1) 
P.j,m;j+1,m + 1)P,j+1,m + 1;7,m)| = —- — - — ’ 
. : ' ; | 2 (27 + 1)(2j + 3) 

iP? (j= m+2)(j5 1 
Ltt ated xt, ww + lh t,e 2 hhabl« Se eh ), 
. . . 2 (27 + 1)(27 + 3) 
where P?= P,’('S, *Ilp) 2. 


The rotation is found by substituting these expressions in Eq. (13) and 
summing over all values of m, the result being multiplied by two, since the 
two upper states, ¢, = +1, contribute equally to the rotation. In this way we 
obtain 

Sreby, P2v? Be~ "i *T ; (m + B16) + m+ Pale) + m+ 1) 


r= - > 


himc*vo(v? — v,*)? . (27 + 1)(2)7 + 3) 


(m — 1)\(7 —m+ 2)(\j —m+1) 
(27 + 1)(27 + 3) 


with », =v(IIp, '), v2 =v(lTp, *I,) (i-e., (01) in our previous notation), for 


the R branch, and an analogous expression for the P branch. Performing the 
summations, 

Lomeby, P2y2Be~ i+! *Tj(j + 1) 

i = =. : 


Shmec*ve(v? — v,*)* 


j here is the j of the “ly level; the minus sign is to be taken for the R branch, 
the plus for the P. Our theory thus agrees with Wood's observations in giving 
positive (normal) rotations for the lines of the P branch, and negative (ab- 
normal) rotations for the R branch. 


Magnetic rotation spectra of the alkalis 


The calculation of the rotations in the alkali vapours is similar. The lines 
observed by Wood and Loomis*® are 'II—'S transitions. Here the Zeeman 


effect is given by the well-known formula W“ (Ajm) =BA?m_ j(j +1). We have 


** Kronig, reference 3. The phases we use are 


Py(Q, 7, m; Q’, 7’, m +1) = + iP (Q, 7, m; 0, j’, m= 1). 


*® F. W. Loomis, Phys. Rev. 31, 323 (1928); Loomis and Nusbaum, ibid. 38, 1447 (1931); 
39, 89 (1932). 
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v(m’; Ojm) = Bm" hj" j'+1), 


iP? (j-1)(y+m)(j+m—1) 


(j-1,m +1; jm) | 


2 M27 —1)(27 +1) 
: iP? (J+2)( 7+ m+2)\(j+m+4+1) 
[P.jomsj+ ims 1)P,j+1, m4 1; jm) | = - ' 
2 (7+1)(27+1)(274+3) 
iP? (j+m+1)(jtm) 
[Psj,mij.m $1) Py(j,m F1sjm)| = + a 
: J#+V) 


Fe a | P(E, HH) | *. 
The rotations are, in terms of the j of the normal state, 
4rev, P22 Be ™s STF) 


Shme*(v? — v7)? 


fp = GF+2) G+) R branch, 


fj =-U-1) i P branch, 
A) = (23 +1) FG +1) QO branch, 
vy, = v('ll, 'S). 


The signs of the rotations of the R and P branches are just opposite those 
found for iodine. These signs are confirmed by Wood's observations of the 
rotation in Nas. 

It may seem a little surprising that, aside from the Boltzmann factor 


W)*T) the rotations of the P and R branches are nearly independent of j 


e 
when } is large. However it must be remembered that although the over-all 
width of the Zeeman pattern decreases as 1 j, the number of components 
increases linearly with j. The angle of rotation is small, only a few degrees, 
hence the intensities, which are proportional to sin® 6, vary with j approxi- 
mately as }f(j)e") *"{*. As the intensities in the absorption spectrum vary 
as (2j+1)e") *", the intensities of lines of large j in the P and R branches of 
the magnetic rotation spectrum are weighted adversely, as compared to the 
absorption spectrum, by a factor e-"./*", (27+1), and so are, comparatively, 
greatly reduced. It is perhaps worthy of remark that the intensity of the ro- 
tation spectrum is proportional to the square of the density, rather than the 
first power. The intensities of the Q branch decrease much more rapidly than 
those of the R and P branches, which accounts for the fact that the Q branch 
is not observed in the rotation spectrum. For small 7 the R branch is con- 
siderably more intense than the P. This, together with the piling up of the R 
lines at the heads of the bands, explains why sometimes the R branch alone 
is observed. The present theory is thus adequate to account for the principle 
characteristics of the magnetic rotation spectrum. 

I am deeply indebted to Professor J. H. Van Vleck, under whose direction 
this work was done. 


* It should be noted that, since the separations between lines increase with j, the density 
of lines on a photographic plate will decrease even more rapidly than the intensities. 
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Experiments on the Nature of Ferromagnetism 


By FRANctIs BITTER 
Westinghouse Research Laboratories, East Pittsburgh, Pennsylvania 
(Received May 31, 1932) 

Irregularities in the magnetization of ferromagnetic crystals were studied by 
means of magnetic FeO; particles about 1 in diameter suspended in ethyl acetate. 
On iron crystals the particles settle along parailel lines spaced about 0.1 mm apart and 
more or less at right angles to the applied field. On nickel crystals the pattern is 
similar, but the lines split up into more complicated structures as the applied field is 
increased. On cobalt crystals straight lines are obtained on some crystals, and spotty 
patterns on others. The spots arrange themselves in rows when the crystals are mag- 
netized. The patterns can be destroyed by relatively slight surface strains. There is 
some evidence that in the straight line patterns the direction of the lines is related to 
the direction of magnetization. An explanation of these effects is not known. 


HE object of the experiments described below was to find any irregulari- 
ties or inhomogeneities that exist in the magnetization of ferromagnetic 
crystals. 
APPARATUS AND MATERIALS 


The method used wes the old magnetic powder method. In order to 
reveal small details, the particles must be small and have as large a magnetic 
moment as possible without causing them to form chains in the direction of 
the magnetic field. They must further be suspended in some medium that 
will allow them to settle according to the pattern of the field. Particles of red 
FeO; having diameters in the neighborhood of 1u were found very satisfac- 
tory.! These were suspended in ethyl acetate, which, because of its low 
Viscosity, allows the particles to settle fairly rapidly and without undue 
disturbance from currents in the liquid. Observations were made by magne- 
tizing a flat sample in the direction of its length and placing a drop of the 
suspension on the surface of the crystal to be investigated. In order to be 
certain that the patterns obtained were determined by the sample investi- 
gated rather than by the interaction of the FeO; particles, observations 
were made with particles of different size, and of different substances sus- 
pended in various liquids, and with no liquid at all. In all cases the same 
patterns were observed. 

A few samples of iron were prepared by elongating small annealed slabs 
about 3 percent and annealing in hydrogen just below the a—y transforma- 
tion point. This material was very pure, having a total impurity before the 
above treatment of about 0.01 percent, chiefly Si. Some samples were ground 
flat and annealed. All other samples were made by melting in alundum 
troughs in an electrically-heated alundum tube furnace through which puri- 
fied dry hydrogen flowed continuously. In order to observe the effects re- 


'F. Bitter, Phys. Rev. 38, 1903 (1931). 
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ported below, it is desirable to have a shiny, smooth, unstrained surface. 
Hydrogen was used to prevent surface tarnish due to oxidation. Considerable 
hydrogen was absorbed while the samples were molten, and given off during 
soliditication, producing cavities in the ingots and locally imperfect surfaces. 
In most cases, however, at least a part of the sample was satisfactory. 

I should like to express my indebtedness to various members of the 
Laboratories for their assistance—to Mr. L. McCulloch for preparing the 
FeO; powder; to Dr. V. Hicks for making an x-ray examination of this 
powder; to Miss M. Ferguson for taking the best of the photographs repro- 
duced below; to Mr. P. G. MeVetty for stretching the strain-annealed 
samples of iron; and especially to Mr. A. A. Frey, who built the furnaces and 
prepared most of the samples investigated. 


RESULTS 
Iron 
In a sample of iron having small grains, whether annealed or not, ir- 
regular patterns are found. The deposits do not follow any obvious structural 
features, such as grain boundaries. If, however, the grains are somewhat 
larger, and unstrained, deposits like that shown in Fig. 1 are formed, which, 
for larger fields, degenerate into irregular patterns. In Fig. 1 some of the 





Fig. 1. Deposit on a sample of 4+ percent Fe-Si with large grains. Field horizontal. Approximate 
magnification X2. Deposits are black in all photographs. 


grains are crossed by systems of parallel lines more or less at right angles to 
the direction of the applied field. The lines are not spaced perfectly evenly, 
but are of the order of 0.1 mm apart. Fig. 2 shows the deposit obtained near 
a grain boundary on an Fe-Si sample, and Fig. 3 shows the slip lines obtained 
on a similar sample after a slight reduction by cold rolling. The similarity is 
rather striking. In both cases there is a tendency for the lines to continue in a 


changed direction after crossing a grain boundary. 

In general, the deposits on iron are indefinite, and it often occurs that 
no regular patterns are observed at all. Better results were obtained with an 
alloy containing 4+ percent Si than with pure iron. An attempt was made to 
determine how the pattern changed as the direction of magnetization was 
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varied. As a rule the pattern was obscured, but on several crystals three sets 
of lines were observed all told. Further observation is required to establish 
the generality of this result. 





Fig. 2. Fig. 3. 
Fig. 2. Deposit on a 4 percent Fe-Si sample near a grain boundary. Field horizontal. Approxi- 
mate magnification X10. 
Fig. 3. Slip lines in a + percent Fe-Si sample. Approximate magnification X10. 


Cobalt 

Typical deposits on cobalt are shown in Figs. 4-6. These illustrate the 
state of affairs before a magnetic field is applied. In Fig. 4 the spots vary a 
good deal in size. Fig. 5 shows grain boundaries, and the tendency for the 
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Fig. 4. Fig. 5. 


spots to form a hexagonal array. Also a grain having a straight line pattern 
is shown. Fig. 6 shows transition patterns between the spotty and the straight 
line types. Just what difference exists between samples showing different 
patterns is not known. One sample containing 20 percent of iron showed 
no pattern whatsoever. Until further evidence is obtained, we may make the 
reasonable assumption that the irregularities are produced by strains, im- 
purities, or other imperfections. In Figs. 7-9 the effect of magnetization on 
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Fig. 7. Deposit on magnetized cobalt near a grain boundary. Approximate magnification X80. 
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the spotty type of deposit is shown, the magnetization being stronger in 
Fig. 9 than in Fig. 8. Fig. 7 shows the distorting effect of the grain boundary. 
The line pattern illustrated in Fig. 5 is affected little if at all by magnetiza- 
tion, and has in no case been observed to change its direction on the crystal 
surface as the direction of the applied field was varied. 





Fig. 8 and 9, Deposits on a magnetized cobalt crystal. Approximate magnification 160. 


Nickel 


A typical deposit on nickel before being exposed to any magnetic field 
other than that of the earth is shown in Fig. 10. On magnetization the pat- 
terns become more definite, as shown in Figs. lla and 11b. Sometimes the 
lines are curved, especially near grain boundaries, and sometimes they show 
a break as is the case for one of the lines in Fig. 10. It quite often happens 
that the spacing of the lines changes progressively from one side of a crystal 
to the other, as can be seen in Fig. 11a. These lines are not associated with 
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fixed places in a crystal. This is illustrated in Fig 
three different deposits on the same region are 


.12 in which photographs of 
superposed. It seems quite 


generally true that the position of the lines changes after demagnetization. 
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Fig. 10. Deposit on a nickel crystal before magnetization. 
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Approximate magnification X80. 


Some samples in the shape of disks were prepared, and magnetized in eight 


directions at 45° intervals. The deposits were photographed and the total 


number of directions in which lines appeared was determined for tive crystals. 
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Figs. 1la and 11b. Typical deposits on nickel with the applied field acting in a horizontal 


direction. Approximate magnification X65, 


In each case the lines appeared in four directions all told. In those positions 


in which two such directions made nearly equal angles with the applied 


field chevron patterns like that shown in Fig. 13 


appeared. 
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If successive patterns are obtained in the same region in gradually in- 
creasing fields without demagnetization between observations, the phenom- 
enon shown in Fig. 14 is found. The wavy marks on the photographs are 
surface irregularities that are of no interest here except insofar as they may 
be used as a reference system fixed in the crystal. In the twelve photographs 
shown the applied field increases progressively, being smallest in 1 and largest 
in 12. The position of the lines as shown in the first photograph does not 
change radically. The lines, however, split up as the magnetizing field is in- 
creased, first into two, then four, and on the original photograph, No. 12, 
each of the original lines is seen to consist of eight components. In addition a 
new line appears between each of the main groups in photograph No. 7 and 
splits into two components also. The smallest splitting here shown is deter- 
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Fig. 12. Superposed photographs of deposits Fig. 13. Peculiar pattern obtained on 
obtained on the same region at different times, a nickel crystal. 


the sample being demagnetized after each 


observation. 


mined by the method, and it seems very probable that much more detail 
exists than is here indicated. To obtain these results it is essential not to 
demagnetize the sample between observations. The sample consisted of 
crystals having linear dimensions of from 1 to + cm imbedded in a large slab, 
27 cmX+ cm X1 cm. It was consequently possible to measure only the 
intensity of magnetization of the aggregate. This was done with a fluxmeter 
and the value of /=350 or about 75 percent of saturation was obtained for 
the conditions under which deposit No. 12 was produced. Magnetization of 
this intensity is confined to the directions of easiest magnetization within the 
crystal. 

Alloys containing 10 percent, 20 percent, and 30 percent of iron were pre- 
pared by melting and slow cooling. That containing 10 percent of iron 
showed patterns typical of nickel, but less distinctly than usual. No patterns 
whatsoever could be found on the other samples. 

The effects are exceedingly sensitive to strain. Some samples that were 
lightly pressed by a C-clamp no longer gave any patterns, nor was anything 


to be found on polished samples until they had been annealed. 
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DiscUssiONn 
If it is quite generally true, as has been found in a few cases, that the line 
patterns can appear in three directions in iron, in four directions in nickel, 
and in only one direction in cobalt, it seems reasonable to assume that they 
are somehow related to the three (100) axes in iron, the four (111) axes in 
nickel, and the single (0001) axis in cobalt. This assumption is supported by 
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Fig. 14. Twelve photographs of the same region, the applied field being 
gradually increased between observations. 


the fact that in each case the axis mentioned has a special significance mag- 
netically—it is the direction of easiest magnetization within the crystal. 
Taking into consideration the above material, and the fact that the lines are 
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sometimes curved or broken, we suggest the hypothesis that inhomogeneities 
occur along planes perpendicular to the direction of magnetization, and that 
the lines observed are parallel to the intersection of such planes with the 
surface under observation. In the fields used the direction of magnetization 
probably differed in no case appreciably from the direction of easiest mag- 
netization. The above hypothesis fits into the observation that in nickel and 
cobalt alloying tends to destroy the pattern, for the effect of a foreign sub- 
stance must, as a rule, be to introduce local fluctuations in the direction and 
intensity of magnetization. This reasoning can be applied to the observations 
on iron if we suppose that it is exceptionally difficult to produce a satisfac- 
torily uniform crystal of iron, and that the addition of silicon tends to re- 
move certain imperfections. The above discussion does not take into account 
the spotty pattern obtained on certain cobalt crystals. It seems likely that 
this is more complicated in its origin, and occurs whenever the hexagonal 
axis is normal to the surface of the crystal. 

From the observations so far made, it is not possible to determine the na- 
ture of the inhomogeneities revealed. They may be changes in direction or 
intensity of magnetization, or both. Bearing in mind that in this discussion 
a change in direction, to be detected as such, must cover a region large com- 
pared to 1p, it is probable that we have to do with changes in direction in 
cases like that shown in Fig. 13, but that in patterns like that shown in 
Fig. 11 there is only a change in intensity. It is most important to settle 
definitely the nature of the inhomogeneities by further and more detailed 
observations. 

The origin of the patterns obtained is of course even more obscure. The 
conception is current that ferromagnetic crystals are made up of spontane- 
ously magnetized regions whose directions of magnetization are distributed 
among the various energetically possible orientations. Magnetization is then 
supposed to result partly from a growth of some of these regions at the ex- 
pense of their neighbors,? partly from changes in the direction of magnetiza- 
tion of each region. It seems a little difficult to fit this picture into the above 
photographs unless the regions are smaller than 1u. The secondary structures 
suggested by Zwicky® offer a different attack on the problem. Finally, we 
must remember that a crystal consists of atoms in different excited states, 
that the various periodic motions which its constituent atoms and electrons 
perform most probably are coupled together in a variety of ways, and that 
in addition there must be a certain amount of electromagnetic radiation 
present. Any of these factors may have an influence in determining the equi- 
librium configuration of a crystal, and it is perhaps the most important aspect 
of experiments such as the above to furnish evidence from which conclusions 
on these points may eventually be drawn. 


2 F. Bloch, Zeits. f. Physik. 74, 295 (1932). 
°F. Zwicky, Phys. Rev. 38, 1772 (1932). 
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Thermoelectric Power of Single Crystal Bismuth 
Near the Melting Point 


By AvoLF Soroos 
Physical Laboratory, University of Iowa 
(Received June 14, 1932) 

By investigating directly the thermoelectric power of bismuth single crystals near 
the melting point a transition region is found in which the thermoelectric power 
changes gradually from that characteristic of the solid to that characteristic of the 
final liquid state. The region extends from the melting point to about 8° above it. 


As an explanation it is supposed that a definite crystal structure persists into the liquid 
state. 


INTRODUCTION 


N measuring the thermal e.m.f. of bismuth single crystals against constan- 

tan as a reference metal, Boydston! found certain irregularities in the 
e.m.f.-temperature curves in the region just beyond the melting point. As an 
explanation he says, “It appears, in fact, as if some crystalline arrangement 
still persists in the molten metal as shown by a transition region between the 
thermal e.m.f. of solid and liquid.” 

Goetz? makes use of a similar explanation for certain effects observed in 
his study of the factors influencing the growth and perfection of bismuth 
crystals. He found that if he used a seed crystal with regions of twinning to 
inoculate a molten rod he would invariably obtain a crystal which still showed 
the twin orientation unless he melted the seed back considerably beyond the 
region of twinning so that this region was at a much higher temperature than 
the melting point. He also found that a polycrystalline rod melted and raised 
to a temperature just slightly above the melting point, recrystallized, on 
cooling, into a rod which had the “same position and same orientation of the 
prominent crystal elements” as the original rod. If he heated the specimen 
considerably above the melting point, however, he found an entirely new 
arrangement on recrystallization. 

Although Boydston observed only slight deviations of thermal e.m.f. from 
what is expected on the basis of an abrupt change at the melting point, it is 
obvious that these deviations must be accompanied by relatively large 
changes in the derived thermoelectric power curve. For this reason it seemed 
advisable to investigate directly the thermoelectric power of bismuth crystals 
over a temperature range extending on both sides of the melting point with 
the hope of definitely establishing the phenomenon foreshadowed by Boyd- 
ston. 

The method devised to do this consists in measuring the thermal e.m_f. 
over a small temperature interval at successive mean temperatures in the 


1 R. W. Boydston, Phys. Rev. 30, 911 (1927). 
* A. Goetz, Phys. Rev. 35, 193 (1930). 
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desired region. The observed e.m.f. divided by the temperature interval gives 
directly the mean thermoelectric power at the mean temperature. By use of 
a small enough temperature interval the resolution can be made sufficiently 
high to justify considering this mean thermoelectric power as the actual. 


APPARATUS AND METHOD 


The experimental arrangement is shown in Fig. 1. The bismuth crystal, 
BB, has its ends fused to two copper blocks, AA. In fusing, the copper surface 
is cleaned with a flux (ZnCl) and then with the end of the crystal pressed 
lightly against the block, the latter is heated until the bismuth just begins to 
melt. A little added pressure with the finger tip at this point insures a good 
contact. The crystal structure is not altered at the junction. Two holes a, } 
were drilled into the ends of each copper block, with the bottoms of the holes 
separated from the surface onto which the bismuth was fused by a thickness 
of copper of about 0.5 mm. Into the holes, aa, were inserted the two junctions 
of a constantan—copper thermocouple which measured the difference in 
temperature of the two ends of the bismuth rod. These junctions were en- 
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Fig. 1. Diagram of apparatus. 


closed in small Pyrex glass tubes, the copper and constantan wires being 
separated by thin strips of mica. To the bottoms of the holes, bb, were silver 
soldered copper wires which completed the bismuth-copper thermal e.m.f. 
circuit. These wires were insulated from the side walls of the holes by small 
Pyrex tubes. This latter procedure was necessary in order to make sure that 
the effective junctions of this bismuth-copper couple should be at the surface 
of separation between the copper block and bismuth crystal. A second copper- 
constantan circuit had one junction placed at C, the other junction being 
maintained at zero degrees. 

The electrical circuits are also pictured in Fig. 1. The reading of the 
galvanometer, G;, is proportional to the temperature difference between the 
ends of the bismuth crystal, the reading of galvanometer, Gp, is proportional 
to the e.m.f. of the bismuth-copper circuit, and the reading of the galva- 
nometer, G3, is proportional to the mean temperature of the crystal. 

The crystals were grown by a method similar to that described by 
Kapitza.’ A groove‘ in an asbestos block was used for a mold. Crystals were 


3 P. Kapitza, Proc. Roy. Soc. A119, 358 (1929). 
‘ The cross section was almost V-shaped, but with a rounded bottom. 
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grown of three brands of bismuth, Mallinckrodt C.P.,5 a special silver-free® 
and Kahlbaum’s best. 

The results given below are for only the first brand, preliminary work? 
having shown no essential difference in the behavior of the others. The 
crystals were about 0.5 cm? in cross section and either 10 or 20 cm long. The 
longer ones were cut in half and a separate determination made on each piece. 

In making a determination the crystal was put into an asbestos holder 
and placed in a horizontal furnace. The temperature gradient in the furnace 
was so controlled by auxiliary coils that the temperature difference between 
the ends of the crystal could be kept at about 1°C. The furnace temperature 
was allowed to pass slowly through the region under investigation. The time 
required to take the crystal through the desired temperature range (260° 
to 290°C) was about one hour. 

RESULTS AND DISCUSSIONS 

Measurements were made on six single crystals ranging in orientation’ 
from 29° to 69°. These results are shown in Fig. 2. All crystals at a tempera- 
ture sufficiently above the melting point, (i.e., in the region 280°-290°) 
showed a constant thermoelectric power against copper of 4.5 microvolts per 
degree (within 1 microvolt per degree). This has therefore been taken as a 
zero point and the ordinates of the curves depict the excess values above this 
constant level. While the curves are plotted against mean temperature, this 
is not obtained directly from the readings of the thermocouple, C, Fig. 1. 
This couple could not be placed in thermal contact with the crystal, and while 
for any one crystal it is believed that its readings give a correct measure of 
the rising temperature these readings may differ from the actual temperature 
of the center of the crystal by a presumably constant and small but undeter- 
mined amount. The observed curves have therefore been so adjusted in 
plotting Fig. 2 that they agree in mean temperature at the point where the 
final constant value is reached. This adjustment does not amount for any 
crystal to more than about 2° difference between the adjusted temperature 
scale and the indication of thermocouple, C, Fig. 1. The location of the melt- 
ing point (271°C) as shown on the figure actually coincides with the reading 
of thermocouple C in the case of the 40° crystal. Moreover, it occurs, for this 
crystal, about 0.5° above the first sharp break which is then interpreted as 
coinciding with the beginning of the melting process. For no other crystal is 
the onset of melting so sharply marked. Even without such an adjustment, 


5 This is the same grade of bismuth that Boydston used and has the same stated analysis. 
Other work in this laboratory has shown that the silver content may be much less than the 
0.04 percent maximum limit. 

® The writer has to thank the U.S.S. Lead Refinery for this material. The analysis given 
with it shows: Pb 0.005 percent, Ag 0.0006 percent, Cu, Fe 0.000 percent. 

7In this preliminary work about fifteen crystals were investigated, the experimental 
set-up being improved from time to time. As soon as a reliable arrangement was obtained the 
set of crystals reported on herein was grown and measured. The preliminary results were all 
in agreement with the final ones in their general characteristics. 

’ Angle between the trigonal axis and the length of the specimen. 
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however, the general character of the results is unchanged. To avoid confu- 
sion, the curves in Fig. 2 are plotted with an arbitrary displacement in the 
scale of ordinates between successive curves. Four of the crystals were suffi- 
ciently long for two determinations to be made, the circles indicating deter- 
minations made on the half onto which the nucleus had been fused during 
growth and crosses indicating the second half of the same crystal. 
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Fig. 2. Thermoelectric power as a function of temperature. 


It is clear from Fig. 2 that bismuth is abnormal in its thermoelectric be- 
havior in the region just above the melting point. For most substances which 
melt at a definite temperature, the thermoelectric power changes abruptly 
at the melting point or the thermal e.m.f. curve suffers an abrupt change in 
slope.*'’ For the crystals shown here, however, the curves before settling 
down to a practically constant final liquid value pass through a transition 
region which lasts for approximately 8°C. This is exactly the effect predicted 
by Boydston. This transition of thermoelectric power, of course, takes place 
entirely in the liquid state since the crystal is completely melted when the 
mean temperature (temperature at center of crystal) is 0.5° above 270°, since 

* E. G. Linder, Phys. Rev. 29, 554 (1927). 

10 T, R. Harrison and P. D. Foote, J.0.S.A. and R.S.I., 7, 389 (1923). 
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the difference in temperature between the two ends of the crystal is only one 
degree. 

On account of the doubt of the exact indication of thermocouple C and 
from the general appearance of the curves in Fig. 2 one might be inclined to 
believe that the melting process was not complete until the final constant 
value of the thermoelectric power had been reached and that the observed 
effect all lies in a region below the melting point and is, therefore, charac- 
teristic of the solid and not the liquid state. That this is not the case was veri- 
fied by an independent experiment. A separate determination was made on a 
90° crystal and as all of the curves show a small constant portion shortly after 
melting, the determination was stopped when that place was reached. The 
holder containing the crystal was then quickly taken from the furnace, 
tipped on its side and jarred a little. There was only one possible place where 
the liquid could seep out through the holder and this was at the position where 
the thermocouple C was inserted. It was found that the liquid had indeed 
forced its way out and on examining the material remaining in the holder, it 
was found recrystallized into a changed orientation" from what it was at the 
beginning. This seems conclusive evidence that the transition region belongs 
to the liquid state. 

It seemed desirable also to make a check on the apparatus to be sure that 
the results obtained were not due to some peculiarity in the apparatus itself. 
To check this point a determination was made on a specimen of polycrystal- 
line lead. The results of this determination are shown in the curve in the 
upper right corner of Fig. 2. Since the thermoelectric power of lead against 
copper is small, a temperature difference between the ends of the specimen of 
4°C had to be used in order to obtain a conveniently large galvanometer de- 
flection. This made the resolution much less than that for any of the bismuth 
crystals. The result obtained for lead, however, is the one to be expected on 
the assumption of no transition region existing. The thermoelectric power 
passes directly from the solid value to the liquid value without a change of 
slope in this region and the temperature interval of this change corresponds 
almost exactly to the temperature difference of the ends of the lead rod. The 
fact that the transition is not shown by a line parallel to the ordinate axis, is, 
of course, due to the finite resolution. 

The above results yield incidently a value of 2.6 microvolts per 1°C for 
the thermoelectric power of solid against liquid lead at the melting point. 

As a matter of interest a plot has been made, Fig. 3 of the thermoelectric 
power of the solid crystals just below the melting point against the square of 
the cosine of the orientation angle. To do this ordinates were read from Fig. 2 
for the value of temperature shown by the dotted line. This is 0.25° below the 
temperature at which the hotter end of the crystals start to melt. The points 
lie on a straight line, confirming the Voigt-Thomson relation. This is in 
agreement with Boydston though at variance with the slight deviation from 


" Doubtless from tipping and jarring. From Goetz’s observation it should remain the 
same if undisturbed. 
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the law later found by Bridgman” and Fagan and Collins. It does, however, 
give the writer considerable confidence in the method used to adjust the 
mean temperature scale. A numerical comparison“ of Boydston’s and the 
writer's data shows that at the melting point the former obtains: 


ey — e; = 70 microvolts/degree!® 
¢: — ¢, = 29 microvolts/degree 


while the writer obtains: 




















ey — €; = 65.5 microvolts/degree 

e: — €; = 31.0 microvolts/degree. 
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Fig. 3. Thermoelectric power as a function of cos? 6. (a) At the melting point. (b) At the knee 
of the curves of Fig. 2. 


The results then show that the thermal e.m.f. curve of bismuth does not 
change its shape as abruptly and in as clearly defined a fashion as in the case 
of zinc, tin, lead and some other substances. The thermoelectric power be- 
fore settling down to a final constant value passes through a transition region 
lasting for a temperature interval of about 8°C above the melting point. 
Moreover, this settling down is gradual. The curves for differently oriented 
crystals preserve their individuality for several degrees above the melting 
point. Although it is only quite close to the final liquid stage that the curves 
become practically identical there is within the transition region a feature 
common to all but one of the curves of Fig. 2, namely a slight knee or prac- 
tically stationary value of thermoelectric power. This seems to indicate a 
definite phase, or intermediate liquid state, existing for a temperature range 
of about one degree. If this is the case it ought to occur at the same tempera- 


® P. W. Bridgman, Proc. Amer. Acad. of A. and S. 63, 351 (1929). 

1H. D. Fagan and T. R. D. Collins, Phys. Rev. 35, 421 (1930). 

4§ Boydston’s values for comparison were obtained by plotting his data for e, —e, against 
cos? @ and then extrapolating the straight line to the ordinates of 0° and 90° crystals. 

1 e=dE/dT, the thermoelectric power, E being the thermal e.m.f., e; =thermoelectric 
power of liquid bismuth, e,,=thermoelectric power of crystal of 0° orientation, e, =thermo- 
electric power of 90° crystal. 
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ture for all the crystals. Although it obviously does not do so its significance 
can be tested in another fashion. This is shown in Fig. 3b in which the value 
of the thermoelectric power at the knee is plotted against cos’ 6. A definite 
intermediate liquid phase might be expected to give either (1) a straight line 
(dotted in the figure) with a positive slope indicating that the liquid is still 
very much like its parent crystal (2) a straight line (solid line in the figure) of 
zero slope, indicating that an identical state has been reached for all crystals. 
Unfortunately the observed data do not decide unambiguously between these 
two suppositions. In spite of this the existence of a knee at all strongly sup- 
ports the idea of an intermediate state. It may be noted that this knee was 
also obtained for most of the crystals which were studied in the preliminary 
work previously mentioned. 

The writer believes that the existence of a transition region extending 
from the melting point to about 8° above it, is established by the foregoing and 
that the explanation which has already been advanced is adequate to explain 
this phenomenon. Up to 8° beyond the melting point the liquid retains, but 
in decreasing measure, the organization of the crystal from which it is de- 
rived. Beyond that point the state of the liquid is independent of its previous 
history. It is practically impossible to think of any state or condition of the 
liquid, other than that postulated, which could give rise to such a transition 
region. It must be noted however, that the thermoelectric effect must be 
particularly well suited to observe this phenomenon as the appearance of a 
transition region for other properties which depend on crystal orientation has 
not been observed, so far as the writer knows. W.L. Webster's study” of the 
change in magnetic susceptibility in solidifying bismuth apparently shows no 
transition region for this quantity of the type shown here, but a sudden 
change occurring at the melting point. It may be, however, that a degree of 
organization suitable for modifying the thermoelectric effect may have little 
effect on susceptibility. Moreover, \Webster was looking for an effect said to 
occur 10° below the melting point and not for one above. 

It may be said finally that only x-ray diffraction can give positive and 
unambiguous proof of the explanation offered to explain the transition region. 

The writer wishes to express his most sincere thanks to Professor E. P. T. 
Tyndall, who suggested the problem and under whose direction the work was 
carried out. 


16 W. L. Webster, Proc. Roy. Soc. A133, 162 (1931). 
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The Acoustic Resonator Interferometer: II. Ultrasonic Velocity 
and Absorption in Gases 


By J. C. HUBBARD 
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The derivation of the equivalent electric network of the acoustic resonator inter- 
ferometer in Part I of this paper has made it possible to develop the theory for the 
current in a simple resonant circuit in which the electrodes of the piezoelectric plate 
of the interferometer are connected to the terminals of the variable condenser of the 
circuit. The special case of this theory in which the circuit is excited at a constant 
frequency determined by the crevasse frequency of the resonator plate in its given 
situation with respect to electrodes and associated circuit, when the acoustic path in 
the interferometer is detuned and the resonant circuit is tuned so that its resonant 
maximum occurs at the same frequency, takes an especially simple form and leads to 
a direct procedure for determination of ultrasonic velocity and absorption in a gas in 
terms only of current in the resonant circuit and path-length in the interferometer, 
all circuit and interferometer constants dropping out. The values of current as a func- 
tion of path-length obtained experimentally are in complete accord with the theory 
and data for ultrasonic absorption in air and in CO: so far obtained are in agreement 
with the meager data available by other methods. The role of the coefficient of re- 
flection at the fluid-reflector surface is discussed. 


I. EXPERIMENTAL OBJECTIVES AND METHODS 


HE vibration of a column of fluid of variable length, one end of which is 

bounded by a vibrating plane face of an electrical vibrator such as a 
piezoelectric plate, the other end being bounded by a plane reflector parallel 
to the vibrating face of the plate, has been considered in Part I' of the 
present paper, and the equivalent electric network of the system was there 
shown to be the same as that of the quartz plate alone with modified re- 
sistance and capacity coefficients. Such a system, associated with an ap- 
propriate driving circuit, may be termed an acoustic resonator interferometer, 
and has already been used for the measurement of compressional velocities 
in liquids,?**> and in gases,® of absorption coefficients in gases,® and of the 
modification of velocity in liquids due to scattering.’ Specific consideration 
of some of the methods of measurement of absorption in gases and of the role 
of the reflection coefficient at the boundary of the fluid will be given in the 


1 J. C. Hubbard, Phys. Rev. 38, 1011-1019 (1931). Part I will be referred to hereinafter 
as I. 

2 J. C. Hubbard and A. L. Loomis, Phil. Mag. 5, 1178-1190 (1928); and A. L. Loomis and 
J. C. Hubbard, J.0.S.A. and R.S.I. 17, 295-307 (1928). 

3 E. B. Freyer with J. C. Hubbard and D. H. Andrews, Jour. Am. Chem. Soc. 51, 759-770 
(1929). E. B. Freyer, ibid. 53, 1313-1319 (1931). 

* E. Klein and W. D. Hershberger, Phys. Rev. 37, 760-774 (1931). 

5 C. R. Randall, Bur. Stand. Jour. Res. 8, 79-99 (1932). 

6 J. C. Hubbard, Phys. Rev. 35, 1442 (1930); 36, 1668-1669 (1930). 
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present paper. Mr. W. D. Hershberger’? has presented independently an 
analysis of the motion of the fluid column and of its reaction on the piezo- 
electric plate, his treatment being somewhat more restricted than that of the 
author, particularly in that he does not consider the effects of coefficients of 
reflection at the fluid-reflector surface of value less than unity. 

Among the methods of measurement which have been used, undoubtedly 
the simplest is that in which the interferometer electrodes are connected to 
the terminals of a variable condenser which is part of a simple resonant cir- 
cuit driven by an oscillation generator of suitable frequency. The measure- 
ments consist in this case of a determination of current in the resonant circuit 
as a function of the length 7 of the fluid column.?*4-° The current readings 
pass through sharp maxima (for gases, minima for liquids) or peaks as r is 
increased, these peaks having been interpreted as marking the value of 7, for 
resonance in the fluid, successive resonant positions of the reflector being one- 
half wave-length apart. An analysis of this method for gases is given in the 
present paper. 

In a second method, far more sensitive than the first, the interferometer 
electrodes are connected to a secondary coil coupled to the primary of an 
oscillation generator. As 7 is varied a cyclical variation of frequency of the 
generator takes place which is compensated by the variation of a vernier con- 
denser in parallel with the capacity of the generator, constancy of frequency 
being assured by using the double heterodyne beat method.*?** From the 
variation of vernier condenser readings as a function of 7 the acoustic be- 
havior of the medium is deduced. The theory of this method will be given in 
a later paper. 

A third method® makes use of the current in the leads to the interfer- 
ometer itself as a function of r, the circuits being otherwise as in the first 
method described above. Klein and Hershberger* report a variety of experi- 
ments in which a number of methods were emploved, the interferometer 
system being treated as an impedance in alternating current measurements. 

It is essential that the power dissipated in the acoustic system be so small 
as to produce negligible temperature effects. In practice it has been found 
possible to reach the limits of accuracy imposed by the sensitiveness of the 
thermocouple and galvanometer used in the first method with amounts of 
power several orders smaller than that which is capable of influencing the 
results. The power necessary in the second method is still much smaller. 

The most important factor causing a departure of the acoustic system 
from that theoretically expected is to be found in the conditions at the elec- 
trodes of the piezoelectric plate. This factor need not be considered if the 
electrodes are of sufficiently small mass and if their attachment to the plate 
is of such nature as not materially to modify the crevasse. As will be seen the 
constants of the crystal including the effect of electrode mounting upon the 
position and shape of the crevasse may be eliminated from consideration in 
determinations of velocity and absorption. 


7W. D. Hershberger, J. Acous. Soc. Am. 3, 263-268 (1931). 
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ACOUSTIC RESONATOR INTERFEROMETER 


Il. Meruop oF FoRCED VIBRATIONS. ULTRASONIC VELOCITY AND 
ABSORPTION IN A GAS 
As shown in I the acoustic interferometer may be represented by an 
electric network consisting of a fixed capacity K, shunted by a series consist- 
ing of the fictitious inductance L, capacity K’ and resistance R’, where 


R’ = R+ ABprP (1) 


1/K’ = 1/K + ABprwO (2) 


Rand K being the fictitious series resistance and capacity respectively of the 
quartz plate itself, A the cross section of the fluid column, B a constant of 
quartz, and p and 7, respectively the density of the fluid and the wave velocity 
in it. 

P and Q have the values, in general, with x=0 in I, Eqs. (12) and (13), 
respectively 


i me y*e— tra oe y(1 — ye ara COs (2rw v) 


P= a (3) 
1 — 2y*e~?"* cos (2rw/v) + y*te"*"" 
y¥(1 + y) sin (2rw 7) 
QO = (4) 
1 — 2y2e7 7" cos (2rw 7) + yte7*"" 
—_ ‘cena maedasts =a os 


Fig. 1. Resonator circuit with interferometer replaced by its equivalent electric network. 


where a is the attenuation factor of particle velocity in the fluid and is equal 
to one-half the coefficient of absorption, and y is the coefhcient of reflection 
of particle velocity between the fluid and the reflecting plate. The distance r 
is that between the vibrating plate and the reflector, i.e., the length of the 
fluid column. 

Let the electrodes of the piezoelectric plate be connected to the terminals 
of the variable condenser C of the resonant circuit L;, Ri, C, Fig. 1. The re- 
sistance R, includes the resistance of the vacuum thermocouple 7 which is 
inserted between C and L; on the grounded side of the interferometer. Let 
the inductance 1, be loosely coupled to the output of an oscillation generator 
which may be kept at any desired frequency, so that in L; we have the in- 
duced e.m.f., E = Ey e’*'. We have for the impedances of the various branches 
of the circuit, putting @=jw: Z=1,0+Ri; Z:=1/K8; Z72=L04+R'+1/K'8; 
Z3=1/C8. 

Solving for io, the amplitude of current in the 1,R,; branch, including the 
thermocouple heater, we have 
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. JS 2 R29 = rg kK’ j \ 2 
ko RK’? + 41 LK‘w*? + w 
) cak 


f 
_ } (3) 

(Pl — (C+ Ki Liew? 
) C+K (1 — LK’w*) — (weer i> meee oi)" 

1 

a — (C+ A))Lyw* K y 

+ R’K’ +4 RA a LK'w?)  - sania 

C+ A, C+ with Se 


Putting J)>= Fy Ry, the maximum value of 7) at resonance with the inter- 
ferometer disconnected, and also putting 


p=1-(C+ Aj)L oe? (6) 

q = 1 — LK'w*, go = 1 — LK? (7) 

@; = RC + Aj)w (8) 

b>: = R'K'w (9) 
and 

g = io/I 


we have 


o:2|o22 + (¢ + K’(C + K)))?| 
o- - (10) 


~ [pg — dda — (1— pK (CH+ KDE + lop togt kK (C + KDE 


Eq. (10) is an exact description of the current amplitude in the thermo- 
couple heater relative to the current amplitude at resonance in the resonant 
circuit with the interferometer disconnected, as a function of the decrements 
of the oscillating circuit and of the interferometer system, of the departure 
from resonance in the two systems, of the pure capacity C+ K,, and of the 
fictitious capacity K’ of the acoustic system. This equation is similar to the 
one deduced by D. W. Dye® for the piezoelectric resonator in a resonant 
circuit but is complicated by the inclusion of the modified R’, K’ and 2 
terms introduced because of the interferometer system. Dye’'s analysis in- 
cludes a capacity Ke between the resonator and the resonant electrical circuit 
to account for the effect of variable electrode distance from the quartz plate. 
This capacity is here suppressed, the electrode conditions being maintained 
constant. 

In view of the exhaustive study which Dye has made of his equation no 
extended discussion of its analogue, Eq. (10), is necessary here. It should be 
pointed out, however, that in general o? as a function of the frequency rises 
to the usual resonance maximum for p=0, and that further, if g+K’/(C+K,) 
= (0 in the same neighborhood there will be a sharp minimum or crevasse in the 
resonance curve. This crevasse is extraordinarily narrow because of the ex- 
tremely small value of ¢@2 as compared to @¢;, and at its center o° falls to a 
small fraction of its value on either side. By repeated adjustment of C and w 
it is possible to find the value of the frequency for which pand q+ K’'/(C+ K,) 


§ D. W. Dye, Proc. Phys. Soc. Lond. 38, 399-458 (1926). 














ACOUSTIC RESONATOR INTERFEROMETER 


1 
bo 
~ 


. 


are zero simultaneously, provided the reflecting plate is set at a distance 
r=(2n+1)d/4, so that P and Q are both equal to zero and K’=K, making 
qg=qo. In practice it is sufficient to find a crevasse for any setting of r, where- 
upon a rough determination of \/2 is made. The setting of r is now made at 
\/4, then by successive adjustment of C and w the crevasse is made sym- 
metrical. The interferometer is then disconnected and C is varied, w being 
kept constant, until 7) has the maximum value Jo, which is determined. The 
interferometer is then reconnected and C is adjusted to its former value for 
a symmetrical crevasse, which completes the adjustment necessary for a 
determination of ¢? as a function of r. A series of observations of 7 taken under 
these conditions is shown in Fig. 2, curves I and II, in which only maxima 
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Fig. 2. Curves I and II, maximum and minimum galvanometer readings as a function of r. 


Curve III, representation of Eq. (17) as a function of the same observations. 


and minima of current are shown. The behavior in the immediate neighbor- 
hood of a maximum for the same case is shown in Fig. 4, Insert a, the abscissa 
scale being magnified twentyfold. In the preliminary adjustments just des- 
cribed, we have made p=0, g+K/(C+K,) =0, and, remembering Eqs. (1), 
(6) and (7),q¢+K’'/(C+K,) =K’A BpvwQ, so that 


. (1+ SP)? +5 
[1+ SP + 1/Rwoi(C + Ki) ]? + S90 











o 


(11) 
where 
S = ABpo/R. (12) 


When the reflector is at a distance r=(2n+1)\/4 from the face of the 
vibrating plate, if 7 is small, P=zero to a very close approximation, and 
O=0.° Let oo be the value of o for this reflector position. We have then 


* See I, Figs. 2 and 3. These figures were drawn for the case in which 7 was assumed to be 
unity. As we shall see later the statement made above is practically always valid for gases. 
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o, = 1, [1+ 1/Rwd(C + K,)| 


or 
Red (C + Ay) = oo (1 — ao). (13) 
The value of ¢ given by Eq. (11) under the prescribed conditions will have 
maximum values, ¢,,, for values of r such that Q=0 and P= P,, given by 


Pm = (1+ ye) (1 = yer), (14) 


! 


We then have 


II 


o» (1+ SP,,) [1+ SP, + (1 — oo) oo], or 


S = ABpy R = [o,,(1 — oo) (1 — ooo — 1] P» = constant. — (15) 


Putting y=e 8, Eq. (14) becomes P,, =e8? [cosh (ra+8/2)]| sinh (ra+s), 
or, toa very close approximation, unless r is large 


Pm = (1 + B/2) (ra + B). (16) 


This value of P,, inserted in Eq. (15) would lead at once to a determination 
of a and 8 if we had the means of accurate determination of A, assumed to be 
the area of the vibrating plate exposed to the fluid. As is well known, only a 
small portion of the surface of a piezoelectric plate as ordinarily cut is in 
Vibration in a given mode!®:!! though by optical means it can be shown that the 
entire cross section of the column of a fluid in an interferometer is set into 
resonant vibration and presumably reacts upon the crystal over the whole of 
its exposed face. Hershberger™ has recently pointed out that the assumption 
that A is to be taken as the area of the plate exposed to the fluid is not valid. 
As we shall see, however, it is sufficient to assume an effective 4A, which is 
constant. In the absence of knowledge of its magnitude we have then in 
Eq. (15) only a means of computing the relative effects of absorption and 
reflection. For, substituting P,, from Eq. (16) into (15), we have 


[1 — o0),/c0] onl — on) — o0/(1 — 00) |(ra 8 + 1) = 8’ (17) 
where 
S’ = S(1 + 8/2) 8 = constant. (18) 


Eq. (17) enables us at once to compute a/8 from a series of observations of 
values of ¢,, and the corresponding values of r, and the mean value of a. 
Fig. 2, curve III, shows the observations of Fig. 2, curves I and IT, in terms 
of Eq. (17), the smooth curve passing through the maxima being the hyper- 
bola given by Eq. (17), where a/8=6.276, and S’=2.548. All types of 
acoustic interferometers which have been used by the method described here 
yield results for gases in entire accord with the theory showing that within 
the limits of errors of observation the quantities entering into S’, Eq. (18), 
may be considered as constant. 
10 Nat. Phys. Lab. Report for 1928. 


1 P, Vigoureux, Quartz Resonators and Oscillators, London (1931). 
" W. D. Hershberger, Physics 2, 269-273 (1932). 














ACOUSTIC RESONATOR INTERFEROMETER 529 


Having determined a/f from a given series of observations, it remains 
to be seen how a and 6 may separately be determined. Putting (1 —¢9)/ao =a, 
combining Eqs. (11) and (13), and solving for 1+.SP, we have 


1+ SP = (a/b)}1 + [1 + 6 — (QSb/a)?]"?} (19) 


where b = (1 —o°)/o*. Remembering Eq. (15), it is seen that in (19) we havea 
relation, free of all circuit and interferometer constants, involving only 
galvanometer readings, values of r, and the values of a and 8. Having de- 
termined the ratio a/8 from a series of observations, using Eq. (17), we may 
assume a value of a, from which we may compute 8 and S’. The corresponding 
value of S is then found by means of Eq. (18). Choosing any value of ¢ on 
the r, 0 curve, preferably at about one-half the altitude of a peak, we com- 
pute the value of b. The corresponding value of 7 is found from the curve for 
use in computing P and Q. These functions involve r both in the hyperbolic 
and circular functions composing them. The argument, 27w/v of the circular 
functions may be found directly if the frequency is known, but may be found 
independently of frequency from the interferometric determination of \/ 2, 
since 6=2rw,/v=4Arz,X, where Ar is the distance from the point for which ¢ 
is chosen to the axis of the peak, or one-half the peak width, since, owing to 
its extreme narrowness, the symmetry of the peak is not affected by damping 
to a determinate degree. It is now possible to compute the left and right hand 
sides of Eq. (19). If the left hand side is larger than the right, the tentative 
value of a which was chosen is too large, etc. Successive trials lead, usually 
in three or four steps, to the region of intersection of the curves of the two 
sides of this equation plotted as functions of a, when by interpolation the 
value of @ satisfying the equation is found. The process may be shortened 
by noting that 1+.SP must lie between the values a/b and (a/b) |1+(1+0)"?] 
and by using the approximations P=4(ra+8)/ |4(ra+8)?+6], and Q 
= 24 |4(ra+8)?+6?| in all but the last two steps. 


Measurement of velocity 


The values of r for successive values of maximum ¢? given by Eq. (11) are’ 
to a degree of precision beyond experimental measurement, independent of 
the damping," except in excessively absorbing media, and are determined by 
the maximum values of P, which in turn occur when cos 27w/v=1, or when 
r=nX/ 2. It then remains to measure the frequency in order to evaluate the 
velocity 7. 

Il]. EXPERIMENTAL TEsTs OF THE THEORY 


Three interferometers have been built for the purpose of testing the 
theory. As the results obtained from all of them differ in no important detail, 
no extensive description of them will be given. 

Interferometer /, Fig. 3a, was designed for operation in the open air. It 
consists of the quartz plate Q resting upon a grounded electrode mounted at 
the top of a Bakelite support. This support is attached to a plate provided 


8 Cf. also C. D. Reid, Phys. Rev. 35, 829 (1930). 
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with adjusting screws so that the face of the plate may be made parallel to 
the plate glass reflector R. This adjustment was made by viewing fringes in 
sodium light passed through R and falling upon Q when the two were close 
together and R was started in retreat from the crystal. The supports for R 
and Q are mounted upon the ways of a micrometer comparator of great ac- 
curacy made by the late Professor H. A. Rowland. The screw is provided with 
a large head divided into 100 parts, and has a fine adjustment making pos- 
sible measurements to 10“ mm. This interferometer was particularly useful 
in testing the effect of various methods of crystal mounting and types of 
electrodes. Its great accuracy and fine adjustment were used both with the 
method described in this paper, the crystal being used as a resonator, and 
with the crystal used as an oscillator, after the methods of Pierce“ in order 
to test the effect of variation of velocity with reflector distance reported by 
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Fig. 3. 


Reid.’ The effect found by Reid was of decreasing amount with increase of 
frequency, being small with the highest frequencies which he employed. As 
the frequencies so far used in this laboratory are much higher the effect to be 
expected might be so small as to escape detection. So far, in the range 400 
to 700 k.c. no such effect has been detected by either the resonator or oscillator 
method, although the interferometer distance has been varied from the face 
of the crystals used to a distance of 10 cm. 

Interferometer II, Fig. 3b, was designed to give the largest possible re- 
action upon the crystal by the fluid column. It consists of a solid cylinder of 
brass in which a concentric hole 1’’ in diameter was bored. The upper end 
of the cylinder was ground plane. Upon it rests the crystal Q, the lower face 
of which is gold sputtered providing the grounded electrode. Contact is made 
with the upper electrode by means of a light spring. The plunger P is close 
fitting, with end R ground plane and parallel to the surface G. P is displaced 
by the micrometer screw S. 

Interferometer III, Fig. 3c, was designed to test the effect of the removal 
of water vapor and CO, from the air and to study other gases at room tem- 
perature. The end of the micrometer screw serves as the reflector, its paral- 
lelism to the surface upon which the crystal rests being tested by sodium 
fringes, with an optical flat in place of the sputtered crystal. The air or other 
gas to be tested enters a small channel at the side of the screw, passing up- 


4G. W. Pierce, Proc. Am. Acad. 60, 269-302 (1925). 
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ward and through a small slot under the crystal and into the chamber above, 
from which it escapes. 

The constant frequency oscillator consists of a dynatron oscillator with 
secondary output connected to the inner grid of a UX222 tube having 45 
volts applied to the outer grid. The output of this tube is amplified in two 
stages of resistance capacity coupling, with nonregenerative tuned output 
in the last stage. Loosely coupled to the output is the coil Z; of the resonant 
circuit which excites the interferometer. Minute variation of the frequency 
of the dynatron oscillator is provided by a precision condenser of large ca- 
pacity in series with the condenser of the dynatron circuit. In addition a 
vernier condenser of very small capacity and with a linear scale of 50 cm 
length is connected in parallel with the precision condenser so that rapid ad- 
justment for symmetry of crevasse may be made. Currents are measured by 
a Western Electric vacuum thermocouple and a galvanometer provided with 
a number of shunts. 

Connection to the interferometer is made through a mercury switch, 
which serves also at any time to connect the interferometer into a Pierce cir- 
cuit so that studies by both the resonator and oscillator method may be 
compared, the crystal being under identical conditions. The results of these 
comparative studies will be presented in another paper. 

The most general conclusion which can be drawn from the experiments so 

TABLE I. Galvanometer and screw readings. 597 k.c. First twelve rows give i», with correspond- 
ing peak number m. Column (1), room air, interferometer II: (2), room air, interferometer III: 
(3), dry air: (4), dry, COs—free air: (5), CO. 








m (1) (2 (3) (4) (5) 
1 0.001484 0.000720 0.000698 0.000703 0.000653 
2 1432 674 642 647 541 
3 1383 634 603 614 484 
4 1344 603 577 583 444 
5 1295 582 557 557 420 
6 1255 563 534 540 394 
7 1211 539 513 514 386 
8 1163 524 498 498 377 
9 1119 508 485 486 

10 1094 504 484 483 

11 1050 49? 467 470 

12 1020 483 462 4506 

rr 22.2 23.3 24.0 24.0 23.3 

A/2 0.02928 0.02944 0.02941 0.02944 0.02389 

10 0.000235 0.000310 0.000278 0.000279 0.000304 

Io 0.002398 0.002329 0.002275 0.002293 0.002275 

recom 0.2343 0.3538 0.3529 0.3533 0.1433 

2 Ar 0.00322 0 .00159* 0.00157 0.00156 0.00200 

1 0.0006 0.000375 0.000350 0.000350 0.000350 

rem 0.3514 

2 Ar 0.00325 


1 0.0006 


* Peak structure not satisfactory owing to rapid drift of frequency and necessity of re- 
peated correction. 
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far carried out is that for every case investigated, in the frequency range 
400-700 k.c. and with all three of the interferometers just described, the 
results for the maxima and minima of the o, 7 curve can be expressed within 
experimental error by the hyperbolic relation given in Eq. (17). The value of 
a/ 8 and of S’ for any set of observations may be determined by the method 
of least squares. Table I gives the original data for five cases at the frequency 
597 k.c. The computed values of a 8 and S’ for each of the sets are given in 
Table II. Fig. 2 is drawn for the results of column (3). Fig. 4, Inserts (a) and 
(b), show the structures respectively of the 12th peak for dry air, and the 6th 
peak for COs. From these and other similar peaks the data for 2Ar in Table I 
are obtained. For the purpose of comparison of different sets of results taken 














| 
¢ . 
i a e? > 
0.1 —-4, a (b) ~ 
.1cm r Loos 
O14) 0.1456 Q1450 


Fig. 4. Reduced observations compared with 1/(ra/3-+1). 
Curve I, room air in interferometer II; Curve I], with circles, dry air in interferometer III; 
Curve III, CO, in interferometer III]. Insert (a), twelfth current peak for dry air; Insert (b), 
sixth current peak for COs. 


with different gases or under different conditions, as with different inter- 
ferometers, we may plot [(1—o0)/o0||o,/(1—@») —o0, (1—o0) |/ S’ asa fune- 
tion of r. This should be compared with its theoretical equivalent, 1/(ra/8 
+1). Fig. 4 shows as an example the first of the above functions for the ob- 
servations of Table I indicated by circles, points, ete., and full curves, I, IT, 
III, for 1/(ra/8+1) for three of the sets corresponding respectively to 
columns (1), (3) and (5) of Table I. It is to be noted that all such curves start 
with ordinate unity at r=0. 


Absorption measurements 


Having determined a/Q, it is possible, as has been explained, to find the 
separate values of a and 8 from some value of ¢ and its corresponding r by 
a method of graphic interpolation. The accuracy with which these quantities 
may be obtained is affected by the constancy with which the frequency and 
amplitude of the driving circuit may be maintained. In the absence of a piezo- 
electrically driven generator of required frequency for driving the resonant 
circuit containing the interferometer, it has been the practise to test the 
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constancy of frequency at the conclusion of readings for each peak by setting 
r at the nearest value of (2n+1)\/4 and noting whether 7 is at the minimum 
value ip. To test whether the amplitude of output has been changed, the 
interferometer is disconnected and a reading of Jp is taken, corresponding to 
the resonance maximum of the interferometer circuit. If either 7 or Jp has 
changed the intervening observations are rejected. This practise has entailed 
a loss of time and labor, making it highly desirable to use crystals for driving 
the generator. However, since the particular crevasse of a given crystal most 
suitable for interferometric work is often not at one of the principal fre- 
quencies of the crystal, the process of accumulation of plates for driving pur- 
poses involves first a careful exploration of crevasses for their suitability in 
the interferometer, and then of preparing oscillator plates, one for each fre- 
quency which can be used. 

The results given in Table I were taken with a rectangular plate. The 
crystal was an excellent oscillator and no other crevasses were observed in the 
neighborhood of its principal crevasse frequency. The peaks were symmetrical 
and showed no evidence at any time of having any satellites. This crystal was 
therefore regarded as being particularly suitable for the present purpose of 
illustration of the theory. The values of a/8, S’, a, and y computed from the 
results given in Table I are shown in Table IT. 


TABLE II. Absorption and reflection coefficients. v=597 k.c. 


Interfer- 


ometer Gas rc \/2 cm a/B 5” a Y 
1 Room air ry Be 0.02928 5.045 16.63 0.130 0.9745 
2 III Room air 23.3 0.02944 6.384 2.262 0 .106* 0.9836 
3 6 Dry air 24.0 0.02941 6.276 2.548 0.110 0.9827 
4 . Dry CO: free air 24.0 0.02944 6.312 2.561 0.108 0.9831 
5 . CO, 23.3 0.02389 94.78 5.705 0.702 0.99261 


* See remark under Table I. 


Data on absorption are very meager. From Neklepajew we have for 
dry air, w=0.00073 /?, which for \/2=0.02941 gives w=0.211. For @ this 
gives 0.105, in fair agreement with the value 0.110 for dry air in Table IT. 
For COs, we have the calculations of Herzfeld and Rice“ from the velocity 
measurements of Pierce“ and by extrapolation from the results of Abello” 
giving respectively for wu, 3X10~'y? and 5.6X10-'v*. We have, for vy =5.97 
«10°, from the former, a=0.54, and from the latter a=1.0. These results 
may be compared with a=0.702 for CO, in Table IT. 


The coefficient of reflection 


An interesting feature of these experiments is the large departure from 
unity of the coefficient of reflection y. If we take as the coefficient of reflection 
of air waves at a surface of brass, y = (Ri— Re)/(Ri+ Re), where R; for brass 


1s N. Neklepajew, Ann. d. Physik 35, 175-181 (1911). 


16 K, F. Herzfeld and F. O. Rice, Phys. Rev. 31, 695 (1928). 
7 T. P. Abello, Proc. Nat. Acad. Sci. 13, 699 (1927). 
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is pv =8.5X3.65 X10°= 3.110%, and for air, Re=40, we expect y to be 
1—1.310% The departure from unity in these experiments is actually 
several orders greater, as may be seen from Table IT. In explanation of this 
we may recall that we are dealing with a resonant system, and that it is to be 
expected that resonance in the reflector system will result in a motion of the 
retlector surface depending upon the dimensions and material of the reflector. 
The coefficient of reflection in terms of motion of reflector surface is con- 
sidered in Part I of this paper. It is to be noted in Table IIT that y changes 
abruptly from interferometer IT to IIT, and from air to COs. Similar changes 
are noted in passing from one frequency to another. It would be of interest 
in this connection to use as a reflector a second piezoelectric plate of quartz 
so that the amplitude of vibration of the reflecting surface could be studied 
directly. The equations developed for this case, where it was assumed that 
the retlector was rigid, I, Eqs. (17) and (18) do not apply in so far as y is not 
unity. From I, Eqs. (12) and (13) we have putting «=r, 


P,= |(A+y)(1—y2e72*) > cos (rw /z) | [1 — 2y2e> 2" cos (Qrw 2) Fyte 4] (20) 
O,=|(1+y)Aty%e 2 er sin (rw v) | [1 —2y2e- 2" cos (Qra 2) fy tral (21) 


These equations, with Eqs. (1) and (2) applied to the reflector, considering its 
network in combination with a suitable thermogalvanometer connected to 
its electrodes, would make practicable a detailed study of retlection phe- 
nomena. Grossmann and Wien'’ have used a piezoelectric plate as reflector 
for the study of frequency variation of a radiating crystal oscillator as a 
function of reflector distance. Grossmann'® has used an arrangement of two 
crystals as source and detector for the study of absorption in COs, but at this 
writing details of his method and procedure are not available. 

It is of interest to compare the effect on the observations of using different 
cross sections of fluid column. From Eq. (12), 4 =SR_ Bpz, and from Eq. 
(13), R=oo/ (1 —a)woi(C+ K,). In changing under the same conditions from 
one interferometer to another, Bpy and w@,(C+A,) remain practically un- 
changed, so that the areas of cross section of the two columns should be in 
proportion to the respective values of Soo/(1—o0), or S’a)3/(1—o0)(1+8/ 2). 
For interferometer II, room air, we have S’=16.63, 8=0.0258, 0) (1—«a») 
= (0.1086, and for interferometer III these quantities have the respective 
values 2.262, 0.0166, 0.1537, implying a ratio of cross section of interferometer 
II to III of 8.03. The actual areas were 5.07 and 0.707 cm®?, having a ratio 
7.18. The agreement of these figures is closer than might be expected in view 
of the remarks following Eq. (16). The particular crystal for which the above 
numbers are cited vibrates vigorously for the mode considered at two points 
near opposite corners of its rectangular area, little trace of motion being 
found over much of its surface. In using such a crystal with interferometer 
III it is necessary to adjust the crystal so that one of its vibrating areas is 
exposed to the gas. The nearness to equality of the ratios given above and 


18 E, Grossmann and M. Wien, Phys. Zeits. 32, 377-378 (1931). 
19 E, Grossmann, Phys. Zeits. 33, 202 (1932). 
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similar values in many other cases point to the conclusion that what we have 
called the effective value of A is nearly equal to the actual area exposed to 
radiation. It seems probable that this equality would be exact with crystals 
cut in the manner described by Straubel*® in which case the absolute value 
of Sin Eq. (15) could be determined and the coefficient of absorption of the 
gas could be evaluated directly from observations of maxima and minima 
of a, using Eqs. (15) and (16), thus avoiding the laborious procedure of 
analysis of peak form. 

A study is in progress, with the methods here outlined, of the ultrasonic 
absorption constants of several gases, considered as functions of frequency 
and temperature. 

My thanks are expressed to the Rumford Committee of the American 
Academy of Arts and Sciences for a grant with which the quartz crystals used 
in this work were obtained, and to Mr. Milton Lipnick, a graduate student 
of the Johns Hopkins University for assistance on many occasions. 


20 H. Straubel, Phys. Zeits. 32, 222 (1931). 
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Mutual Impedance of Grounded Wires above the 
Surface of the Earth 


By RONALD M. Foster 
American Telephone and Telegraph Company 


(Received July 27, 1932) 


FORMULA already established for the mutual impedance of any 

grounded thin wires lying on the surface of the earth' has now been 
extended to include wires lying in horizontal planes above the surface and 
grounded by vertical wires at their four end-points. With the same basic 
assumptions and the same general method of derivation as before, the ex- 
tended formula is found to be: 


d*Q(r, H, h) 
Z2 = ff + cose V(r, H,h) |dS ds, 
dS ds 
where 


Or, H,h) = Qolr) + Qilr, H+ h) — Qor, HW — hy) ), 
V(r, H,h) = No(r) + Nir, H +h) — Nio(r, | H—- hh), 


Oo(r) = p/2rr, 


me ergk tae Gl + Py" — 2) 

O(r, hk) = A : . aia 
todo Va wt Lt $y? + aS 
lwp ( r- + hh?) 1/2 ‘+ h 

Qo(r, h) = h log : 


a 


“a 


— (r+ fh?) 24+ r|. 


r 


p 
No(r) = ——[1 — (1+ Pret], 


2rr 


lwy x (u? =o [?)! 7 ose m 
Ni(r, h) =- | (l-—e | . r - [rotrnda, 
dr e/ i) (u- + lr)! = — M 


iwi 1 1 
Nr, h) = - | — - | 
drLr (r? + h?)'/? 


The integrations in the iterated integral are extended over the two wires 
S and s, lying in planes at heights /7 and h, respectively. The elements dS and 
ds are separated by the horizontal distance r and include the angle € between 
their directions The resistivity of the earth is p; the inductivities of the air 
and of the earth are each assumed to be equal to that of free space, v. The 
propagation constant of plane electromagnetic waves in the earth is , which 


'R. M. Foster, Bull. Amer. Math. Soc. 36, 367-368 (1930); Bell System Tech. J. 10, 408- 
419 (1931). 
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ecuals (iwv/p)'’?. All distances are measured in meters, Z;. in ohms, and p in 
meter-ohms; v has the value 1.256 microhenries per meter; @ is equal to 27 
times the frequency. 

For zero heights, the functions Q and N reduce to Qo and No, which are 
the values previously obtained for wires on the surface. For any heights but 
with low frequencies the formula gives the direct-current mutual impedance 
as found by G. A. Campbell.? For any heights but with one wire of infinite 
length the formula gives the mutual impedance as published independently 
by F. Pollaezek,’ J. R. Carson,‘ and G. Haberland.’* 


24. A. Campbell, Bell System Tech, J. 2, 1-30 (1923). 
3 F, Pollaczek, Elekt. Nachr. Tech. 3, 339-359 (1926). 

4 J. R. Carson, Bell System Tech. J. 5, 539-554 (1926). 
5 (G. Haberland, Zeit. angew. Math. 6, 366-379 (1926), 
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Prompt publication of brief reports of important discoveries in physics may 
be secured by addressing them to this department. Closing dates for this depart- 
ment are, for the first issue of the month, the twenty-eighth of the preceding 
month; for the second issue, the thirteenth of the month. The Board of Editors 
does not hold itself responsible for the opinions expressed by the correspondents. 


The Frequency of Resonance Radiation 


It is of considerable importance in many 
astrophysical studies to know whether the 
atom emits as resonance radiation exactly the 
same frequency which it absorbs. If one sup- 
poses, following Eddington, that it does emit 
this precise frequency, then an atom which 
absorbs light of frequency slightly off “reso- 
nance” will also emit off resonance. A prelimi- 
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nary test was made by illuminating mercury 
vapor with a line (2536A) having a narrow 
and sharp reversal, and observing the reso- 
nance radiation. This was accomplished by 
using a water-cooled quartz are shining 
through a quartz absorption cell onto the 
resonance chamber. The temperature of the 
absorption cell could be regulated electrically 
and with ice, and consequently the vapor 
pressure of the absorbing column of mercury 
vapor could be controlled. With the absorp- 
tion cell cold (0°C) resonance radiation was ob- 
served; and as the temperature of the cell was 


raised the resonance radiation became fainter, 


until at 45°C the resonance could barely be 
detected on a long exposure photograph. This 
absorption cell produces a narrow and sharp 
reversal in the center of the emission line of 
the arc, the amount of the reversal being of 
course determined by the temperature of the 
cell. Yet at no time was any reversal of the 
line in resonance observed. Careful tests were 
made for stray scattered light. The spectro- 
graph used was of the quartz Cornu prism 
type with a dispersion of about 3A per mm at 
2536. 

The physical interpretation of this effect is 
not fully understood. If one supposes a slight 
change in frequency to take place bet ween the 
acts of absorption and emission, then energy 
which is absorbed in the wings would be re- 
radiated in the center, as observed. The dif- 
ference in energy between that of the ab- 
sorbed and reradiated quantum presumably 
appears as thermal energy of motion. It will 
be noted that this change in energy, though 
apparently small, is far greater than the 
extremely minute Doppler width of mercury 
vapor at room temperature could possibly ac- 
count for. If on the other hand we postulate 
detailed balancing and require the exact fre- 
quency absorbed to be reemitted, then the 
emission of a nonreversed line becomes some- 
thing of a puzzle unless the upper state is more 
sharply quantized than the ground state— an 
improbable eventuality. 

It is proposed to continue these experiments 
with a spectrograph of greater dispersion and 
to attempt to resolve the resonance radiated 
line more completely. 


S. A. KorFr 
Carnegie Institution of Washington, 


Mt. Wilson Observatory, 
July 10, 1932, 
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Use of Argon in the Ionization Method of Measuring Cosmic Rays 


In designing apparatus for a cosmic-ray 
survey, we had occasion to study argon under 
pressure as material for the ionization cham- 
ber. The results of preliminary measurements 
showed that this gas was much more strongly 
ionized than either air or nitrogen, which have 
been used hitherto. A study was made of the 
ionization-pressure curves of argon and _ air. 
The source of ionization was about one milli- 
gram of radium at a distance of one meter. 
Even at the highest pressure (95 atmospheres) 
of argon that was used the curve was still 
rapidly rising. 

The following current intensity ratios were 
obtained from our own measurements and 
from measurements of others on air and nitro- 
gen, for a pressure of 95 atmospheres: 


I argon)/ Tair) = 3.22; 
Tnitrogen)/ Lair) =1.48 


I argon) I nitrogen) =2.17 


When the ionization chamber was first filled 
with argon to a high pressure, the ionization 
current increased with time. The value on 
fresh filling was only about 50 percent of the 
final steady value; after 3 hours it was about 
87 percent and the final steady value was 
reached in about a day. This effect depends 
obviously upon how much the gas is stirred 
in filling the bomb. In view of this effect our 
ionization-pressure data have been taken in 
the direction of decreasing pressures. 

When no electric field was applied across 


the ionization chamber, ions produced in 
argon at high pressure under exposure to 
gamma-rays remained uncombined for many 
minutes. This is in accord with the small re- 
combination coefficient that had been found 
for argon at low pressures. This effect was also 
found in air but to a less degree. 

An initial recombination of ions in air and 
nitrogen at high pressure before they are 
separated by diffusion has been used to ac- 
count for the fact that in these gases the ef- 
fective ionization is not proportional to the 
pressure. It is reasonable to suppose that in 
argon the low recombination rate makes im- 
probable such initial recombination, with the 
result that the ionization continues to increase 
with pressure over a much wider range. 

This small speed of recombination in argon 
allows a relatively large concentration of ions 
to form under the action of an ionizing source 
before equilibrium is reached. If now a poten- 
tial is put onto the ionization chamber, the 
initial current is many times the value of the 
steady final current, which is ordinarily used 
to measure the rate of formation of ions. This 
cumulative action suggests the properties of a 
photographic plate, and indicates its possible 
application in the measurement of weak 
ionization currents. 

ARTHUR H,. Compton 
Joun J. HopFietp 
University of Chicago, 
July 18, 1932. 


Strain and Magnetic Orientation 


Previous work' has shown that tension 
establishes the axial as the preferred direction 
of magnetization in a 15 Ni 85 Fe wire, but 
this direction is unique in other respects such 
as symmetry of internal strains and partial 
crystal orientation. More recently further ex- 
periments with large Barkhausen discontinui- 
ties have given striking and direct confirmation 
of the view that a large elastic strain estab- 
lishes a preferred orientation independent of 
other factors. 

The NiFe wire, lying with its axis in the OH 
direction, Fig. 1, was strained by twisting, 
with the result that, referred to the surface of 


1K. S. Sixtus and L. Tonks, Phys. Rev. 37, 
930 (1931). 


the wire, the maximum elongation lay in the 
OP direction. In this state the wire exhibits 
self-propagating changes of magnetization 
under the influence of longitudinal magnetic 
fields as described previously, an important 
feature of the phenomenon being the limiting 
(“critical”) field below which reversal cannot 
proceed. The critical field values for 1, 2, and 
3 turns of one end of an 80 cm long, 0.019 cm 
radius wire are shown in Fig. 1 at Ho, Ho2, Hos, 
respectively. In the present experiments, how- 
ever, a transverse field was added by passing 
current through the wire itself. Only the value 
H, of this field at the surface of the wire will 
be considered. 

In the presence of this field a new set of 
critical longitudinal fields was found. For 
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instance, with J/;,=//,', 77 was increased so 
that the terminus of the vector representing 
the total field traversed RQ from left to right. 
No jump was found for 1 turn; the jump first 


became possible at / for 2 


turns, and for 3 
turns the discontinuity occurred against an 
opposing longitudinal field! 

Confining attention to the 2-turn case, we 


have in Off, and Ollg: two vector critical fields. 
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maximum elongation is solely effective in 
causing the large Barkhausen jump. 
Characteristics were also obtained for a 
wire twisted while under strong tension. For 
zero torsion 77; was without effect, the char- 
acteristic being a straight line parallel to the 
H, axis. For increasing torsion the lines in- 
clined toward the //; axis, the angle of inclina- 
tion agreeing well with that calculated for 
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Fig. 1. Vector critical fields in twisted NiFe wire. 


By using other values of #7; the curve joining 
he and Hq was filled in. This curve is the same 
whether /7; is constant and #/ is increased in 
obtaining it, or vice versa. 

The significant features of the curve are 
that (1) it approximates to a straight line and 
(2) it makes 45° with the axes. The parallelism 
of the 1- and 3-turn curves with the 2-turn 
curve shows that this is no accident. As it is a 
straight line, the components, taken along the 
direction of its normal, of all the vector critical 
fields are equal, and as it is at 45°, this direc- 
tion coincides with OP, the direction of a 
principal strain. It follows that the com- 
ponent of the applied field in the direction of 


the deviation of the principal strain from 
the axis. 

These results lead naturally to the view that 
within the range of magnetic fields used and 
with sufficient strain applied, the magnetic 
intensity lies in the direction of maximum 
elongation, and that transverse fields are in- 
effective in aiding the reversal of magnetiza- 
tion under these conditions. 

Lewt Tonks 
K. J: Saures 

Research Laboratory, 

General Electric Company, 
Schenectady, New York, 
July 29, 1932. 


Electrons as Cosmic Rays 


Professor A. H. Compton's fundamentally 
interesting experiments on the variation of 
cosmic-ray intensity with latitude raises once 
more, as he suggests, the possibility of cosmic 
rays being electrons. In addition to the bearing 


of such a hypothesis upon variations with alti- 
tude and with latitude, certain other matters 
suggest themselves as interesting for con- 
sideration. 

In 1924, in connection with a theory of the 




















LETTERS TO THE EDITOR 541 


absence of ionization! by electrons of suf- 
ficiently high energy, the writer pointed out 
that an electron with velocity only 40 meters 
per second less than the velocity of light, i.e., 
with energy corresponding to 10% volts, could 
not approach theearth’s surfaceat the equator 
nearer than a point at a distance of 8 earth's 
radii from the earth's center when projected 
from infinity towards the earth's surface in the 
magnetic equatorial For the 
favorable direction of projection in that plane 


plane. most 
at infinity, an electron would have to possess 
10'° volts velocity in order to reach the sur- 
face. A. Bramley and the writer? having given 
evidence from the wave mechanics supporting 
the conclusion that electrons of energy com- 
parable with this amount may be incapable of 
ionizing, as formerly surmised by the writer! 


°o. 





more closely comparable with that character- 
istic of photons. 

The actual experimental observation of 
secondary electrons of energies comparable 
with 10° volts, by Professor Millikan and his 
collaborators, raises many interesting ques- 
tions. Such electrons would describe, in the 
earth’s magnetic field, circles or helices of 
radius 107 cm; and their energy, on the basis 
of 16 volts per ion and 50 ions per centimeter 
of path at atmospheric pressure, would enable 
them to travel distances comparable with the 
whole thickness of the homogeneous atmos- 
phere, and indeed distances forming at least 
appreciable fractions of the whole circum- 
ference of the circular orbits into which they 
are bent. The conclusions following from the 
foregoing considerations may best be illus- 


Fig. 1. 


for energies of the order 10° volts on the basis 
of the classical electrodynamics. Such a con- 
clusion would render more reasonable an 
explanation of the maintenance of the earth's 
charge by the influx of such electrons, since it 
would relieve us of the difficulty following 
from the large ionization which the 1500 elec- 
trons per square centimeter per second, neces- 
sary to account for the earth’s charge, would 
otherwise cause. Another point of interest in- 
volves the idea that for such electrons, con- 
sidered as primary cosmic radiation, the 
action of ionization would have to occur 
through the intervention of secondaries shot 
out only occasionally by the primaries. The 


operation of the primaries would thus become 


1\W. F.G. Swann, Phil. Mag. [6] 47, 306- 
319 (1924). Very exhaustive treatments of the 
orbits of electrons in the earth's magnetic field 
have been given by C. Stormer on the theo- 
retical side and by K. Birkeland on the experi- 
mental side, in papers extending as far back 
as 1904. These have to do chiefly with the 
positions of the auroral zones. 

2W. F. G. Swann and A. Bramley, Phys. 
Rev. 41, 393 (1932). Read before the New 
Haven Meeting of the Am. Phys. Soc. on 
June 23, 1932. 


trated in limited space by reference to Fig. 1, 
where in each case the circle represents the 
complete orbit of the electron if there were no 
absorption in earth or air (for simplicity, we 
consider orbits in the magnetic equatorial 
plane). The dotted line represents the earth's 
surface, and the point O, the point of emission 
of a secondary in the atmosphere. Even for the 
case of symmetrical emission of secondaries by 
primaries in all directions in the lower hemi- 
sphere, there will not be symmetrical emission 
of secondaries with respect to a plane passing 
through the magnetic poles. Thus, consider a 
case where the electron path (shown by the 
thick line) is one eighth of the circumference. 
Corresponding to an electron striking the 
earth at 45° from east to west as in Fig. 14, 
there will be an electron striking at 45° from 
west to east as in B. But corresponding to an 
electron striking tangentially from east to 
west as at C, there will be no electron striking 
tangentially from west to east. Such an elec- 
tron would have to have come from below the 
surface of the earth as at D, and to have been 
emitted partially upwards. A lack of sym- 
metry in favor of vertical emission of second- 
aries would enhance still further the consid 
erations here cited. 
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Since the orbits of these high energy elec- 


trons are determined entirely by their rela- 
tivistic mass, and are otherwise independent 
of the velocity since that velocity is practi- 
cally equal to that of light, and since the mass 
is proportional to the energy for these elec- 
trons, the orbits will be the same for protons 
as for electrons of the same energy, except 
that they will be described in the opposite 
sense. 


It is of interest to observe that in the case 
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of observations on a high mountain, radiations 
coming from electrons which have passed 
their lowest point should be directed slightly 
upwards. It is also of interest to note the 
shielding effect to be expected from a moun- 
tain M for such conditions as are represented 
in Fig. 1 E. 


W. F. G. Swann 


Bartol Research Foundation of 
The Franklin Institute, 
August 1, 1932. 
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The Partial Differential Equations of Mathematical Physics. H. BareEMAN. Pp. 522+-xxii. 
Figs. 29. Cambridge University Press, 1932. Price. $10.50. 

To quote from the preface, “in this book the analysis has been developed chiefly with the 
aim of obtaining exact analytical expressions for the solution of the boundary value problems of 
mathematical physics.” The volume is hence a valuable and comprehensive compendium of the 
mathematical technique which has been developed for solving the equations of classical physics. 
Professor Bateman is preeminently qualified to write a monograph of this character, since his 
ability is well known in finding new closed solutions of complicated equations. Among the sub- 
jects treated are Sommerfeld’s method of handling diffraction problems, conformal representa- 
tion, use of Green's functions, toroidal coordinates, etc. In practically all cases, physical prob- 
lems are cited to which the mathematical analysis is applicable. For instance, the first section 
of the chapter on non-linear equations bears the caption “Riccati’s equation, motion of a 
resisting medium, fall of an aeroplane, bimolecular chemical reactions, lines of force of a 
moving electric pole, examples.” This is typical of the diversity of physical illustrations, al- 
though no attempt is made to draw upon the quantum mechanics of atomic structure. As the 
author intimates in the introduction, some of the differential equations and functions which are 
treated are, nevertheless, encountered in quantum mechanics; i.e., the confluent hypergeo- 
metric function. The great generality of the mathematical methods permits a high degree of 
comprehensiveness, but, coupled with the varied and scattered character of the physical 
illustrations, perhaps makes the style somewhat difficult in places. The author index of over 
four hundred names shows careful documentation of modern mathematical investigations on 
boundary value problems. The volume as a whole is focused upon exact solutions, but contains 
an introduction of eight pages on methods of approximation for the solution of differential 
equations, such as the Ritz method, least squares, etc. 

J. H. VAN VLECK 
University of Wisconsin 


Die Relativitaétstheorie. Lupwic Horr. (Volume XIV of Verstindliche Wissenschaft.) 
Pp. 148 +-viii. Figs. 30. Springer, Berlin, 1931. Price RM 4.80. 

This first physical volume to apear in Springer's series of semi-popular scientific expositions 
shows to what extent even such an abstruse discipline as the theory of relativity can be pre- 
sented in a form readily understandable by the educated layman. An illuminating survey of the 
fundamental viewpoints of classical mechanics and electrodynamics is followed by a popular 
exposition of the special theory of relativity which is, in the opinion of the reviewer and within 
the scope of his experience in this field, the best of its kind; the only mathematics involved 
(elementary algebra and geometry!) is contained in several scattered paragraphs of fine print, 
one of which begins with the explanatory sentence ‘‘Auch hier sei eine kurze Rechnung fiir den 
Leser, der den Pythoagoriiischen Lehrsatz und den Begriff der Wurzel noch Kennt, einge- 
schaltet.” The remainder of the book is devoted to the more general theory, and although it 
falls somewhat behind the previous portions in clarity and fitness, is perhaps as successful an 
attempt as could be expected under the severe restrictions imposed. 

H. P. ROBERTSON 
Princeton University 


Die gruppentheoretische Methode in der Quantenmechanik. B. L. vAN DER WAERDEN. 
(Volume XXXVI of Die Grundlehren der mathemitschen Wissenschaften in Einzeldarstel- 
lungen.). Pp. 157+-viii. Figs. 7. Springer, Berlin. 1932. Price RM 9.-RM 9.90. 

The appearance of van der Waerden's excellent exposition of the applications of group 
theory in quantum mechanics will be welcomed by those who wish to survey the field without 
subjecting themselves to the more strenuous task of digesting the more extended treatises of 
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Weyl or Wigner. This little book is characterized by a directness of attack which offers the 
interested reader a shortcut to those portions of the mathematical theory which are of im- 
mediate application in the theory of atomic and (diatomic) molecular spectra. 

A short introductory account of the Schrédinger formulation of quantum mechanics is 
followed by a brief exposition of those theorems of the theory of groups and their representa- 
tions which are to be employed in the sequel. A shortcut to the important theorems concerning 
the reduction of a representation into its irreducible constituents is obtained by restricting the 


proof of preliminary lemmas to the case of an “additive (Abelian) group with operators.” These 
sections are followed by an explicit derivation of the complete set of irreducible representations 
of the rotation and Lorentz groups, and of the selection and intensity rules resulting from the 
former. The Lorentz group and its representations are discussed with the aid of the so-called 
spinor analysis, in a form immediately applicable in the following (fourth) chapter of the 
Pauli-Darwin and Dirac theories of the spinning electron. After a discussion of the relativistic 
invariance and the solution of Dirac’s equations, the author drops back to the nonrelativistic 
Pauli-Darwin theory for a derivation of the Hund theory of atomic spectra — insofar as it does 
not involve the Pauli exclusion principle. This exclusion principle, which is of course bound up 
with the theory of the representations of the symmetric permutation group, is dealt with in 
Chapter V from the Dirac-Slater points of view which enable one to avoid much of the general 
theory of this group; the Hund vector model is then completed by the incorporation of these 
considerations into the scheme derived in the previous chapter. The sixth and last chapter is 
devoted to a brief development of the corresponding theory for the spectra of diatomic mole- 
cules. 

The scope of the present book transcends that of Wigner’s by the inclusion of some twenty 
pages on molecular spectra, and has therefore to obtain its end by the sacrifice of much of the 
illuminating detail contained in the latter. On the other hand, it is narrower in scope than 
that of Weyl—to which it bears a greater resemblance in notation, method and choice of 
material—but insofar as it does deal with the same topics it does so by presenting only that 
portion of the theory which is absolutely necessary for the applications under consideration. 
While it may therefore send the otherwise uninitiated reader to other treatises for a more 
thorough background, van der Waerden’s logically self-contained and on the whole praise- 
worthy account of the applications of group theory to quantum mechanics should be of con- 
siderable assistance to the physicist or chemist who wishes to gain an insight into the field for 


himself—and who is not wholly untutored in the ways of mathematicians. 
H. P. RoBERTSON 
Princeton University 
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